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ABSTRACT 


A new approximate method is developed for calculat- 
ing the temperature distribution in an atmosphere where 
local mechanical -energy dissipation is balanced by local 
net radiative loss in one or more spectral lines for which 
the atmosphere is optically thick. The method uncouples 
the equations of radiative transfer for the individual spec- 
tral lines from the energy-conservation equation and from 
each other. This permits a solution of each individual 
transfer equation for the radiation field in each line. 
These results, used in the energy -conservation equation, 
yield the temperature distribution if the optical depth 
ratio of the spectral lines at each point is known. The 
method is applied to selected strong resonance lines in 
the solar chromosphere. 
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Chapter I 

INTRODUCTION 

A fundamental problem in solar physics is to understand the origin of the high temperatures 
in the outermost region of the solar atmosphere. Assuming only radiative transfer of energy in 
this region, we should expect no temperature minimum above which temperature increases with 
height. In fact, however, Edlen (1943) has shown that the outermost region, the corona, has a 
temperature on the order of 10® °K, which is at least two orders of magnitude above the minimum 
temperature. 

In seeking sources of energy to explain the temperature rise, we can look outside the solar at- 
mosphere as well as beneath it. We shall see that no satisfactory sources have been found outside 
the sun. Thomas (1948a) suggested that spicules might heat the corona if they were supersonic 
jets, but current estimates of ambient gas temperatures for most of the region through which spic- 
ules move require that spicules be subsonic. 

The only mechanical source that has proved promising to date is granulation,as first suggested 
by Biermann (1946, 1948) and, independently, by Schwarzschild (1948). They suggested that the 
granules are a manifestation of turbulent convection; this could generate acoustic waves that would 
build up into shock waves as they propagate out into the chromosphere and corona. This trans- 
formation of acoustic waves into shock waves was deemed reasonable because there is a rapid 
decrease of mass density with height in the chromosphere-corona region, where they proposed the 
shock dissipation should take place. Rough estimates of the energy available for shock dissipa- 
tion do agree in order of magnitude with rough estimates of the radiative energy loss from the 
chromosphere-corona, although (as we shall see) estimates of the available mechanical energy 
and of the radiation losses could both be in error by as much as an order of magnitude. 

Much work has since been done on the problem of what wave modes are generated by the con- 
vection zone, how they propagate, and by what mechanism their energy is dissipated into the 
atmosphere. We shall review this work in Chapter II and, after considering various wave modes 
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and dissipation mechanisms, we shall argue in favor of shock dissipation of acoustic waves as the 
most likely heating mechanism in the chromosphere {to which the calculations of this report are 
restricted). 

Having decided upon the most likely dissipation mechanism, we then consider by what means 
the dissipated energy is redistributed through the atmosphere and eventually lost from the atmos- 
phere in radiation. For most of the chromospheric region studied in this report, conduction proves 
to be negligible. Thus, we shall treat only shock-wave heating and radiative redistribution and 
cooling. By ignoring the transient phenomena associated with the passage of a shock wave through 
a region and the subsequent radiative relaxation, and instead working with time-averaged quantities 
for an atmosphere in a statistically steady state, we can use the following local energy-balance 
condition: 


Ein = Eout ' (1) 

where (in line with convention) E^^is the dissipated mechanical energy, and E^^ is the net radia- 
tive loss. The problem is to derive functional forms for E;^ and E^^^ . If this can be done, the 
temperature distribution in the atmosphere can then be calculated from Equation 1. The main ob- 
ject of this entire study is to exhibit a new approximate method for determining the functional form 
of E^^j in a mechanically heated (i.e, , nonradiative equilibrium) atmosphere, so that the tempera- 
ture distribution in the atmosphere can be calculated. 

The term E^^ can be obtained from shock- wave theory for the calculations reported here. The 
term E^^^ is a function of the radiation fields and associated source functions for various spectral 
lines and continua. In general, these quantities can be derived only from equations of transfer for 
the various lines and continua. 

What we need is an approximate method to imcouple the equations of radiative transfer from 
Equation 1 and from each other, so that we can solve the separate transfer equations independently 
and use the results in the equation, which can then be solved at various heights in the atmosphere 
to give the temperature distribution. Such an approximate, iterative method is given in Chapter IH. 
The so-called non-LTE or general approach to radiative-transfer problems provides the basis for 
everything we shall do. Two relevant references which will be used throughout are Physics of the 
Solar Chromosphere by Thomas and Athay (1961) and Some Aspects of Non- Equilibrium Thermo- 
dynamics in the Presence of a Radiation Field by Thomas (1965). Hereafter, we shall refer to the 
above by PSC and NETPRF, respectively. All of the non-LTE work used in the report will be de- 
rived in detail, following the notation of NETPRF. 
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Chapter II 

THE SOLAR ATMOSPHERE 


A. The Classical Problem of a Stellar Atmosphere 

To construct a model for any stellar atmosphere, including the sun, we must make some 
physical assumptions in order to obtain a number of equations equal to the number of unknown 
quantities. More precisely, we must first specify what quantities we need to describe the at- 
mospheric gaseous ensemble, and then introduce the number of physical assumptions that will 
let us derive these quantities from observational data. We call that set of parameters required 
to describe the physical state of the gas a "complete set.” By physical state of the gas, we mean 
the populations of all internal atomic or molecular energy levels and all kinetic degrees of free- 
dom for all individual particles, as well as the photon distribution, in terms of energy and direct- 
tion. These physical-state parameters must be known at every point in the gas. Our problem is 
to derive this complete set of parameters from the observational data, which we call the "observa- 
tional set." Since the observational set will not, in general, be sufficient to give the complete set, 
we must introduce physical assumptions to supplement our observations. We discuss several 
physical assumptions which astronomers often use, and evaluate these assumptions for the solar 
atmosphere. 

Assumptions commonly used in constructing atmospheric models are local thermodynamic 
equilibrium, radiative equilibrium, and hydrostatic equilibrium. Hereafter we shall refer to them 
as LTE, RE, and HE, respectively. We define them and discuss their applicability. 

LTE, as we shall use it here, means that in each small volume element of gas all the well 
known thermodynamic-equilibrium distribution functions are valid at some local temperature. In 
particular, the Maxwell, Boltzmann, and Saha equations describe the energy state of the matter 
gas; the Planck function describes the energy state of the photon gas. Local temperature is the 
only parameter required to describe the energy state of a small gas element; a complete set of 
parameters to describe the physical state would be temperature, density, and chemical composition. 

LTE differs from thermodynamic equilibrium TE in that temperature gradients are permitted 
in LTE, as long as the energy fluxes resulting from the gradients are too small to disturb the 
local energy distributions given by the TE distribution functions evaluated at the local tempera- 
ture. Obviously, a star cannot be a strict thermodynamic equilibrium enclosure, since a TE en- 
closure is defined as a system isolated from its surroundings by an adiabatic wall through which 
no fluxes are permitted, and a star must lose energy to surrounding space or we could not see it. 
Thus, the LTE assumption has been proposed as the simplest assumption we could use (subject to 
verification) to describe conditions in part of the atmosphere. 
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There is one region of any atmosphere where LTE cannot be valid, i. e., the outer region, 
which itself produces negligible radiation but through which passes all the radiation that escapes 
the star. In LTE as defined, the radiation field must be isotropic— a condition obviously violated 
in the above example. To provide a less restrictive assiunption (which we shall see is usually 
still to restrictive) for handling these regions where the radiation field is anisotropic, some as- 
tronomers have introduced the assumption LTE-R. This means that the TE distribution functions 
for the matter are assumed to apply at some local temperature, but that the radiation field is not 
assumed to be either isotropic or equal to a Planck function at that temperature. We must obtain 
the radiation field by solving the equation of radiative transfer for intensity l^, where the restric- 
tion of LTE-R on the matter distribution functions requires that 

s„ = K (Te) , lte-r . (2) 

where is the source function, and (t^) is the Planck function, being the local value of the 
electron kinetic temperature, which is the temperature we shall use throughout this report. 

Note: Not all authors mean the same thing by LTE. For example, Aller (1963) uses it to de- 
note what we mean by LTE-R. To avoid confusion, the reader is urged to keep in mind the defini- 
tions in this section. 

Radiative equilibrium, RE, is the second commonly used assumption we consider here; this is 
defined by the statement that the net integrated flux 

F - [ 

Jo 

remains constant with depth in the atmosphere, F^, being the monochromatic flux given by 

^ 2tt [ l^(/j.)/xdM, (3) 


where m = cos5 and 0 is the angle between an outward normal to the atmosphere's surface and the 
direction of the radiation intensity Assumption RE obviously implies that 

f = 0. (4) 

where h is the depth parameter. It an atmosphere is "in RE" then, in the absence of local, non- 
radiative sources of energy, the radiation field carries a constant fraction of the total energy trans- 
ported through the atmosphere. In the absence of convective or conductive energy transfer, this in 
thrn means that the radiation field carries all the energy transported through the atmosphere. 
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Hydrostatic equilibrium, HE, is the third commonly used assumption we consider. This means 
that all velocity fields are thermal, unless the pressure, P, used in the following equation for HE 
includes the turbulent pressure caused by small-scale isotropic turbulence, in which case this too 
can be treated under the HE assumption. Any large-scale mass motions and all motions resulting 
in a net mass transport across any arbitrarily oriented surface in the atmosphere represent depar- 
tures from HE. Where HE applies, we can use the equation 

dP /,-x 

dh ^ 

where P is the sum of the gas and radiation pressure at a point, p is the mass density, and g is the 
acceleration caused by gravity. 

Consider, in turn, the applicability of LTE-R, RE, and HE to the solar atmosphere. Insofar as 
these assumptions apply, they simplify our problem of deriving the temperature distribution. Where 
they do not apply, we must introduce less restrictive assumptions that can be justified. 

Before evaluating these assumptions, let us clarify what we mean by photosphere, chromo- 
sphere, and corona— subregions of the solar atmosphere we shall refer to constantly. We define 
the photosphere as that region of the sun we can observe anywhere in the visible spectrum lying 

below the following zero point for height, measured positive outward in a radial direction: h = 0 

0 

is that depth in the solar atmosphere at which the radial optical depth at ^ = 5000 A is 1.0 x 10"^. 
This value corresponds to the zero point on the chromospheric height scale used in Chapter 6 of 
PSC. The chromosphere and corona, defined in terms of the radiation characteristic of each, must 
necessarily lie above this h = 0 surface. 

We are now ready to evaluate the assumptions of LTE-R, RE, and HE. The first question we 
should pose is: At what depth in the solar atmosphere does the source function equal the Planck 
function so that LTE-R is valid? In NETPRF it is shown that for "collision-dominated" lines an 
approximate measure of the optical depth above which the source function differs from the Planck 
function is e" * , where € can be considered the ratio of collisional deexcitation to radiative deexci- 
tation (€ is defined more precisely by Equation 27). Thus, for those lines for which e « 1, the 
optical depth above which (t^) will be large compared with unity. We shall see that e -10“ 

is a resonable value for several strong spectral lines treated in this report. So the optical depth 
at frequency corresponding to the center of these lines above which LTE-R is invalid, is 

10.“. For these lines, only for 10“ can we be sure that LTE-R applies. On the other 

hand, we shall see that (t^) is a good approximation for the solar continuum source func- 

tion in the photosphere. Each spectral frequency regime must be tested independently for the va- 
lidity of LTE-R. Where LTE-R is invalid, a more general expression for the source function than 
S^, = B„ (t^ ) must be used. 

Assumption RE provides a photospheric temperature distribution close to the empirically de- 
termined one we adopt, but RE is certainly violated in the chromosphere and corona where the 
presence of nonradiative energy sources is necessary to maintain the high temperatures. Thus, 
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Equation 4 is not valid for these regions. What we need for the temperature calculation there is 
an analogue of Equation 4 for the case of departure from RE caused by the nonradiative heating 
mechanism. To derive this, we start with the trai^fer equation for a plane parallel atmosphere, 

dl. 

^ = I. - S. . (6) 

with the usual notation, as given in Kourganoff (1963). Source function is a composite of line 
source function and continuum source function S^. Introducing ; 

d-7-c 

= dTl ' (7) 


where the subscripts c and I refer to continuum and line, we can write composite source function 
as 




( 8 ) 


Substituting Equation 8 in Equation 6, multiplying both sides by dr^^/dh, operating on both sides with 


II 


■ • • djj. dv , 


and recalling the definition of monochromatic flux. Equation 3, we obtain 


dh 


4t7 




&T9 . 


dTc 

(j. - Se) ^ 


(9) 


where the sum extends over all N lines which are important energy sinks in the region where we 
apply Equation 9. We shall transform equation 9 into a more convenient form for our calculations 
after discussing non-LTE theory. 

We assume HE is our calculations, which extend from h = 0 to h - 1250 km. In the region lying 
above h 1250 km, this assumption is invalid, owing to spicules. It is possible that HE is never a 
good assumption for any region of the solar atmosphere. This question is complex and we shall de- 
fer a more detailed discussion to section E. 

Our general procedure, then, will be to develop an approximate method for calculating the 
chromospheric temperature distribution for a gas in which HE is assumed but departures from 
RE and, when important, LTE-R are taken into account. 
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Another assumption we shall make is spherical symmetry. Obviously, this assumption is in- 
valid in the upper chromosphere where the spicules are observed. We must consider atmospheric 
inhomogeneities in greater detail in the photosphere and lower chromosphere to see if we can as- 
sume spherical symmetry there. We do this in Section E. Finally, the thickness dR of the atmos- 
pheric region we shall treat is very small compared to the radius R of the Sun; this justifies the 
plane parallel atmosphere approximation, which will be made throughout. 

B. The Photosphere and the Temperature Minimum 

The photosphere has a two-fold importance for our chromospheric energy balance problem. We 
shall see later in this chapter that, in the lower (deeper) region of the photosphere, the gas is con- 
vectively unstable and the associated turbulent motions provide the most probable source of me- 
chanical energy for heating the chromosphere. Also, the bulk of the solar radiation field in the 
continuum originates in the photosphere, and it is this background continu;un that will be important 
for the source functions of certain types of spectral lines. This section is devoted to a brief dis- 
cussion of this continuum radiation field and the associated temperature distribution in the photo- 
sphere. We are particularly interested in the temperature distribution in the outer region of the 
photosphere and in the lowest part of the chromosphere. In this region, the temperature goes 
through a minimum, whose position and value we want to know as accurately as possible. 

The most notable feature of the photosphere is the continuous radiation of the Sun. The Sun's 
high surface gravity produces a rapid density increase in the photosphere. This produces a rapid 
increase in both the volmne emission and absorption coefficients, restricting the emitting region 
to depths above h - -400 km on our scale of heights. We can estimate the total power output of 
the photosphere by measuring the radiant energy from the Sun at a position on Earth. All values 
reported in Minnaert (1953) lie at or below 2.00 cal cm"^ min"', which corresponds to a net emer- 
gent flux of 6.50 X 10'“ erg cm'^ see” ' . We shall see that this is at least three orders of magni- 
tude larger than current estimates of radiant energy losses from the chromosphere and corona. 

Of even more interest to us is the spectral distribution of continuous radiation. Early work by 
Abbot (1920) used light from the entire solar disk and this gave the true monochromatic intensity. 
Later work by Peyturaux (1951) used light only from the center of the disk, thus giving the radiation 
intensity along the normal. We compare the two, using a reduction factor derived by Minnaert (1953) 
who tabulates the results. We have graphed these data in Figure 1, along with black-body radiation 
fields at 6500°Kand 5800°K. The curves demonstrate that the solar continuum radiation approximates 
very well a black-body or thermodynamic equilibrium radiation field at the given temperatures. 


Photographs of high spatial resolution taken in white light reveal a significant intensity fluctua- 
tion, with a scale for the brighter features of the order of 1000 km. These bright features are, of 
course, the granules, which will be discussed in Section E. The values for the emergent intensities 
discussed so far are averages, including radiation from the granules and the material between 
them. 
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The temperature distribution in the photosphere has been obtained theoretically and empir- 
ically. The most elementary theoretical temperature distribution that reasonably fits empirical 
distributions in the solar case is the grey-body distribution, given by 

V = fTetf , (10) 


where q(r) is a complicated function satisfying 1/2 < q < 1.0 and tabulated in Chandrasekhar (1960), 
and (the "effective temperature") is defined by 



Figure 1— Spectral distribution of solar radiation based 
on observational data and comparison with Planck 
function. 



Figure 2 — Comparison between observed and grey-body 
temperature distribution in the photosphere. 


F = . (11) 

where cr is the Stefan- Boltzmann constant. From 
our value for F of 6.50 x 10'° erg cm"^ sec"', 
we obtain = 5825°K, very close to the value 
providing a good fit between the Planck function 
and the observed intensity in Figure 1. The 
assumptions underlying grey-body distribution 
are LTE-R, radiative equilibrium, and a grey 
absorption coefficient (i.e., one independent of 
frequency). Figure 2 shows how the grey tempera- 
ture distribution compares with values from 
several empirical studies. 

The closeness of the fit suggests the as- 
smnptions underlying simple grey-body distribu- 
tion are good ones over most of the photosphere. 
It is already clear that gives a good 

representation of the observed solar intensity, so 
it is not surprising that LTE-R should be a good 
assumption. As for RE, it is shown in Section E 
that the deeper layers of the photosphere around 
Tjggg = 1 are convectively unstable, but that above 
this region RE should apply. The grey-atmos- 
phere assumption also proves reasonable be- 
cause the H" negative ion of hydrogen turns out to 
be the dominant source of continuous absorption; 
moreover, studies of the H" absorption coeffi- 
cient for both bound-free and free-free transi- 
tions show that has only a very weak frequency 
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dependence over a frequency band that includes the strongest part of the solar continuum. This 
dominance of H" was first suggested by Wildt (1939); the ion was first observed in the laboratory 
by Branscomb and Smith (1955); a recent reference for calculations of the absorption coefficient 
(H" ), which contains reference to previous calculations, is Doughty et al. (1966a, 1966b). 

CXir discussion of the grey-atmosphere temperature distribution has provided some useful in- 
sights into physical conditions in the photosphere, but has been of no help to us in fixing the tem- 
perature minimum mentioned before. We can show from Equation 10 that decreases monoton- 
ically to a "boundary" value of (t = 0) = 0.8544 , so we learn nothing about the solar atmos- 

phere for T « 1.0 from Equation 10. Can we, perhaps, learn more about the temperature minimum 
from models based in part on observational data? PSC contains a discussion of two temperature 
models based on limb -darkening measurements: the Minnaert and the Pierce and Waddell models, 
which appear in Figure 2. The Pierce and Waddell (1961) model appears the more satisfactory of 
the two, based on Thomas' criticism in PSC. We add for reference the Heintze (1965) model, which 
is a slight modification of the Utrecht Reference Photosphere, and the Bilderberg Continuum Atmos- 
phere, which is discussed by Gingerich and de Jager (1968). The temperature curves from these 
latter two models lie on either side of the Pierce and Waddell curve, which we adopt for our dis- 
cussion here. 

Pierce (1954) concluded that limb -darkening measurements in the continuum for ^ < 0.2 are 
not very reliable for temperature determinations. Accordingly, the Pierce and Waddell limb dark- 
ening measurements end at about this value of m. However Dunn (1959) analyzed observations near 
the solar limb in the continuum at X6600 A, carrying his limb darkening measurements I„ (0 ,m)/ 

(0,1) right out to M = 0. His data for m 5 0.2 are more reliable than those of previous investi- 
gators because he uses a narrow-band interference-filter to reduce the scattered light, thus re- 
ducing the uncertainty Ai^ (0 ,m) in (0 ,m). Dunn's results out to m - 0.03 fall very close to an 
extension of the Pierce and Waddell data curve. Thomas in PSC concludes from these two studies 
that we can have the following (''■5000) ^sooo'^ 0 . 2 : 


^50 00 

0.1 

0.05 

~ 0.02 

Te 

5200 ± 50° 

4900 ± 100° 

-4750 


We accept these values out to = 0.02. 

We now wish to proceed outward to those regions of the atmosphere where limb-darkening data 
are of no avail. We begin by discussing the line-blanketing effect, to estimate roughly the gross 
effects of the Fraunhofer lines on the temperature distribution in this region where ^5000 0.02. 

We ask for the effect of line blanketing on the so-called "boundary" temperature. 

The gross effects on the boundary temperature (or temperature minimum, for our purpose) 
give us our first glimpse of the importance of non-LTE effects on the solar atmosphere. First, 
we recall that for LXE-R, (t^). We now anticipate a result of the next section, and state 



that for many spectral lines the source function can be written as 




1 


+ (tJ 


1 + e 


( 12 ) 


where e is effectively the ratio of collisional to radiative de- excitation. In NETPRF, Thomas 
shows that, if we write the volume absorption coefficient as 

K = K, + Ke + ■ (13) 


where K^, K^, cr^ are absorption coefficients for true continuous absorption, line absorption, and 
scattering, respectively, we can rewrite the transfer equation in a form that yields, (when we oper- 
ate on both sides with / • • • d/Ltdi/ and apply the condition of radiative equilibrium of Equation 4) the 
following two equations for : 


. ^ N 

jK^B„(Tjdv = jK,J„d.-2^K. 

Bvj (Te)- 

L _ 

for the case where is given by Equation 2, and 




c c ^ 

Jk,B, (Tjd^ = J 

j 


(Te)- L^vd^ 

^ J 

e 

1 + f 


(14) 


(15) 


for the case where S.j^ is given by Equation 12. The summation is over the number of lines in the 
spectrum. Equations 14 and 15 illustrate the difference between the LTE-R case and the non-LTE 
case of Equation 12 for spectral lines. The only difference in form is the inclusion of the factor 
e .y/(l + e . ) which multiplies the second term on the right in the non-LTE case. We have already 
noted that for many lines in the solar atmosphere e ^ « 1. Thus, the effect of these lines in lower- 
ing the temperature is less than in the LTE-R case, and we would expect the temperature minimum 
to have a higher value than if we assumed (t^) as in LTE-R. The larger value of the first 

term on the right in the non-LTE case, due to less absorption by the lines (physically, the scatter- 
ing term fj^ dv predominates over (t^) in Equation 12), also raises the temperature. How- 
ever, a detailed calculation is required to determine the extent of this effect, because the change in 
the expression (t^) -fj^ 4^^. dismay reduce the difference between the two cases. 

What we need is a reasonable estimate of the temperature. In PSC, Thomas argues for a value 
not lower than 4750°K, occurring in the interval 10'^ 5 Tj < 10"^ , with the possibility of a somewhat 
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lower value at a slightly smaller He bases this conclusion on the previously mentioned work 
of Dunn and his own model for the chromospheric region 0 £ h (km) £ 1300 which we shall describe 
in Section D. 

A slightly lower value for the minimum temperature is obtained by Curtis (1965) from an anal- 
ysis of the sodium D lines. The analysis is based on limb -darkening data on the sodium D lines 
collected earlier by Waddell (1962). Following a suggestion by Jefferies, Curtis assumes that the 
line source function is equal at a given physical depth for the two components Dj and of the 
sodium doublet and estimates the value and location of the minimum temperature. The derived 
temperature distribution reaches a minimum of about 4600°K at = 10"^ and retains this value 

out to ^5000 ~ 10'“'. This result is not inconsistent with the above-mentioned conclusions of Thomas 
and is also very close to the temperature distribution of the previously mentioned Bilderberg Con- 
tinuum Atmosphere in the range 10'^ > ■^5000 - 10 "' • We shall adopt the value = 4600° for our 
temperature at h = 0, which lies in the above optical depth range. Unfortunately we cannot say with 
certainty that we have not overestimated the minimum temperature. For example, Athay and 
Skumanich (1968) conclude from their analysis of the Can K line that a minimum temperature of 
£ 4200°K is required to match their computed K, profiles with observations. We hope that this dis- 
crepancy will be cleared up soon. 

Now, let us discuss the chromosphere and corona to see what physical conditions exist there 
and what temperature models for these regions have so far been derived. It would be advisable to 
set down beforehand the basic ideas and relations of non-LTE theory, which will be used in dis- 
cussing the outer atmosphere and in the derivation of the equations for the chromospheric tempera- 
ture calculation. 


C. Non-LTE Theory for Solar Spectral Lines 

First, let us emphasize that non-LTE does not mean we are forcing the source function to de- 
part from the Planck function. Instead, it means that we assume the more general condition of a 
statistically steady state and derive a source function which, under the special conditions for LTE-R, 
will equal the Planck function. We shall follow NETPRF quite closely, but necessarily be much 
more brief. 

The non-LTE theory grew from an attempt to express the source function for spectral lines 
formed in an atmosphere where simple estimates readily show that collisional processes do not 
control the internal, atomic- energy level populations. This is certainly the case in the atmospheres 
of the Sun and most stars. The radiation field which, therefore, plays an important if not dominant 
role in fixing local energy-level populations must be considered. This radiation field is produced 
in regions of the gas where excitation conditions are usually not the same as conditions where we 
want to know the source function. Through its radiation field, the source function depends on a 
range of excitation conditions— not just those conditions at the point in question. 

To develop a more general expression for the source function under these conditions, we must 
assume a more general physical condition to exist locally than LTR-R. We presuppose a 
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statistically steady state (abbreviated SSS), given by 


dOk 

“dT = 0 , SSS, (16) 

where Ok is the number density of a given chemical species in the k-th energy level. The non-LTE 
theory based on Equation 16 applies at each instant of time only when the time for a transient per- 
turbation is long compared with the time for a SSS to be established. It may be objected that a 
shock-heated gas does not generally satisfy these conditions. However, in this study, we are con- 
cerned only with time-averaged quantities. We ignore how the gas relaxes after a shock wave pas- 
ses through it, and how a chromosphere- corona develops out of an atmosphere in radiative equili- 
brium after a source of mechanical energy dissipation in the atmosphere has been "turned on." 

We proceed to derive the source function for a spectral line in an atmosphere in a SSS, fol- 
lowing NETPRF. Thomas shows that the source function for a transition between an upper-level u 
and a lower-level L can be written in the following general form: 


where and are called "profile coefficients" for spontaneous emission, absorption, and 

stimulated emission, respectively; g„ and are statistical weights for levels u and L, respectively; 
and the other symbols have their usual meaning. The profile coefficients are normalized, so their 
integral over v and solid angle is Air. If </>j, = >/'„ , the above reduces to an expression written more 
compactly as 


s 


V 


^ / 

T (Te 


(18) 


where 




2hy3 / guUi, _ y ' 

7 


(19) 


The advantage of Equation 18 as a form for the source function is that it separates into two fac- 
tors, giving the frequency dependence of s„ on frequency and the occupation numbers, respectively. 
In LTE-R, 

Bu "l r hv 1 
= ""nWe) ' 

so - Tj. or B,, (Tex) ^.s must occxir for the non-LTE theory to be valid. 
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The factor is evaluated by Thomas (1957) for the central core of resonance lines on the 

assumption that chromospheric conditions apply. He concludes that the profile of the absorption 
coefficient is Doppler out to some three Doppler half -widths from line center, and that the ratio 
j equals 4 at the most, over the same Doppler core. Since the number of absorbers for the 
Doppler (Gaussian) line profile varies as exp|- (Ay/Avjj)2| (where is the Doppler frequency half- 
width) Ai^/Azvjj =z 3 corresponds to a variation in the absorption coefficient of e® - 10“* over the Dop- 
pler core of the line. Since our chromospheric temperature calculations of Chapter IV will require 
source functions for three resonance lines (Ly a, and the resonance doublets of CaH and MgH), all 
of which have their cores formed in the chromosphere, we conclude that Thomas' conclusion on 
applies; furthermore, we ignore frequency dependence altogether by estimatii^ 

S. - (T^ J . (20) 

It is easy to derive a general expression for (t^x) ; See Appendix A. Our starting point is 
Equation 16, expressing the statistically steady state. Following Thomas in NETPRF, we introduce 
a "net rate bracket" notation for treating the net rates of various processes. If P - R refers to the 
net rate of process p (Ris the reverse process), then for P (radiative ; u , L) we define the net radia- 
tive bracket (NRB)^^^ as 


P(rad; U, L) - R(rad; L, U) = (NRB)^^^ 


(21) 


Likewise we define the net collisional bracket (NCB)nj as 

P(col; L. U) - R(col; U, L) = C^u (NCB)j^y . (22) 


InEquations 21 and 22 the Einstein coefficient for spontaneous emission, and is the ine- 

lastic collision rate, which we assume for collisional excitation by free electrons is given by 


c 


LU 


n 


e 



^■lU fe (v) vdv , 


''lu = ■ 


(23) 


where is the electron density, is the inelastic-collision cross section for these collisions, 
(v) is the electron- velocity distribution function (which we assrnne is Maxwellian) and is the 
minimum velocity for collisional excitation obtained from 1/2 , where Xlu " Xl ~ Xu 

the difference in ionization potentials between levels L and u. 

Using the net brackets defined by Equations 21 and 22, we can write a completely general 
expression for the equation of statistically steady state (Equation 16). For a level k, this 
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expression becomes 


7", Z] ") J A^k"k"e (NRB)^^ dx ■■■ 

^>k i<x 


" (NKB)k, + Z] 

-e>k J 


Subscript k refers to states in the continuum; n. is the ion density. According to convention, we 
include collisional processes that have an upper level in the continuum in the last term of Equa- 
tion 24. To obtain the desired expression for , we write Equation 24 for levels u and L of 

the transition. We solve these two equations for then substitute this expression for in 

Equation 19, thus arriving at the desired expression for given below. (Refer to Appendix 

A for details.) The result is 


B (t ) 

V \ ex/ 


jj. di. + eB, (tJ +TJB- 


1 + e + 7) 


(25) 


where 


J. 




(M)dM , 


(26) 




(27) 


Bl ^ Fu 


cr A K K 

^UL * L U 


Ih _ 

Bl ^ ^ 


(28) 


7)B' = 


2hv^ Bl Sj Fjj 


:2 guV^^L'^fu 


(29) 


Quantities S^, S 3 , , and F^ have not yet been introduced; Appendix A gives their detailed 

functional forms. Here, to save space, we merely state what they represent physically. These 
four quantities equal zero, if we consider our atom to be made up of only two levels, U and L. In 
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this case, (T^*) becomes 




+ €B^ (tJ 


1 + e 


(30) 


where collisions and the radiation field in the line itself control (Te,). In general, Sj, S 3 , Fj^, 
and Fy express the dependence of B„ (lex) collisional and radiative processes that couple u and 
L to levels other than u andL. In particular, we see from Appendix A that Sj represents terms that 
depopulate L through all levels except u , 83 represents terms that depopulate u through all levels ex- 
ceptL j F^^ represents terms that populate L through all levels except u , and F^ represents terms 
that populate U through all levels except L . 

When we use b^(t^^) in Chapters m and IV, we find that the line source function of Equation 
30 is probably often a good approximation; in all other cases, we consider the dependence of 8 ^, 

S 3 , Fj^, and Fjj on the continuum only. Here, we work with only two level atoms, or with atoms ap- 
proximated by two levels plus a continuum. When these simplified treatments are not justified, 
we must use Equation 25, including as many additional levels as are needed. 

Now it is evident why LTE-R is a poor assumption for strong chromospheric lines of the sort 
we shall treat in the energy-balance problem. In the two-level atom source function (Equation 30), 
Bv(Tex)"* only if / j,, dv - b^ (tJ, or if e >>1 and eB^ (t^)» /j^<#>^dv. We noted in Section A 

that, for the lines of interest here, e 10 The second alternative above is ruled out immediately. 
To see from what depth in the atmosphere the radiation field escapes at the outer "boundary" = 0, 
and hence how far the influence of the boundary penetrates to producing an anistropy in I,,, negating 
Jj„ i B^ (t^), we use the following fact: e' '■ provides a rough measure of the optical thickness 
of a region through which a photon will move, once it is created and before it is destroyed. Above 
optical depth ' the radiation field is anisotropic, as there are not enough sources lying above the 
point in question to provide an inward directed 1 , 7 to equal the outward directed l^. 

Thus, for Tg < 10 for the lines involved in our final calculations, LTE-R is certainly invalid. 

A rough criterion for LTE-R then becomes: » 1. A more rigorous development in NETPRF 

yields the same result, as well as a very lengthy criticism of the misuse of the LTE-R assvimption. 

Earlier in this section, we mentioned that we would want source functions for Lya and the reso- 
nance doublets of Call andMgn. We shall also find Ha an important chromospheric Une, though we 
shall not solve a transfer equation for it. Jefferies and Thomas (1959) offer some interesting 
general results on Ha and ionized metal lines; Morton and Widing (1961) make a detailed study of 
Lya. A summary of the relevant findings of these authors appears in NETPRF. More detailed 
studies of these lines have recently been made, but they do not affect these general comments. 

In an effort to find the proper form of the source function for different types of lines, the 
general source function of Equation 25 can be written in several different ways, each corresponding 
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to different physical mechanisms of line formation. Jefferies and Thomas considered the following 
two possibilities: 




■» 

(tJ 

B. (T, J = 


(31) 


j dv + tjB* 

— 


(32) 


where v and tjB* are given below. Equations 30 and 31 are identical in form and represent the two- 
level atom. We call fj^ the "storage term, "as it represents scattering which neither enhances 
nor depletes the radiation field in the line. We call eB^ (t^) a "source term," along with r]B* in 
Equation 32. We call e in Equation 31 and v in Equation 32 "sink terms." The source terms en- 
hance the radiation field in the line, while the sink terms inhibit the build-up of this radiation field 
by representing direct (e) and indirect (v) ways of populating the lower level at the expense of the 
upper without producing line radiation. Note that Jefferies and Thomas chose source functions 
whose source and sink terms are colli si on- dominated in the case of Equation 31, and dominated by 
radiation fields in the free-bound continua coupled to the upper and lower levels of the line transi- 
tion in Equation 32. We can see this for the latter case, if we take Equations 28 and 29 for v and 
ignore all levels except U and L of the line and the continuum, and consider only radiative 
transitions between the continuum and u and L. We obtain for rj and B*; 


Tj - (atomic parameters) * 


2hi/3 


W p ^ 

/ Yj-Y^ _ 


- i). 


(33) 


- .)■ 


(34) 


where is the dilution factor for the Lyman continuum, Y,, = Xk/*^lrk > is a radiation tem- 

perature for the mean radiation field in the /< - k free-bound continuum such that (x - k) 

- (irk)- The atomic parameters, in general, vary from line to line. 

The important conclusion of the Jefferies and Thomas study is that there is an apparent one- 
to-one correspondence between type of source function, or mechanism of line formation, and type 
of emitter, or class of line-forming ions in the solar atmosphere. Using source functions given by 
Equations 31 and 32 to solve the transfer equation for representative models of the chromosphere, 
they obtain line profiles typical of ionized metals from Equation 31 and line profiles typical of 
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neutral metals and the Balmer lines from Equation 32. The ionized metals' source functions are 
called "collision-dominated," and the neutral metals' and Balmer lines'isource functions are called 
"photoionization-dominated," after the dominant source and sink terms. The line shape of the two 


types appears is Figure 3. 

The Morton and Widing study considers Ly a, 
which is a strong emission line with a central 
reversal 35 percent below the intensity maxima 
on either side. Evaluating e, (t^), v , and B* 
from Equations 27, 33, and 34, Morton and 
Widing find that eB^, (t^) -»77B* and.T?»e. Thus, 
for Lya, which becomes a line of the "mixed 
type" of source function, 

+ £B„ (tJ 

B. = • (35) 




Figure 3 — Line profiles in cores of (A) collision-dominated 
lines and (B) photoionization-dominated lines— asymmetries 
ignored. 


Equations 31, 32, and 35, to the first approximation source functions for three types of solar 
lines, cover all interesting possibilities. Although we need detailed studies of individual lines (for 
which good observed profiles are required) to specify more completely the details of the line- 
formation mechanism, this rough scheme for classifying solar lines will serve for our preliminary 
calculations. 

There are two more points to be made, before leaving the subject of the source function. The 
ionized metal lines and Lya have source functions for which the source terms exhibit a first-order 
dependence on the local electron temperature . This is not true of the Balmer lines or the neu- 
tral metal lines. This will be of some importance to our discussion of the energy balance in Chap- 
ter IV. Also, the -dependence of the first named lines demonstrates the existence of a tempera- 
ture rise in the chromosphere, discussed in detail by Jefferies and Thomas (1959, 1960). 

The final matter we wish to consider is a convenient form for the net brackets defined by 
Equations 21 and 22. Here we follow NETPRF exactly. We can write Equation 21 in the form 

P(rad; U,L) - R(rad; L.U) = B„j^ f B^„ f dv , (36) 


where B^^^ and are the Einstein coefficients for stimulated emission and absorption, respec- 
tively, and the other symbols have been defined. Combining Equations 21 and 36, we obtain 



"lBlu ( 

B„^ ( 


(NRB)ul = 

^ "u ^UL J 


J,. 


-I UL 


( 37 ) 


Now let us assume that - 4 ^^, as we did to derive Equation 19 for (t^,) ; and then, after using 
the relations for the Einstein coefficients, 


2hv3 

-2 ’ B, 


gu Bul 


(38) 


let us combine Equations 19 and 37. Then we obtain 


(NRB)^ 


jj. 


cf) du 




(39) 


To obtain a simple expression for (NCB)j^y, we proceed in a similar way. We can write Equa- 
tion 22 as 


P(col; L. U)-R(col; U, L) = ’ (“^^) 

where the notation is obvious. Equations 40 and 22 give 


(NCB)^„ 


1 



LU 


(41) 


which can be written as 


(NCB)j,„ 


1 


Su ^UL 



(42) 


where we use the Boltzmann equation for the ratio in LTE, the * indicating LTE, and in- 

troduce so-called "departure co-efficients" bj^, defined by the LTE Saha equation 


n 


k 



(43) 


which is equivalent to 


b 


k 



(44) 
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since the quantity in square brackets in Equation 43 is the LTE value of at local values of 
"e > "i > • Also X^u = XLu/'^^e’ Saha equation, U; is the partition function for the ion. 

Other symbols should be self-evident. Now, in LTE, 


« />• * 

'^L ^LU 


« • r* * 
^UL 


c,*. 




LU 


(43a) 


as (NCB)nj = 0 from detailed balance and = 1. Remembering that we assume a Maxwellian 

velocity distribution for electrons, and noting that and then assume their LTE values, we 
obtain 


(NCB)lu 



(44a) 


This is the desired expression for (NCB)^^. 

Equation 39 for the (NRB)^^ appears in the energy- conservation equation we shall solve for 
our final chromospheric temperature distribution. Equation 41 for the (NCB)^y appears in the ex- 
pression we derive for the apportionment of mechanical energy among all spectral lines formed in 
the same region of the atmosphere. 

Now let us consider current models for the chromosphere and corona— especially the chro- 
mosphere. We shall find that the detailed model of the low chromosphere described in Section D 
is based on studies of emission in line-free regions of the continuum. Thus, at first we need not 
use everything developed in this section, but we must still use the general principle of a SSS con- 
tained in Equation 24. We shall find that significant departures from LTE-R do occur in the low 
chromosphere. 


D. Solar Chromosphere and Corona 

In Section A, we defined the top of the photosphere by the surface which lies at normal optical 

O 

depth 1.0 X 10“^ at X5000 A. This defines the base of the chromosphere and the zero point for our 
vertical height scale, measuring positive upward. The chromosphere and corona comprise the 
entire region lying above h = 0; hereafter we refer to this composite region as the "outer solar 
atmosphere." The chromosphere and corona are usually defined by the types of emission charac- 
teristic of each (van de Hulst, 1953). In general, then, we describe as chromospheric the "flash" 
spectrum seen at the limb in eclipse. This spectrum consists of emission lines of hydrogen, neu- 
tral metals, singly ionized metals, neutral helium, and ionized helium. There is also a background 
continuum of the H“ ion and the lower free-bound continua of hydrogen. Much of the continuum is 
scattered radiation. A more detailed consideration of chromospheric lines and the question of how 
much of the continuum emerging from the chromosphere is scattered photospheric radiation will 
appear in Chapter IV. 
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There is some question about the origin of certain lines from ions of ionization stages inter- 
mediate between chromospheric values (neutral, singly ionized metals) and coronal values (same 
metals - Ca, Fe - many times ionized). For ejample, strong lines of CIV, A. 1548, and OVI, ^1038, 
and M032, are seen in emission in UV disk spectra (Detwiler et al., 1961). We shall not attempt 
to identify these lines as either chromospheric or coronal. Their existence does illustrate a vex- 
ing problem, however; the nature of the transition region between the chromosphere and corona. 

The corona is characterized by emission lines of highly ionized metals, superimposed on a 
background of partially polarized continuum radiation. The corona has a mean density several 
orders of magnitude lower than the chromosphere. As a result, the net energy losses from the 
corona, estimated in Chapter IV, fall about two orders of magnitude below our net-loss estimates 
for the chromosphere. Also, the chromosphere's location just in and beyond the temperature mini- 
mum gives it a basic role in the energy-balance problem, since the source of mechanical energy 
seems to lie in the subphotospheric convection zone. We shall restrict our calculations to the 
chromosphere in this report, subject to some uncertainty as to where the boundary between chro- 
mosphere and corona lies. 

A general feature of chromospheric emission lines discussed in PSC is that their intensity 
scale height is much larger than the pressure-density scale height H = (where R=8.31 x 10^ 

ergs°K“ ^ mole'i is the gas constant, and m is the mean molecular weight) for an isothermal atmos- 
phere, evaluated at our minimum temperature = 4600°K. This strongly suggests an outwardly 
increasing temperature and possible departures from LTE-R, both of which we find in the chro- 
mosphere. An increase in temperature will increase the local pres sure -density scale height. The 
exact role of the non-LTE effects will become apparent only after a detailed study of each emitting 
ion. 


The model we shall adopt to begin our calculations for the first 1300 km above our zero point 
is the model derived in chapter 6 of PSC. Table 1 gives the quantities derived for this model. We 
choose this model, based on data obtained at the 1952 Khartoum eclipse expedition, because the 
method of analysis is highly self-consistent. We cannot accept those earlier models, typified by 
the one given by van de Hulst (1953), which neglect all non-LTE effects in their analysis. The 
Utrecht Reference Photosphere (Heintze et al., 1964) of the photosphere and low chromosphere ex- 
tends out only to h - 230 km on our height scale, and the Bilderberg Continuum Atmosphere (Gin- 
gerich and de Jager, 1968) ignores the conflicting results from line-profile analysis. Thus, for our 
calculations based largely on features of the chromospheric line spectrum, the PSC model of the 
region 0 < h < 1300 km seems best. 

The PSC model for 0 fh 11300 km, a region we shall hereafter call the "low chromosphere," is 
obtained from data in the optical continuum. The choice of continuum data is prompted, in part, by 
a desire to treat that spectral region where strong self-absorption plays a negligible role. A thin 
atmosphere analysis, where a SSS assumption in the form of Equation 24 takes account of non-LTE 
effects, is all that is required. No transfer equation need be solved, as all relevant radiation fields 
are known, except for the radiation field in the Lyman continuum discussed below. 
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Table 1 


Model of Low Chromosphere from Chapter 6, PSC, Using Tables 6-9 and 6-10. (Subscripts 1 and 2 

refer to levels in atomic hydrogen.) 


h 

(km) 

T,(“K) 

(cm" 3) 

10“'^ 

(cm"3) 

Wn„- 
(cm" 3 ) 

10‘® Uj 

(cm " 3) 

b, 

T rad 
>^5000 

T. a 
Ly 

0 

5050 

3.78(15) 

6.36 

38.8 



1 X 10"^ 


100 

5250 

1.60(15) 

5.05 

11.5 



4 X 10“* 


200 

5460 

6.99(14) 

4.65 

4.10 



2 X 10“* 


300 

5680 

3.15(14) 

4.58 

1.62 



6 X 10"® 


400 

5910 

1.46(14) 

4.63 

0.670 



3 X 10" 5 


500 

6150 

6.95(13) 

4.72 

0.290 


3.4 

1 X io"5 


750 

6650 

1.21(13) 

4.04 


9.98 

6.1 


lx 10^ 

850 

6840 

6.28(12) 

3.69 


8.38 




950 

7040 

3.19(12) 

3.38 


7.02 




1000 

7150 

2.27(12) 

3.24 


6.44 

12.6 


2 xlO® 

1100 

7400 

1.10(12) 

2.97 


5.42 




1200 

7700 

5.05(11) 

2.75 


4.61 



3 X 10 = 

1250 

7901 

3.19(11) 

2.65 


4.28 



2 X 10® 

1300 

8310 

1.71(11) 

2.60 

. 

. 

- 



7 X 10“' 


The model is based on observed continuum radiation at V3646, which is the head of the Balmer 
continuum, and at ^4700. Rayleigh scattering due to neutral hydrogen is taken into account; also 
Thompson scattering, free-bound and free-free scattering of H" , free-bound transitions in the 
Balmer continuum at x3646 and in the Paschen continuum at ^4700, and free-free transitions of 
hydrogen at ^4700. The observations give two equations in four unknowns: ,T 3 ,nj^, and 0 ^, An 

assumed hydrogen/metals ratio for h < 500 km or the assumption = Op for h > 500 km adds a 
third relation. For the region h > 500, an expression based on the SSS relation (Equation 24) for 
bj , T^) , defined by Equation 43, provides a fourth equation. The bj (n^,T^) relation is derived 
assuming detailed balance in the Lyman continuum, and a correction scheme is introduced to allow 
for some departure from detailed balance. Spherical symmetry is assumed. 

Actual construction of the model is complicated and involves selecting good data points to 
provide boundary conditions to start the integration, as well as invoking other relations to sup- 
plement estimates made in regions where the data are considered unreliable. For example, hy- 
drostatic equilibrixim is assumed, to get the results for h > 750 km. The data below 500 km are 
not reliable, and begin to enter the region where self-absorption is not negligible as we approach 
the limb. Thomas notes that in the range 0 < h < 500 km the values of the model are an upper 
limit, which is -reassuring in view of the great discrepancy between the (h = 0) of this model 
and our adopted minimum temperature of Section B. 
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What are the conditions beyond the low chromosphere? Hiei (1963) concludes that spicules 
extend at least as far down as h = 1500 km, and there is evidence that they may extend down farther 
(House, 1961). The region 1300 < h <1500 km is certainly unknown territory, and the situation is 
little better from h > 1500 km until the corona is reached. The unanswered questions include: 

What is the state of the gas in spicules? 

What is the state of the gas in the interspicular material? 

What fraction of the surface area at a given height h is covered by spicules? 

How far does the corona penetrate into the interspicular material, if it penetrates at all? 

There is no universally accepted answer to any of these questions. 

In Chapter V of PSC, Athay considers some general stability arguments which he uses, along 
with an analysis of emission in the region h > 1500 km, to produce a model for this region. A 
stability criterion is obtained by requiring that an increase in energy input must cause an 
increased energy output along with the associated initial rise in temperature, such that 


dE 


dE. 


dT. 


- dT. 


(45) 


Thus, at any point in the atmosphere, we require an increase in net radiative losses to counter a 
small temperature increase, and vice versa, or the atmosphere is unstable. Consider a point in 
the atmpsphere where has reached a value at which there are no efficient radiators to provide 
the necessary E^^^ to satisfy inequality 45. In the presence of an Ej^, the atmosphere will not be 
stable at this . Obviously, the atmosphere temperature will rise until, through some process 
such as ionization, a new efficient radiator will be produced and condition 45 can be fulfilled again. 
The atmosphere will be stable at the new . The process is presumably repeated at each stage 
of stable equilibrium, as the temperature increases as height increases until the efficient radiator 
becomes inefficient, at which point the atmosphere finds a new stable equilibrium at a higher tem- 
perature. The process stops when E.^ is exhausted, presumably somewhere in the solar corona. 
Beyond that point the process might be thought to reverse as decreases with height, but other 
complications enter and we wish to consider only the upper chromospheric regions. 

An idealization of an atmosphere, viewed from this simple stability viewpoint, would be a 
succession of temperature "plateaus” connected by infinite temperature gradients. In reality, 
conduction arising from these infinite temperature gradients and radiative energy redistribution 
between and within the plateaus should smooth out the temperature curve, but the gross feature of 
very shallow rises connected by steep rises will presumably remain. As Athay argues in PSC, 
there is evidence from spectral data that this is true for the outer solar atmosphere. 

Athay notes that if observed stages of ionization of common elements are plotted against cor- 
responding ionization potentials, there are gaps in which no stages of ionization of certain common 
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elements are observed in the solar atmosphere, even though a corresponding stage of ionization 
should be expected on chemical- composition grounds alone. The inference is clear: the ion we 
expected to observe is not present over a significant height range. The most direct explanation, 
in the light of our comments above, is that very little atmosphere gas is at the temperature at 
which we expect to find the ion (considering non-LTE effects, of course). This temperature lies 
on one of the steep rises of the temperature curve. The reason is that not enough of the ion is pre- 
sent to provide a sufficiently strong radiator to satisfy inequality 45, so instead of forming a plateau 
at (ion), continues to rise. 

On the basis of these considerations, and after studying emission from hydrogen, Hel, and Hell, 
including the Lyman continuum of hydrogen and Hell, Athay concludes that the hydrogen emission 
comes from the spicules, and the helium emission from the interspicular gas, in the region h > 1500 
km. However, conditions at the top of the spicular component are found to blend with those of the in- 
terspicular matter, suggesting a rapid temperature rise near the top of the spicules. 

We shall not digress to show how Athay constructed a model for the chromospheric region 
1500 < h <4000 km. We only note that his results show a strong dependence on the values of elec- 
tron density used in these regions, and the uncertainty in may be as high as an order of mag- 
nitude. Table 2 gives the resulting model of this region 1500 < h <4000 km, which we shall call the 
"upper chromosphere." We see that he calculates the fractional area covered by the cold compo- 
nent (spicules), and also the and values for both components. 

Table 3 gives the fractional areas covered by both cold and hot components, from Table 2. This 
can be used to construct a temperature model for a spherically symmetrical atmosphere whose ex- 
citation conditions lie between those of the cold and hot components. We can use such a rough tem- 
perature distribution to begin our calculations. 

To construct this temperature model, we assume pressure equilibrium between the hot and 
cold components. Thus, for = P,,, we defines 


eh 

Tec 


(46) 


Table 2 


Table 3 


Model of Upper Chromosphere from Chapter 7 of PSC 
where a. ^ = Percent of Surface at h Covered by Cold Component. 


Estimates of a^ Based on Table 2. 


h(km) 

ac 

r , " " 

logic "ec 

logic "eh 

T 

e c 

T.h 


h(km) 

ac 

ah 

1500 

4-20 

11.9-12.4 

11.1-11.4 

7500 

50000 


1500 

20 

80 

2000 

4-8 

11.4-12.4 

11.0-11.3 

9000 

50000 


2000 

8 

92 

3000 

0-4 

11.6-11.9 

10.8-10.8 

10000 

50000 


3000 

4 

96 

4000 

<60 

10.6-11.3 

< 11.0 

50000 

>50000 


4000 

0.1 

99.9 
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where and are particle densities for the two regions. We now form our "weighted tempera- 
ture" in the following manner: 
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(47) 
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where the fractional areas and appear in Table 2, and in Table 4. 

It is appropriate now to consider some of the 
shortcomings of Athay's upper- chromosphere 
model. Athay knew that it predicted much higher 
intensities of radio emission (due to free-free 
electron transitions in a proton field) than are 
actually observed. To remove this discrepancy, 
he made a complicated alternation, in PSC. A 
more radical suggestion was made by Zirin and 
his co-workers (Zirin and Dietz, 1963 ; Zirin and 
Howard, 1966; Zirin, 1966): that the interspicular 
region is coronal gas, and all the helium emission 
hitherto ascribed to this interspicular gas, as well 
as the radio emission and hydrogen emission, actually arises in spicules. While there is some ad- 
ditional evidence for coronal excitation in the upper chromosphere, no quantitative treatment of 
this problem has so far appeared. We are left with no satisfactory model of this region, and ter- 
minate our calculations in Chapter IV at h = 1250 km. 

Regardless of what model for this region finally emerges, electron thermal conduction is cer- 
tain to play an important role somewhere in the transition region between the chromosphere and 
corona. We do not consider a mechanical flux of energy down from the corona due to conduction 
in the spherically symmetric model we construct, as this conduction is demonstrably negligible 
over the greater part of the region covered by our calculations. 

This completes the background material dealing with radiation-matter-interaction aspects of 
the solar atmosphere. Now we shall consider macroscopic velocity fields and possible mechanical 
energy sources for heating the outer atmosphere. Sections E and F are devoted to these subjects. 


Table 4 

Weighted Temperature for Upper Chromosphere 
Based on Table 3 and Equations (46) - (49). 


h(km) 

R 

A 

C 

Ah 

T 

e W 

1500 

6.67 

0.625 

0.375 

23,440 

2000 

5.55 

0.328 

0.672 

36,550 

3000 

5.00 

0.173 

0.827 

42,130 

4000 

>5 

0 

1 

>50,000 
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E. Velocity Fields and Inhomogeneities 


S, 


’ - 'V' 


So far, we have not questioned our assumptions of HE and spherical symmetry in the chromo- 
sphere in the light of ambient velocity fields and inhomogeneities. Also, observed velocity fields may 
provide indirect support for the energy- generation mechanism and the mode of energy transport dis- 
cussed in the next section. The suggestion for acoustic-wave generation, propagation, and shock 
dissipation at least must agree with these fields. This section is devoted to the observational study 
of velocity fields and inhomogeneities in the solar atmosphere. Some current theoretical interpre- 
tations of these observations will be discussed in the next section, together with the actual mechani- 
cal heating of the atmosphere. 

Photospheric granulation is the first phenomenon we consider. Granulation is an aggregate 
of cells which are brighter than their surroundings when seen in continuum radiation and appear 
roughly circular near the center of the solar disk. Since optical depth unity occurs somewhere 
around h = -350 km in the continuum at X = 5000 A, depending on the photospheric model (Minnaert, 
1953), this is the region of the photosphere where we see disk granulation. A detailed study of 
granulation has been made both by Bahng and Schwarzschild (1962) and by Edmonds (1961). Both 
studies are based on observations obtained from a balloon-borne telescope launched under Sch- 
warzschild's direction (Project Stratoscope). In spite of differences in their interpretations, both 
studies agree that disk granulation represents elements with diameters ranging from a few hundred 
to 1800 km, a mean lifetime of six to nine minutes, and a mean temperature excess of several hun- 
dred degrees above the intergranular matter. 

The best current account of the nature of disk granules is given by Bdhm-Vitense (1958) and 
Biermann, et al., (1959); their quantitative calculations support what had already been conjectured: 
that disk granules are a manifestation of convection at the top of the Sun's subphotospheric convec- 
tion zone. Indeed, they are thought to be the convective elements in the observable region of the 
convection zone. That the innermost region of the observable photosphere is convectively unstable 
has long been known, from applying the familiar Schwarzschild stability criterion to the region's 
temperature gradient. A simple application of Prandtl mixing length theory of turbulent convec- 
tion to this convectively unstable region permits us to examine granulation data to see if the 
granules could be convective elements. This approach was developed by Vitense (1953). Using 
her equations and assuming a mixing length equal to the local density scale height, we estimate 

typical diameter = scale height = 200 km, 

AT < 1000°K, at is temperature excess, 

t 5 minutes, t is mean time, 

which are not in exact agreement with the observations but come close enough to provide some 
support for the convection hypothesis. 
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If the bright granules are convective elements, they should be rising, as they are at a higher 
temperature than the surrounding gas and therefore less dense, presuming that the two are in pres- 
sure equilibrium. Evidence indirectly confirming this is offered by Evans and Michard (1962b), who 
find that 70 percent of the bright continuum features they observe are associated with violet shifts 
in several weak spectral lines. This evidence is of marginal value, however, because the correla- 
tion is low and because the continuum and lines refer to different heights in the photosphere. 

The granules observed near the limb of the disk are higher in the photosphere than those ob- 
served near the center. Whitney (1963) suggests that limb granulation may be evidence of acoustic 
waves in the upper photosphere. 

Next, consider the determinations of upper photospheric and chromospheric velocity fields 
made by Evans and Michard (1962a, 1962b, 1962c) based on an extensive series of observations at 
the Sacramento Peak Observatory. From measurements of the Doppler shifts in a large number of 
spectral lines— which extend far in the direction of the spectrograph slit and appear "wiggly” in 
the photographs— the authors derive a r.m.s. value for what is called the "random turbulent ve- 
locity"^. Then they assume that the stronger lines are formed farther out in the atmosphere, and 
conclude from their observations that ^ is increasing with height over the region where all the 
lines are formed. 

To see the relevance of these observations for this study, we must inquire just how isotropic 
is the velocity fields . Where it is highly isotropic, there is turbulence that can generate acoustic 
waves. Where, on the other hand, the velocities are predominantly radial and not too large, these ra- 
dial motions may be quite consistent with the picture evolved in the next section, where acoustic 
waves with a net flux vertically upward pass through an atmosphere which, to a first approxima- 
tion, is in hydrostatic equilibrium. The evidence of Evans and Michard is not conclusive, but it 
does support the view that the velocity field ^ is not isotropic. Center -to-limb observations of 
each line for /u, < 0.6 yield decreasing values of When we consider that s is thought to increase 
with height, as already mentioned, the only explanation, other than one based on instrumental- 
seeing problems, is that the horizontal component of velocity is less than the vertical component. 

The evidence is weakened by the authors' belief that atmospheric seeing effects may, in fact, pro- 
duce these observations and thus render the above conclusion invalid. Nevertheless, these obser- 
vations provide the best available clue to the extent of photospheric-low chromospheric turbulence; 
their evidence is against isotropic turbulence in this region. Table 5 gives the results of the Evans 
and Michard study. 

Unfortunately, the question of turbulence has not been settled, as the work of Holweger (1967) 
shows. He concludes, after an analysis of many weak lines, that the horizontal-velocity component 
exceeds the vertical one and that the micro-turbulent velocity ^ decreases with increasing height 
in the photosphere, contrary to Evans and Michard. 

An interesting feature of the solar atmosphere, regarding wave propagation, is the oscillatory 
vertical velocity field near the temperature minimum. The heights in the atmosphere at which these 
oscillations are observed are inferred from depth-of-formation estimates of the spectral lines 
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Table 5 



Center-to-Limb Variation of ^ (km sec-i) from Evans and Michard (1962a). 
First Set: Tangential Slit 


Line 


Fel 5324.193 


1.00 

0.37 


0.76 

0.34 


0.65 

0.37 


0.52 

0.26 


0.29 


0.26 


Line 


Roland 


1.00 

0.64 

0.38 

0.16 

Intensity 

Mgl 5172.700 

0.61 

0.46 

0.34 

0.25 

20 

Fel 5171.612 

0.39 

0.40 

0.34 

0.24 

6 

Til 5173.751 

0.31 

0.43 

0.30 

0.25 

2 

Nil 5168.665 

0.31 

0.41 

0.34 

0.20 

1 

Fel 5170.770 

0.35 

0.34 

0.27 

0.26 

0 



Second Set: Radial Slit 



Line 


M 



Roland 

1.00 

0.60 

0.46 

0.30 

Intensity 

Mgl 5183.621 

0.47 

0.35 

0.41 

0.29 

30 

Mgl 5172.700 

0.43 

0.37 

0.39 

0.19 

20 

Fel 5171.612 

0.32 

0.29 

0.38 

0.20 

6 

Til 5173.751 

0.32 

0.35 

0.46 

0.27 

2 


Roland 

Intensity 


studied. The papers by Evans and Michard (1962c), Leighton, et al. (1962), and Noyes and Leighton 
(1963) are all detailed studies of this phenomenon. A summary of the salient results of both shows 
the following; 1) There is a periodic, quasisinusoidal variation of vertical velocities in medium- 
strength lines with periods from 200 to 300 seconds that tend to decrease with increasing line 
strength, i.e., with increasing height in the atmosphere. In addition, the r.m.s. velocity amplitudes 
increase with line strength, i.e., with height. The velocity amplitudes are always less that 1.0 km 
sec" ^ 2) Further support for the velocity-field variations comes from intensity fluctuations in the 
cores of medium- strength lines, with a period roughly the same as for the velocity field and decreas- 
ing with increasing line strength. 3) The observed oscillations disappear as one scans the disk from 
center to limb, indicating that the oscillatory velocity field is largely normal to the solar surface. 

The next section takes up the theoretical interpretation of the above three types of observa- 
tions; their interpretation in terms of various kinds of wave propagation will be particularly direct. 
We wish to relate some of these observations to our previous considerations on isotropic turbulence. 
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There is evidence that the non-thermal velocity fields studies by Evans and Michard are aniso- 
tropic. Further observation of oscillatory Doppler shifts and brightness fluctuations in the cores 
of medium-strength lines led Evans and Michard (1962c) to remark, "A sizable fraction of the 
kinetic energy is contained in the distinctly organized non-turbulent motions of many unrelated 
small elements." What fraction of the kinetic energy of the photospheric velocity field is tied up 
in isotropic turbulence is still an open question, but some kind of non-thermal anisotropic motion 
does seem to be coupled to and lying above the disk granulation. That there might be more nearly 
isotropic turbulence in the deeper layers of the photosphere is suggested by the results of Evans 
and Michard (1962a) which show that the ^ derived from the weaker lines actually increases when 
M changes from 1.00 to about 0.6. Adding the deduced oscillatory velocity field to the upper part 
of the region gives a picture of the photosphere -low chromosphere we refer to in the next section 
in discussing the generation and propagation of various wave modes through this region. 

On a much larger scale than the granulation is the recently discovered "supergranulation net- 
work," reported by Leighton, Noyes, and Simon (1962). This is a system of large cells easily seen 
in the Ca >^6103 line. They are approximately circular near the center of the disk, have a mean 
diameter of 1.6 x 10‘* km, horizontal motion of matter from center to rim with a mean velocity of 
0.5 km/sec, cell lifetime in the range lO"* - 10® sec (difficult to estimate due to the Earth's rota- 
tion), and an average spacing between cell centers of 3.0 x 10'* km, giving some 5 x 10® cells over 
the solar surface. A similarity in appearance was revealed between the cells observed in Ca \6103 
and the solar surface seen in the central emission component of the Call K line. 

The supergranulation network will not play any direct role in the calculations of this report, 
not because it is unimportant, but because its role and, hence, its importance in the-chromospheric- 
energy-balance problem is difficult to evaluate. The super granulation appears to be a manifestation 
of convective motion that could originate deep in the convection zone, with material moving essen- 
tially upward in the center of the cell and down along the rim. This motion appears to be more 
ordered than the rather random motions of the granulation, which exhibits horizontal as well as 
vertical motions. Exactly what coupling might exist between the two, and hence what role the super- 
granulation may play in mechanical energy generation, is an unsolved problem. However, there is 
one atmospheric phenomenon with which the supergranulation appears morphologically connected: 
the spicules, already mentioned in Chapter I and in the discussion of the chromosphere and corona. 

Two comprehensive studies of spicules are those by Athay (1959), based on data obtained by 
Dunn (1956), and Beckers (1968). The reader is referred to PSC, where the Athay study is re- 
ported in context with many other earlier studies of spicules. The Beckers article is a review of 
more recent work and offers a comparison with the older studies. Although a fairly consistent 
picture has emerged correlating high-resolution spectroheliograms with Ha bright mottles 

O 

seen at line center between ±0.2A, and some information of a statistical nature has appeared re- 
gardii^ the numbers, rise-times, etc., of spicules, their dynamics are not understood, apart from 
some qualitative theories discussed in the Beckers review. For this reason we shall not discuss 
spicules any further here, except to mention that the above correlation between the chromospheric 
network in and H a (presumably originating in spicules seen on the disk) also correlates with 
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the supergranulation network and regions of high chromospheric-magnetic-field strength. Even- 
tually a dynamical theory will have to correlate all these observed phenomena. 

We return to the assumptions of HE and spherical symmetry mentioned at the beginning of 
this section. We shall not explicitly use the HE assumption, but since we rely on the PSC model 
to start our calculations in this region, and since HE was assumed in constructing the PSC model, 
we implicitly assiune HE. The evidence presented here suggests that we need not consider iso- 
tropic turbulence in the low chromosphere. However, oscillatory motions do occur in the vicinity 
of the temperature minimum, with velocity amplitudes always less than 1.0 km sec"^ . This con- 
trasts with a thermal velocity for hydrogen of the order of 10 km sec ' *. On this basis, we shall 
ignore the pressure perturbation of the oscillatory velocity and retain the assumption of HE. 

There was some evidence in the studies cited that the velocity field would lose its "standing 
wave" character around the temperature minimum and move outward with increasing velocity. In 
the next section it is argued that these running waves are acoustic waves in the chromosphere. 

This indicates that to a first approximation HE applies to the state of the low chromospheric gas, 
which is perturbed by a progressive sound wave that builds up into a shock wave. 

Spherical symmetry is also assumed for the PSC-model low chromosphere. Barring signifi- 
cant penetration of convective elements from the low photosphere into the chromosphere, or sign- 
ificant penetration of spicules below h = 1300 km (for which there is no direct, observational evi- 
dence), spherical symmetry in the low chromosphere is probably a reasonable assumption. On 
the other hand, House (1961) concludes from a study of emission in Mgl that departures from 
spherical symmetry must occur somewhere between 500 and 1000 km. This may be correct; the 
effect would be to produce deviations from the PSC-model low chromosphere we have adopted. 
However, in the absence of a complete chromospheric model based on possible deviations from 
spherical symmetry, we shall retain the PSC model for h < 1300 km, along with the implied as- 
sumption of spherical symmetry. For h > 1300 km, we do not invoke spherical symmetry, which 
is quite invalid in the upper regions, except in our final calculation for T^(h), which must be in- 
terpreted as some suitable average temperature in this region. 

One final idea should be mentioned before we try to tie the observations of this section into 
our mechanical- energy generation and dissipation model for the next. Bhatnagar, Krook, Menzel, 
and Thomas (1960) point out that the large extension in space of many stellar atmospheres can be 
seen as a momentum transfer or an energy transfer. These concepts can apply simultaneously. 
There could be momentum transfer due to the chromospheric velocity fields, which would push the 
gas radially outward, and energy transfer due to some effective dissipation mechanism, which 
would cause extension by producing higher temperatures. Some coupling should exist between the 
two. Which process is dominant in producing the observed extension of the solar atmosphere? The 
authors give evidence for energy transfer. A rapid rise of temperature with height would favor 
energy transfer, and go far toward eliminating the need for isotropic turbulence to explain the widths 
of many spectral lines. This would also support the picture of anisotropic, preferentially radial, 
mass motions in the chromosphere, as suggested in this section. 
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F. Generation and Dissipation of Mechanicai Energy 


In the introductory statement of the problem it is stated that several sources of energy for 
heatii^ the outer solar atmosphere have been proposed. We consider these sources in turn, and 
finally select the most probable. First, we should estimate the total energy required to balance 
the net radiation losses from the chromosphere and corona. In chapter IV we discuss the esti- 
mates that have been made to date; that discussion indicates the need for a flux of 10^ ergs cm"^ 
sec"^ at the temperature inversion. Practically all the radiant energy loss is due to the chromo- 
sphere alone. 

Hoyle (1949) was the first to calculate how the corona, but not the chromosphere, might be 
heated by accretion of interstellar particles as the Sun moves through the galaxy. Billings (1966) 
has reviewed the development of this idea to date; he concludes that, for the corona, the accretion- 
heating hypothesis is weak though not entirely dead. He points out that no accretion-heating theory 
proposed can explain the very high-temperature, low-lying, yellow-line regions of the corona and 
that no provision has yet been made for the efflux of matter due to solar wind. Considering these 
difficulties and the fact that most of the radiant energy loss occurs in the chromosphere (a fact not 
considered in the accretion studies and of primary importance to us here), we reject the accretion 
hypothesis. 

Cayrel (1963) suggested the photospheric radiation field as a possible means of raising the 
temperature in the low chromosphere. He claims that, in an atmosphere in radiative equilibrium 
which is optically thin in the continuum and is subjected to a strong radiation field (the photosphere), 
a departure of the continumn source funtion from the Planck Function can produce a temperature in 
the outer region equal to the effective temperature of the continuum -producing region. Unfortun- 
ately, I shall show in my treatment of the H" problem in Chapter V that this mechanism will not 
work in the low chromosphere. So the Cayrel treatment must be modified to allow for some me- 
chanical heating. The relevant photoionization reaction is 

H" + hi^ - H + e . (51) 


As for the dominant collisional reaction, Pagel (1959) showed that estimates by Dalgarno (in a 
private communication to Pagel) of the reaction rate k = 3.0 x 10"*°cm ^ sec"^ for the associative 
detachment reaction 


H‘ + H - + e (52) 

would cause reaction 52 to dominate the ionization equilibrium of H" in the solar atmosphere out 
to ■^5000 - 0.002. More recent theoretical calculations (Dalgarno, 1967) and laboratory measure- 
ments (Schmeltekopf et al., 1967) of the rate constant k for reaction 52 give k = 1.3 - 1.9 x 10"® cm^ 
sec"' for a temperature range of 300 - 8000°K. Thus, Pagel's conclusion remains valid. That is. 
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the reaction 52 now predominates out to Tjoqq - 0.0003. Table 1 shows that the collisional reaction 
predominates out to about 200 km. Beyond 200 km, photoionization may predominate. In Chapter 
IV we shall consider the role of collisions and photoionization in fixing the temperature. 

Solar activity can be ruled out for many reasons, the most obvious being that conditions in the 
steady- state chromosphere (as here defined) are essentially independent of time on the solar cycle, 
whereas solar activity is strongly dependent on this time. 

We have mentioned that spicules seem an unlikely cause of coronal heating, because of their 
subsonic speeds and low temperatures. Yet they might play a role in heating the upper chromo- 
sphere. From PSC or the Beckers' (1968) study it appears that 30 km sec”^ is a typical spicule 
velocity. Following Thomas (1948a), we note that this velocity exceeds the adiabatic speed of sound 
for all temperatures below 50,000°K, depending on the ratio of specific heats y which we use to cal- 
culate from the adiabatic sound speed formula , where R = 8,31 x 10’ ergs °K~i 

mole"' is the gas constant and /j- is the mean molecular weight and using = 30 km sec"'. Thus, 
even if spicules are not found below h = 1300 km, they may play a role in heating the upper chro- 
mosphere by the mechanism of shock dissipation. Also, we must remember that spicules may 
originate at some height below h = 1300 km. 

If the spicules are supersonic jets, heating the ambient chromospheric gas through shock dis- 
sipation, then the problem lends itself to a shock -wave treatment, which will be developed pre- 
sently. Since the origin and dynamics of spicules are not yet well understood, we shall treat the 
problem no further from this viewpoint, but merely keep in mind that spicules may play some role 
at least in upper chromospheric heating. A qualitative theory of spicule formation has been devel- 
oped by Kuperus and Athay (1967), in which thermal conduction downward through the chromosphere 
corona transition region is the source of energy that generates the spicules. Only a detailed study 
can settle this question. 

We must look for mechanical sources lying in or beneath the photosphere. This source must 
satisfy two conditions: 1) Enough energy must be provided to balance the net radiation losses of 
the outer atmosphere, already mentioned. 2) We must show how this energy is generated, pro- 
pagated, and dissipated. 

We tentatively identify the energy source with the turbulent motions of the gas at the top of the 
convection zone. That some kind of turbulent motion are present is supported by observations of 
granules and by theoretical considerations; both should reveal something as to where these motions 
occur and as to their nature, isotropic or anisotropic. 

Unfortunately, continuum observations have so far not been analyzed to determine a 
"random turbulent velocity", because resolution has not been sufficient, even in the Project Stra- 
toscope work discussed in the previous section. Thus, in the region where disk granulation is ob- 
served at 300 - 400 km below the h = 0 surface, many current estimates of the energy generated 
are based on the assumption that an isotropic turbulent velocity field is present, when this has not 
been established empirically. On the other hand, inferences of a random turbulent velocity, made 
from studies of weak spectral lines, invariably refer to higher regions of the photosphere (which 
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may be important for mechanical- energy generation and which we shall discuss later). At present 
the center of interest lies in the deeper region; thus we are seeking evidence of random turbulent 
velocity at levels where h ~ -350 km. 

Empirically, only vertical velocities have been inferred for this region. (This type of observa- 
tion is discussed in Section E.) We assume that there is some kind of turbulence here because of 
the locally high "Rayleigh number" (a dimensionless parameter defined below). Following Unsold 
(1955) and Schwarzschild (1959), we note from laboratory experiments that the Rayleigh number has 
a critical value 10® below which, in the range of values 10® - 10®, convective motion in the form of 
laminar flow in stationary, vertical, hexagonally cross -sectioned Benard cells occurs and above 
which some form of non- stationary convection occurs. The Rayleigh number Ra is defined by 


g^4 -y - I 

“ i^k T ’ (51a) 

where -£ is the density scale height across which the density decreases by a factor of e, is the 
kinematic viscosity coefficient, k is the radiative conductivity, and A(VT) is the excess of the tem- 
perature gradient over the adiabatic gradient. The other symbols have their usual meaning. Using 
the same values for these parameters as Schwarzschild (1959), which are 

T = 8200°K, = 300 km , 


k - 2 X 10^^ cm^ sec”', v - 2 x 10® cm® sec 


(52a) 


y - 1 

■ j- A(VT) 


1 

4000 


km” ' , 


which gives Ra - 10®. Thus stationary, laminar convection is ruled out. On the other hand, as 
Schwarzschild points out, solar granulation does not seem to be a structureless, entirely random 
(i.e.,, isotropic) phenomenon. He considers this inconsistent with the very high value of Ra ~ 10®, 
which implies a coihpletely random turbulent motion. Arguing that one should use "turbulent vis- 
cosity" in place of kinematic viscosity in Equation 51, he reevaluates Ra, only to find a new value 
too low for any convective motion at all, much less non- stationary convection. He concludes that 
the true Rayleigh number should lie slightly above 10®, because kinematic viscosity should be used 
in Equation 51 and the "Reynolds number" (see below) has a high value appropriate to turbulent flow. 
He adds that the turbulence is probably not random, on the grounds that such turbulence would re- 
quire using turbulent viscosity in Equation 51, which leads to a ridiculously low value of Ra. His 
final conclusion provides a consistent picture both from observational data and from arguments 
based on the Rayleigh number. 
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Although some anisotropy is highly probable, we shall assume here that these motions in the 
lower region of h ~ -350 km are isotropic. 

Further evidence for highly turbulent motions is provided by estimating values for the Reynolds 
number in this region. The Reynolds number is defined by 


Re 



(53) 


where is a characteristic length (taken to be the order of the density scale height as in Equation 
51), V is flow velocity, ij. is dynamic viscosity (related to kinematic viscosity v hy v = y-/p), and p 
the mass density. Minnaert (1953) suggested the following values for these parameters at the top 
of the convection zone; 

't = 10^ km , 

= 10”^ dynes sec cm' 

which differ somewhat from the values in Equation 52 but not enough to affect the conclusion below, 
based on a Reynolds number Re ~ 10*^, which follows from Equations 53 and 54. Minnaert's con- 
clusion is that such a high Reynolds number ensures turbulence at the top of the convection zone, 
because laboratory experiments show that transition from laminar to turbulent flow occur when 
Re > 10^. On this basis, we assume that a highly turbulent flow exists at the top of the convection 
zone. 

Before reviewing the theory of acoustic-energy generation by this turbulence, let us reconsider 
some previous statements. First, the choice of the characteristic length -i in Equations 51 and 53 
is somewhat arbitrary. In the original laboratory experiments where these two numbers Ra and Re 
were used, I referred to the thickness of the fluid layer and the radius of a pipe through which the 
fluid was flowing, respectively (Huang and Struve, 1960). Since the element's (granule's) size, the 
distance (Prandtl mixing length) through which a typical element moves before losing its identity, 
and the density scale height at the top of the convection zone are all of the same order of magnitude 
in the range 10^ - 10^ km, according to theories of this region already considered (Vitense, 1953; 
Bbhm-Vitense, 1958; Biermann, 1959), it is felt that a choice for I in this range is the best we can 
make; but there is some uncertainty. 

Consider the Lighthill (1952) theory for acoustic-power generation applied to an isotropic, ho- 
mogeneous, turbulent velocity field. It is assumed that all the acoustic energy produced will be 
radiated away with no feedback into the turbulent velocity field; this is justified as long as the Mach 
number for the turbulent motions is much less than unity. Theoretical estimates of turbulent ve- 
locities (Bbhm, 1963), based on convection theories which are admittedly rough approximations and 


(54) 


= 10" s 


gm cm 
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treat vertical velocity components only, give maximum values of about 2.0 km sec"*. This value 
is sufficiently small compared to the local velocity of sound in the low photosphere for Lighthill's 
assumption on the turbulent Mach number to be valid in this region. 

Lighthill starts with the basic equations of continuity and motion of fluid dynamics (Courant and 
Friedrichs, 1948), and derives an expression for density perturbation in a turbulent fluid as a func- 
tion of the stresses associated with the turbulent velocity field. This expression for density per- 
turbation has the form 


P (X. t ) - /^o " DI + QU , 


(55) 


where Di and QU represent dipole and quadrupole acoustic radiation fields, respectively, because of 
their formal similarity to the dipole and quadrupole integrals in the expansion of the electromagnetic 
potential; Pg is the equilibrium density of the gas at point x. The intensity of sound is 


i(x) 



s 




({p (z. t) 



(56) 


where t) " mean- square fluctuation in p and the other symbols have been de- 

fined already. The total power output is obtained by integrating the intensity over a large sphere 
of radius lxl , which encloses the region of sound generation. Dividing by the volume of the sphere 
gives the mean power density. 

Proudman (1952) performed the integration of the quadrupole term obtained by combining Equa- 
tions 55 and 56 for the case of isotropic, homogeneous turbulence. Lighthill evaluated the dipole 
term. Kuperus (1965) compares the results of the two integrations and shows that 



P 

q 


M' ^ , 


(57) 


where and P^ are the power densities for dipole and quadrupole acoustic radiation, respectively, 

Mj. = is the Mach number for turbulent motions with v^. the velocity of turbulence, and <?^is a 

number, evaluated by Kuperus, weakly dependent on the assumption concerning the correlation func- 
tion f ( r ), given by 


f(r) = (v; (r,) ■ V. (r,)> 


(58) 
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Here, Vj is the velocity component in the i -th co-ordinate direction, and ‘ is the 

mean value of v. ^Xi) ' ''i(x 2 ) subject to the condition r = Ixi ~ X 2 I . Specific forms for Equation 
58 are available for the Heisenberg correlation and the Gaussian correlation of velocities. Kuperus 
uses Proudman's work to show that a difference by a factor of only 2.86 exists in the choice of 
depending on whether we choose the Heisenberg or Gaussian correlation; he claims that experimen- 
tal evidence favors the Heisenberg. Evaluating </>with the Heisenberg correlation function for f(r), 
Kuperus concludes that 0.18 for = 0.15, making the quadrupole term dominant. We ten- 
tatively accept this conclusion with one reservation; the assumptions of isotropic turbulence and a 
Heisenberg velocity correlation have not been verified to date for the lower photosphere. Unno and 
Kato (1963) have also evaluated the terms of Equation 56, using the Lighthill theory with some mod- 
ifications. They note the sensitivity of the results to these unverified assumptions, including the 
possibility that sufficient anisotropy in the turbulent-velocity field might cause the dipole term to 
be as important as the quadrupole term (Unno, 1964). 

Lighthill (1954) supplemented the above calculations by deriving a more general form of Equa- 
tion 55 which includes the effect of anisotropic convection (this, presumably, "drives” the turbulent 
motions by agitating the gas at the top of the convection zone) on the total power radiated. He shows 
that the effect on the quadrupole radiation field is to multiply the QU term of Equation 55 by the fac- 
(1 - cos S )~^, where is the Mach number for convective motions with v^, the convective 

velocity, and S is the angle between and the direction of sound emission. The intensity due to 
quadrupole emission must then be multiplied by (1 - cos 8)'®. This correction is not included 
in the total mechanical flux estimate we give below, because turbulent velocities are assumed equal 
to convective velocities as a first approximation. A more detailed study should distinguish between 
the two and include the effect of the convective motion with a preferred direction. The effect of con- 
vective elements that penetrate beyond (above) the radiatively unstable lower photosphere will be 
discussed later in this section. 

Now we can consider some estimates of the net flux of mechanical energy generated by the 
Lighthill (1952) mechanism in the turbulent, radiatively unstable lower region of the photosphere. 
The effect of magnetic fields is ignored. The formula Proudman obtains for the acoustic power 
generated per unit volume by the QU term of Equation 55 is 


^0 <^t) 


3/2 


(59) 


where c is a number dependent upon our choice for the velocity correlation f(x) of Equation 58, 
is a characteristic length of the turbulent-velocity field, and the other parameters have already 
been defined. The variation in c is only a factor of 2.86 (as with 4 > of Equation 57), depending on 
whether we choose the Heisenberg or Gaussian correlation. Again we choose the Heisenberg cor- 
relation and, like Proudman, obtain c = 37.7. 
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Kuperus (1965) numerically solves Equation 59 for two different models of the convection zone, 
to obtain a rough estimate of the mechanical-energy flux generated. We use his estimates here. He 
chooses L = 500 km, which is reasonable in the light of our arguments for choosing a characteristic 
length-^ to evaluate the Rayleigh and Reynolds numbers. Thus Equation 59 becomes 


= 7.55 X 10- ’ M/ , 


(60) 


where Pq and are obtained directly from the convection zone models, and 

is calculated from these models, which specify T^,m, andy. The models chosen by Kuperus to cal- 
culate Pq are theoretical ones calculated by Biermann et al. (1959) and Heintze et al. (1964). The 
Biermann model is practically identical with the Vitense (1953) model, mentioned in the last section 
to establish similarities between observed granules and calculated properties of convective ele- 
ments. The Heintze (1964) model is called the "Utrecht Reference Photosphere" and is similar to 
the Heintze (1965) model, for which the temperature distribution is plotted in Figure 2. Isotropic- 
turbulence velocity is assumed equal convective velocity as a first approximation. 

Kuperus assumes that all acoustic power generated will ultimately propagate upward from the 
convection zone. Since half of the power initially generated will propagate downward, according to 
the early Lighthill (1952) theory, a reflection is needed to change the direction of those waves that 
initially propagate downward. Kuperus does not consider what fraction of these downward propagat- 
ing waves will "bounce off" the denser material below and propagate upward before appreciable at- 
tenuation. This omission is probably allowable because 1) some reflection can be expected, 2) the 
already mentioned correction factor (1 - cos S)"® shows us that more sound will be initially 
emitted upward than downward for a convective velocity directed predominately upward, and 
3) the maximum uncertainty of a factor of 2 implicit in this assumption is no greater than the un- 
certainties upon which the convection zone models or Equation 60 for the power generation are 
based. 

The calculation of the mechanical energy flux is reduced to an integration of the acoustic 
power output per unit volume over the thickness of the generating region. Thus, Kuperus numer- 
ically solves 


F 


m 


P, (h) dh 


(61) 


where Pq(t>) is obtained from Equation 60, and hj and are derived from the convection zone model 
chosen. The lower boundary, hj , occurs where the decrease of (due to decreasing and in- 
creasing inward) makes Pq(h) negligible. The upper boundary, occurs where radiative equili- 

brium is reestablished in the photosphere as one proceeds outward. The region of appreciable 
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acoustic-power generation is a thin layer about 100 km thick as shown by a simple estimate based 
on either model Kuperus used Jor the upper convection zone. 

The results of Kuperus' calculation of the flux F^, obtained from Equation 61, are as follows: 

= 2.9 X 10® ergs cm-2 sec-i from the Biermann model of the convection zone; ^ ^ 

2.7 X 10’ ergs cm~^ sec“^ from Heintze model, depending on whether upward or downward convec- 
tive velocities are used, as the Heintze model is an inhomogeneous model with columns of both 
rising and falling material. Kuperus simplifies the calculation by selecting a mean temperature 
~6250° for evaluating F^. Since does not vary strongly within the 100-km layer where the energy 
is generated, this simplification is justified. 

These three estimates of F^ are all equal to or greater than the various estimates of the net radi- 
ative loss from the outer atmosphere, given in Chapter IV. Although there are many uncertainties 
in the estimates of both quantities, we conclude that there is probably enough mechanical energy 
generated by the Lighthill (1952) mechanism to balance the net radiative losses of the outer 
atmosphere. 

That magnetic fields in the low photosphere probably do not affect the above results, we see 
from the following considerations. Kulsrud (1955) evaluated the quadrupole term of Equation 55 
for the more general case with a magnetic field. The generated power is enhanced by the turbu- 
lence components of the magnetic field. This enhancement is exhibited as a function of the Cowling 
number S, defined here as the ratio of the turbulent magnetic-energy density to the turbulent kinetic- 
energy density. We have 



where is the turbulent Cowling number and is the turbulent magnetic-field strength. For 
< 0.1, Kulsrud shows that the magnetic field plays a negligible role. We do not know the value 
of , but if we assume that it is of the same order as the steady-state, photospheric magnetic- 
field strength of about 1 gauss (Babcock and Babcock, 1955), we can estimate s^. Suppose we as- 
sume, for the low photosphere. 


Vj - 1 km sec”* , - 1 gauss, 


(63) 


- 10”’ gm cm' 


This gives 10”'* « 0.1, well within the limit for neglecting the magnetic field in the quadrupole 
term of Equation 55. Following the method of Kulsrud, Osterbrock (1961) evaluated the effect of a 
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steady-state magnetic field on the dipole term of Equation 55. He shows that H 2:10^ gauss is re- 
quired before the corrected dipole term will complete in magnitude with the quadrupole radiation 
term of Equation 60. Thus, ignoring the effect of magnetic fields on energy generation is justified 
outside active regions. It is probably not justified in the plage regions, but, since they are solar- 
cycle dependent and cover a small fraction of the photosphere, we do not consider them here. 

This completes our review of possible sources of mechanical energy for heating the outer solar 
atmosphere, except for a discussion to follow on the not yet well understood effects of penetrative 
convection. We conclude that acoustic waves generated by a turbulent-velocity field in the low pho- 
tosphere provide the most likely source in view of current observational data and theoretical models 
of the photosphere and the underlying convection zone. Now it is necessary to consider how these 
waves propagate. 

The term "acoustic" has so far been used to describe the waves generated. "Pressure" waves 
(generated in the presence of a gravitational and a magnetic field) would be as valid. The resulting 
energy flux obtained from the Lighthill theory, when generalized by Kulsrud for magnetic-field ef- 
fects, is seen to be almost independent of gravity and magnetic fields in the low photosphere. Grav- 
ity independence is due to the fact that only the weaker dipole term of Equation 55 is gravity- 
dependent. We have already advanced arguments for magnetic -field independence. The important 
point is that even though the generation of these pressure waves is practically independent of solar 
surface gravity and the photospheric magnetic field, the same might not be true of wave propaga- 
tion. Therefore, we must examine the problem of the propagation of an infinitesimal disturbance 
through a compressible plasma with gravity and a magnetic field. 

This problem has been studied by Ferraro and Plumpton (1958) for an isothermal atmosphere of 
infinite electrical conductivity. If we further assume zero viscosity for the solar atmosphere, and 
orient our co-ordinate system so that g— a vector terminating at (0, 0,-g)— lies along the z-axis,we 
can start with the basic equations of magnetohydrodynamics (the so- called Navier- Stokes equations 
with a magnetic-force term) and the adiabatic law and, following Ferraro and Plumpton, derive the 
wave equation governing the displacement £ of a particle from its equilibrium position. The equation 
is 


ill 

at2 


((,^2 + c^2) y 2 y + (y- 1) g(V • 7) - V(gr J + • • • 


[Ho ' (64) 


where all the old symbols have their usual meanings, the subscript "0" refers to the equilibrium 
value of the quantity, is the Alfven speed, and is the z- component of r. 
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Equation 64 is rather complicated, but by assuming plane wave solutions of the form 


0 = 


exp 




cot + k 


+ k y 
y ■’ 


+ k 




( 65 ) 


and substituting these solutions in Equation 64, we obtain the dispersion relation for the dependence 
of angular frequency co on wave numbers k^, k^ , and on physical conditions in the atmosphere 
through g, , and . From this dispersion relation we can determine what wave modes are likely 
in the photosphere and low chromosphere. The dispersion relation is itself quite complicated, but 
if we assume Hg = (0, 0, Hg) for a vertical magnetic field and let k^ = 0 for two-dimensional pro- 
pagation, we still obtain all the physical behavior of the more general case with simpler mathe- 
matical expressions. Consider the following special cases. 

Case 1 is one-dimensional propagation, obtained if = 0. The dispersion relation then be- 
comes either Equation 66 or 67: 


^2 -^ 2^2 _ ^ Case la (66) 

Here, our phase velocity Vp = w/k^ “ “ H is the Alfven speed. The wave is a transverse 

wave travelling vertically along the magnetic field. Also possible is 

^,2 _ ^2 ^2 _ yg]i ^ i = 0 , Case lb (67) 


This is a sound wave travelling vertically. Solving Equation 67 for k^, 




c^2) 1/2 , c; 
s / ' s 



( 68 ) 


Only for w > will the wave number k^ have a real part. Thus, only for " > will propagation 
of a pure sound wave be possible. We shall return to this point later. 

Case 2 occurs if k^ f 0, but g = 0. In this case, we neglect gravity, but include magnetic- 
field effects. The dispersion relation becomes 


j'^4 _ (^^2 + c^2) + k^2)<2;2 - C^2 (k^2 + k^2)J = 


(69) 
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which describes magnetohydrodynamic waves in an isothermal compressible plasma of infinite 
electrical conductivity. If the first factor in Equation 69 is zero, case la occurs. If the second 
factor is zero, two new types of independent modes occur. These are called the "fast" and the 
"slow" mode, depending on the choice of the ± sign in the equation for the phase velocity, which is 

V " i “ 4 c/c^2cOs2 , (70) 

where 0 is the angle between the wave vector and the z-axis. 

Case 3 occurs if we put q 0, thus neglecting the magnetic field, but put g 0. The disper- 
sion relation becomes 



s 


0 , 


(y - g _ 


(71) 


Equation 71 yields real only where cu > or , assuming that , which is true in the 

photosphere. Ti cu < the waves are called "gravity waves." 

We are now ready to discuss wave propagation in the photosphere and low chromosphere in 
terms of the above. Osterbrock (1961) took Equation 69 and showed that, in both the photosphere 
and the low-chromosphere outside regions of stroi^ magnetic fields associated with plages, the 
fast mode is practically a sound mode, in that which makes the second factor in Equation 

69 reduce to the g = 0 form of Equation 67 for a pure sound wave. This is evident if we note that 
is equivalent to , From Table 1, we see that the place in the low chromo- 

sphere where this condition is least likely to be met is the higher part, owing to the rapid decrease 
in density with height. Testing the above condition at h = 1250 km, with p - 1.5 x 10” '^gm cm- 3 
and y 7500^K, and assuming H = 1 gauss and y = 5/3, we obtain y 0.05. Thus Osterbrock's 

contention about the fast mode is surely correct below this height. 

Equation 70 with the minus sign shows us that, for »c/, so the slow mode is prac- 

tically an Alfven mode. We shall review arguments for the rapid damping of the Alfven mode in the 
photosphere presently when we come to study dissipation mechanisms. Here, we use the result due 
to Osterbrock's discovery that the Alfven mode is so rapidly damped in the photosphere that effec- 
tively no flux of mechanical energy passes through the photosphere in the form of Alfven waves. 

From this we conclude that the entire flux of energy through the photosphere is in the fast mode or, 
in effect, sound waves, provided that the effect of gravity is negligible. 

Whitaker (1963) argues that the effect of gravity cannot be ignored. Accepting Osterbrock's 
conclusion that the magnetic-field influence is negligible over most of the solar surface, he argues 
in favor of gravity waves with the dispersion relation given by Equation 71. Since sound waves of 
angular frequency cc < cannot propagate outward but are reflected back down into the photosphere 
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at the temperature minimum where ^^3 attains a maximum value, Whitaker argues that some of the 
energy in sound waves cannot propagate into the chromosphere, where the energy is needed. To 
show that this is a serious problem, let us evaluate for conditions around the minimum tempera- 
ture. Using 


g = 2.74 X lO"* cm sec" ^ = 4600°K , 


(72) 


7 = 1.2 , 




1 , 


after Whitaker (whose y lies between the adiabatic and isothermal values) the second Equation 68 
gives £^3 = 0.0243 sec"* . This corresponds to a critical period of = 259 sec. We saw in 
Section E that a characteristic period for the observed photospheric oscillatory velocity field is 
300 seconds. We also saw that a typical granule lifetime lay at or slightly above this value. As- 
suming that these times typify the periods of sound waves in the photosphere (since we lack theo- 
retical calculations of the frequency spectrum of the acoustic waves generated by turbulence), we 
see that much of the energy in the sound mode may not be able to pass through the temperature- 
minimum barrier. Thus, we must inquire anew whether sound waves with dispersion relation given 
by Equation 67 with g 7^ 0 are still the most likely mode for transmitting energy to the chromo- 
sphere, or whether the gravity mode may do this more efficiently. 

The gravity mode favored by Whitaker and described by the dispersion relation (Equation 71) 
also encounters difficulties. These two-dimensional waves cannot propagate outward when co > 
where is given by Equation 71. If we evaluate for the same conditions (given by 72) for which 
we evaluated ^ we find that = 0.0181 sec"*, corresponding to a critical period for gravity 
waves of = 348 sec. This value lies above the characteristic period of 300 seconds which seems 
to dominate oscillatory processes in the photosphere. Thus there can be no vertically propagating 
gravity waves of T y 300 sec; the same type of argument used to question the efficiency of sound- 
wave transmission through the photosphere can be applied to gravity waves. In one respect, the 
situation is more severe for gravity waves. Whereas the temperature minimum, where and 
are evaluated, offers the greatest resistance to sound-wave propagation, it offers the least resist- 
ance to gravity-wave propagation, if we hold y and g constant and consider the dependence of "g on 
. Thus, the critical frequency above which gravity waves cannot propagate becomes lower on 
either side of the temperature minimum, excluding a broader range of frequencies and raising the 
question whether an appreciable flux of gravity waves generated in the low photosphere could reach 
the temperature minimum. 

The fact that raises a further problem. For "3 = 0.0243 sec"* and "g = 0.0181 sec"*, 

the corresponding periods are T3 = 259 sec and = 348 sec. Thus, periods between and 
define a range of non- outwardly propagating disturbances. Yet this is the very range of periods 
that dominate oscillatory processes in the photosphere. What, then, are the actual facts about this 
range? Moore and Spiegel (1964) suggest that these photospheric velocity-field oscillations are a 
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direct manifestation of non-propagating, standing waves in the region of the temperature minimum. 
They further suggest that convective elements, penetrating above the upper boundary of the radi- 
atively unstable region to the h = 0 region, may be able to generate these oscillations. Such pene- 
trative convection is to be expected (Moore 1967) from the presence of an inertial term in the equa- 
tion of motion, which the mixing-length theory used for our convective zone models ignores. 

Lighthill (1967) argues that one would expect this penetrative convection to generate gravity 
waves in the upper photosphere, assuming that these motions can be represented as a highly ani- 
sotropic, turbulent-velocity field. Indeed, he estimates that gravity-wave generation in the upper 
photosphere should be so efficient as to produce as much energy in gravity waves as the isotropic 
turbulence mechanism produces in acoustic waves in the low photosphere. However, he does not 
say what fraction of this energy flux can propagate through the temperature-minimum region. Thus, 
these gravity waves may be standing waves; then, Lighthill's suggestion recalls that of Moore and 
Spiegel. 

Uchida (1967) solved an eigenvalue problem to determine what waves can be trapped in a re- 
sonant cavity consisting of the temperature "trough" region around the minimum temperature and 
bounded on the sides by the "walls" of the supergranulation elements, which are reinforced by 
strong magnetic fields around the rim. He concludes that a standing-gravity-wave oscillation will 
result in a horizontal wavelength of 2000 - 3000 km and a period of about 300 sec. This conforms 
surprisii^ly well with observations of the oscillatory velocity field described in the last section. 

On the evidence presented so far, the general features of gravity waves and oscillatory veloc- 
ity fields in the solar atmosphere might seem to be understood, suggestii^ more precise ques- 
tions such asr 

1) What fraction of the mechanical energy generated by the convection zone goes into the grav- 
ity modes? 

2) What fraction of this energy is available in propagating frequencies to heat the chromosphere? 

Unfortunately, the picture may not be so simple. Souffrin (1966) studied the general problem of 
the response of an atmosphere to an underlying turbulent convection zone. He introduces radiative 
dissipation and describes the response of an isothermal atmosphere as a function of distance from 
the convection zone. For the solar photosphere, he concludes that gravity waves are unlikely, as 
the radiative relaxation time of this region is too short. 

It is clear that the related problems of penetrative convection, gravity waves, and oscillatory 
velocity fields have not yet been solved. So we shall not consider what role gravity waves might 
play in heating the chromosphere, as it is still not certain that there is a flux of gravity waves into 
the chromosphere. 

In spite of the uncertainty as to how much energy will pass through the temperature minimum 
in sound waves, they seem the most likely agent for carrying energy from the convection zone to 
the outer atmosphere. Now we must find a means for the dissipation of this energy. 
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As already stated, Alfven waves are so rapidly damped in the photosphere that they become 
negligible. Osterbrock (1961) defines a damping length C as the distance in which the energy flux 
drops by a factor e . For the damping of Alfven waves due to electrical conductivity and viscos- 
ity, he uses a damping length $i, from van de Hulst (1951), 
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where v is the wave frequency, cr the electrical conductivity, and m the kinematic viscosity. For 
"frictional" damping due to collisions between positive ions and neutral atoms, he uses a damping 
length , from Piddington (1956), 
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(74) 


where v is the wave frequency, v the ratio of the mass density of neutral atoms to the mass density 
of ions, and the mean collision time for an atom to lose its momentum in collisions with ions. The 
net damping length for Alfven waves is given by 


_L - _L J_ 

" ^2 


(75) 


Osterbrock obtains the electrical conductivity for electron collisions with positive ions fromSpitzer 
(1956) and for electron collisions with neutral atoms (of importance in the photosphere) from his 
own previous work (Osterbrock, 1952). He uses Edmonds' (1957) table of viscosity coefficients, 
extrapolated from the convection zone through the upper photosphere and low chromosphere, as- 
sumes H = 2 gauss outside plage regions, combines this with p from Minneart's (1953) model pho- 
tosphere and the van de Hulst (1953) model chromosphere to obtain , and takes = 1.2 

X 10"^ sec“‘ for the wave frequency. He finds that Sj « $2 photosphere and low chromo- 

sphere and that ~ 1.0 km in this region. For reasons already given, we cannot accept the van 
de Hulst model chromosphere; however, the damping length is so small that Alfven waves 

can certainly be ignored in the photosphere, where wave generation occurs. In the chromosphere, 
the fast-mode hydromagnetic wave (which is essentially a sound wave outside active regions) will 
be the only survivor among the modes generated in the photosphere with dispersion relation given 
by Equation 69. 

Osterbrock has also estimated damping lengths for sound waves in the photosphere and chro- 
mosphere. He adopts a result due to Lehnert (1959), which considers the effects of viscosity, 
electrical conductivity, and the frictional damping already mentioned. Since the viscous effect is 
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predominant for sound waves, the damping-lei^h formula Osterbrock uses is 


3c/ P 
167t^ fi 


( 76 ) 


where is the damping length, and the other symbols have already been defined. Application of 
Equation 76 to the photosphere and chromosphere readily shows that even viscous dissipation is 
negligible. We obtain 10® km, using v 10“^ sec, the Edmonds' viscosity tables, and the (p,t) 
values from any of the atmospheric models we have mentioned. Thus, sound waves will not be ef- 
fectively damped in the photosphere or low chromosphere by any of the above, linear, damping 
mechanisms. We are led to consider shock dissipation of these sound waves. 

If shock dissipation of acoustic waves is the mechanism for heating the chromosphere, two ques- 
tions arise: 1) Under reasonable assumptions for the depth of generation of acoustic waves, will they 
build up into shock waves over a short enough distance from the generation zone to produce the tem- 
perature inversion? 2) Can we express the dissipation of shock waves as a function of local atmos- 
pheric parameters so that we can solve the final energy conservation equation for temperature dis- 
tribution, subject to a successful solution of the radiation-transfer problem? 

First, consider how the shock wave builds up. Any finite -amplitude disturbance traveling in a 
compressible fluid builds up into a shock wave because the velocity of particles in the vicinity of 
maximum compression (the condensation phase) exceeds that of the particles elsewhere (the rare- 
faction phase), leading to the familiar saying, "the crest of the wave catches up with the trough" 
(Courant and Friedrichs, 1948). Using this criterion to define the appearance of a shock wave, 
Kuperus (1965) estimates that a sound wave travels approximately 500 km in the photosphere before 
the crest and trough coalesce. Bird (1964c) has obtained numerical solutions of the nonlinear equa- 
tions of continuity, momentum, and energy for a spherically symmetrical fluid, and obtains a forma- 
tion distance of the order of 500 km for an isothermal photosphere at T = 6000^K. 

In view of these results, and recalling that the sound waves for which we give mechanical flux 
estimates are generated in a region about 100 km thick centered at h ~ -350 km, this implies 
the existence of fully developed shock waves above h ~ 100 km in the chromosphere. If acoustic 
waves are also generated by penetrative convection in the upper photosphere, the resulting shock 
waves would presumably develop higher in the chromosphere. Since does not seem to rise be- 
yond the minimum value until '^5000 = 16'^ is reached, according to Curtis (1965), and since 
= 10" lies at 150 £ h < 250 km for the models of the low chromosphere discussed, the shock wave 
first appears predictably, somewhat below the level where the temperature rise is thought to begin. 

It is convenient to introduce at this point a new parameter, "shock strength "^s," defined in 
terms of density "jump" ~ Po across the shock front, where p„ is the maximum value assumed 
by p during a cycle and Pq the value of p in front of the shock wave; 
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From the RanMne-Hugoniot conditions (Courant and Friedrichs, 1948), it follows that shock 
strength is related to shock Mach nmnber by 







(78) 


Thus, either Vs or suffice to indicate the shock strength. 

To include shock dissipation in our final energy conservation equation, we need an expression 
giving energy loss per unit volume and time as a function of local conditions in the atmosphere. 
Bazer and Ericson (1959) derived such an expression for so-called "weak" shockwaves, defined 
here as < 1.0. They give a general derivation for hydromagnetic shocks; their results are given 
here, in the notation we shall use consistently hereafter; 
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where D is the mechanical energy input in ergs cm"^ sec" * , S is the period for "shocking" the gas 
( a; = f" ' , where f is the frequency with which shock waves pass through the local region), m is the 
mean molecular weight, and the other symbols have meanings already given. The above expression 
is for one-dimensional propagation along a magnetic field; otherwise, the term becomes 

more complicated. In fact, this second term is important only where magnetic -field effects are 
important. Since we ignore magnetic field effects (a resonable procedure in the low chromosphere 
where « 1.0), the expression for D that we use is 
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Finally, if we choose y = 5/3, S = 300 sec, and = 1, we obtain 

D = 1.025 X 1Q5 p 77/ . (81) 

It might be argued that these values for and m should be different. For example, y may 
vary across the shock wave, particularly where ionization occurs. Also, we have chosen for our 
constant y the adiabatic value for a monatomic gas. Weymann (1960) argues that radiative cooling 
behind the shock front occurs so rapidly above = 1 that one should there use the isothermal 

value y = 1. If we accept his result, we must consider the variation of y over a "shocking" cycle. 
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We do not do this, because 1) the extreme isothermal value 7=1 chaises Equation 81 by a factor 
of 2. 2 only, and 2) any mean r we might use in that equation must lie well above the isothermal 
r = 1 to be consistent with the Rankine-Hugoniot condition 


7+1 

~P^ - 7- 1 ■ 


and our estimates of 



which are tabulated in Table 6 and obey the condition < 1. Rather than conduct a detailed analy- 
sis for a resulting mean 7 (which will not significantly affect Equation 81), or arbitrarily 4)ick some 
value 5/3 > r > 1, we merely use the adiabatic 7 = 5/3. We justify our choice for 2 j and m on the 
grounds that they are representative values, averaged over the region we are studying. The error 
we introduce by ignoring their variations will not strongly affect our final (h). 

The estimate of is the weakest point in our development. Neither , nor equivalently, M^, 
has yet been directly measured and we have only weak theoretical grounds for estimating it. We 
shall base our estimates of on the work of Osterbrock (1961), which is the most thorough treat- 
ment of this problem available for the low chromosphere. He derives a differential equation for 
the mechanical flux gradient in the form 


^ . (82) 

where d(?]^ ) is the energy input given by Equation 79, with the additional effects of refraction included. 
He also has an expression for of the form 

= f(^s) - (23) 

which includes refraction effects. Assuming an atmospheric model for (p, t) and an initial value 
for T]^ (or F^) at h = 0, we can solve equations 82 and 83 for F^ (h) and (h) . From [p, 7, v^), 
using Equation 79, we obtain D(h) . 

Osterbrock assumes an initial value for F_^ of 3.0 x 10’ ergs cm"’ sec"* at h = 0 (which lies 
in the mid-range of Kuperus' three estimates given earlier in this section) corresponding to (0) 
is 0.32. Osterbrock uses the van de Hulst (1953) model chromosphere for his {p, T). We 
have rejected this model, as it was based on LTE; indeed, the {p, T) variation is quite different 
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from that of the PSC model we have adopted. While T (h) for PSC is greater than T (h) for van de 
Hulst, the scale height in the van de Hulst atmosphere is somewhat larger than in the PSC atmos- 
phere, so p(h) for van de Hulst is greater than /s(h) for PSC. When h has increased to 1000 km, the 
difference in D is an order-of magnitude because the density inequality is greater than the tempera- 
ture one, but the general variation in^;^ (h) with height h obtained by Osterbrock is still valid. 

We illustrate Osterbrock's for magnetic -field strengths of h = 0.5, 2.0 gauss in Table 6. 

Our own choice for , also given in Table 6, is guided by the following considerations: 1) 
should grow from a small value 2 0 at h = 0 to some value < 1.0 in the low chromosphere. This 
initial trend is also predicted by Bird (1964c). 2) Osterbrock admits that there may be some dis- 
sipation for Vs less than his minimum value of 0.32 (where he concludes the shock waves are first 
fully developed). Frisch (1967) concludes that Osterbrock overestimates in these low regions 
of h < 300 km. Combining these observations, we conclude that it would be better to choose v^ (0) 

= 0 and let v^ grow rapidly with height. 3) The condition (0) < 1.0 must be fulfilled, or our use 
of Equation 79 for D becomes invalid. 

The calculations of both Osterbrock and 
Bird support this view. On these bases, we 
have chosen the followii^ v^ (h); 


V^ = (l - e , P - 400 km . (84) 

These values appear in Table 6. They exhibit 
the proper trend and remain within the proper 
limits. This is as reasonable a choice for 
as we can make; note, however, that the esti- 
mate oi Vs is uncertain. 

In the region we call the upper chromosphere, where h > 1300 km, we could assume that v^ - 
according to Equation 84 or, in view of the sharp transition between the low chromosphere and upper 
chromosphere, we could examine the upper part separately. Bird (1964a, 1964b) argues for a con- 
stant Vs in this region, and in the upper part of the low chromosphere for that matter. He derives a 
condition for the gradient AT^/dh in a spherically symmetrical atmosphere with gravity, starting 
with the basic hydrodynamic equations and arriving at 


Value of Shock Strength Parameter Vs • 
(Tj^are from Osterbrock, 1961.) 


h (km) 

"7s 

be = 0-5 

77b = 2-0 

0 

0 

0.32 

0.32 

250 

0.46 

0.41 

0.43 

500 

0.71 

0.62 

0.61 

750 

0.85 

0.71 

0.70 

1000 

0.92 

0.78 

0.78 

1250 

0.96 

— 

— 


da _ g /2g/g _ _L_ 

dh " 4(1+ a) dh ^ ’ 


(85) 


where we have restricted his result to plane parallel atmospheres, (a = - 1 is a new measure 

of the shock strength, and the other symbols have their usual meaning). From Equation 85, Bird 
argues that, for a negative quantity in square brackets, the large dT^/dh > 0 will cause a to drop 
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rapidly and make the dissipation also drop, until dT ^ /dh has dropped to a value where da/ dh = 0. 
The converse is also true, he notes. Thus, he predicts da/dh = o is the criterion for stability; 
thus a (and hence and v^) is constant. 

Bird's argument is weakened by the fact that his criterion cannot be applied rigorously at a 
point, but only shows a trend. 

To see this, note that assuming da/ dh = 0 in Equation 85 gives 


dh 


2ga 

R 


(h in cm) 


( 86 ) 


For fj. - 1 and the solar surface gravity, we obtain dT /dh - 65°K km~'. This is much too large for 
the low chromosphere. In fact, for reasonable values of (T^, dT^/dh), Equation 85 tells us that 
da/dh > 0 in the low chromosphere, as we have shown already. For h > 1250 km, however, the 
temperature may rise at a rate not appreciably different from 65°K km' ' . Thus, in the upper chro- 
mosphere, the condition da/dh = 0 may apply. 

We shall assume that a is constant for h > 1250 km. Actually, we have no reliable model for 
the upper chromosphere, so it would be useful to consider a range of values for ■ We shall com- 
pute models for the upper chromosphere for = 0.5 and = 1.00, realizing that as approaches 
1.00 the expression for D in Equation 67 becomes questionable. 

There is one aspect of the shock-wave heating problem we have completely neglected. This is 
the tendency of shock waves to travel along so-called "ray" paths which diverge from the vertical 
direction. Osterbrock used only a linear theory (geometrical acoustics, analogous to geometrical 
optics) for his estimates and concluded that there would be no energy left for heating the corona 
because of diffraction. Kuperus, and also Pikel'ner and Livshits (1964), pointed out that inclusion 
of the nonlinear terms would reduce this diffraction effect and thus permit coronal heating by shock 
waves generated low in the chromosphere. We are not studying coronal heating in this report; so 
we shall not consider the details of shock propagation. We consider only one -dimensional propaga- 
tion. 


Finally, there is the importance of electron thermal conduction from the corona to the upper 
chromosphere. This problem has been studied in detail by de Jager and Kuperus (1961) and Livshits 
(1964). We shall include a term for conduction in the general energy- conservation equation, de- 
rived in Chapter HI, but will not include conduction in the calculations, as its effect is demonstrably 
small except in the outermost part of the chromosphere. 
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Chapter III 

GENERAL METHOD FOR CALCULATING THE TEMPERATURE DISTRIBUTION 


It is stated in Chapter I that the aim of this study is to develop a method to solve an energy- 
conservation equation and several initially coupled transfer equations simultaneously to determine 
the temperature distribution in the atmosphere. This chapter is concerned with the development 
of this method. 

The final temperature distribution is obtained by solving the energy- conservation equation 
which we shall now derive. We begin with Equation 9, which expressed departure from RE: 


477- ^ 

N 

z: 

r i r 

dh di^ +J 

dr 

(J. - Se) ^ dz. 


j 


_ 


We desire a different form of ejqjression for the term 


(9) 



“dfT ■ 


Consider for the j line. In NETPRF, t e optical depth can be written as 


dr 

V 


4T7 






dh , 


(87) 


where the parameters are all defined in Chapter II, Section C. Ignoring stimulated emission and 
using the relations between the Einstein coefficients given by Equation 38, we obtain 


1 

d^. = - 8^? i;; ^JL ^0 "l dh . 

Now we introduce the definition 



for the j*’' line. Combining Equations 9, 88, and 89 gives 


dh 


h -LJ S-ti] 

j 



( 88 ) 


(89) 


(90) 
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where we have also used Equation 39 for (NRB). , and the approximation Equation 20, (Tex)- 

Equation 90 is the required energy-conservation equation. We must further develop the left-hand 
side, which is the radiative flux gradient that equals the mechanical dissipation in ergs cm'^ sec"*. 

First, let us exhibit a general form for the dF/dh term, which includes the effects of shock 
dissipation in the presence of magnetic fields in one dimension, and of electron thermal conduction. 
The shock dissipation term is given by Equation 79 of the previous section. The conduction term 
is given by 


F 




dh ’ 


where a good estimate for K^, according to Billings (1966), is 

- 10'® ergs cm”' sec"^ 


(91) 


(92) 


Thus, a complete expression for the mechanical dissipation term becomes 


dh 


y 

42; 


y + 1 



■ P T, 



(93) 


where dF^/dh is obtained by differentiating Equation 91. In the calculations that follow, we shall 
Ignore both magnetic field effects and conduction, for reasons given in Chapter II. Taking the 
same values of y, S, and m used to obtain Equation 81, we arrive at our final form of the energy- 
conservation equation: 


dh 


1.025 X 10® T,l 



[A. (NRB). S^j] 


477 


• dr 

(^c ” J v) dh 


dv 


(94) 


We shall show in the last part of this chapter how the right-hand side of Equation 94 depends on 
through the line source function S.^. . 

Now we can consider the radiation-transfer problem. Suppose that we ignore the last, con- 
tinuum term on the right-hand side of Equation 94. The resulting equation couples various spectral 
lines through the mechanical dissipation term which constitutes the left-hand side. We wish to 
uncouple these lines to solve N independent transfer equations for the N lines important in any 
given region of the atmosphere. Only if we know how much of the total mechanical dissipation D 
goes into exciting the upper level of the j line can we solve an independent transfer equation in 
that line. We develop an approximation for this purpose. 

Let us assume that we can estimate the fraction of the mechanical dissipation D that excites 
continuum transitions, without solving a transfer problem. This, in effect, eliminates the con- 
tinuum term on the right-hand side of Equation 94. Calling this fraction , and using subscripts -I 
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and c to denote line and continuum, respectively, we have 


D-t = (l-A,) D . 


Equation 94 assumes the form 



h "Li s^j 

j 


(95) 


(96) 


Let us further assume that we can estimate the fraction of the dissipation that excites transi- 
tions in the j line, without solving a transfer problem for that line; let us call this fraction . 
For those transitions in which the mechanical energy input equals the net radiative loss, we can 
write 


D. 

J 


Lj 


(NRB). S^. 


(97) 


Recalling modified Equation 39 for the (NRB) . , we can solve Equation 97 for to obtain a very 
useful ejqjression for the line-source function s^. in this special case: 




I. 


i 4> dv 

J X, 



(98) 


The problem is reduced to deciding if we can use Equation 98 for to solve the equation of 
transfer for the j'*’ line. In particular, is the two-level atom approximation implicit in equating 
mechanical input and net radiative loss justified? Is there a function for the ejqpression (2/A.) (D./r^j) 
that we are justified in using because it gives a good representation of (2/ A.) (d. ) as the latter 
varies with height, and that permits us to solve the equation of transfer ? We ejqject the two-level 
atom approximation to offer a good first approximation for strong chromospheric lines of the 
collision-dominated type discussed in Chapter II. We shall see in Chapter IV, where we calculate 
the height variation of (2/A.) (D./n^j) for the lines of special interest that a single negative ex- 
ponential term of the form c. exp {-/ 3 . will provide a good representation for the height range 
of interest, i.e.. 


_2 

A. 

} 





Thus, the source function of Equation 98 can be approximated by 




T <h dv + C . e 

v 1 


0 j 


(99) 


(100) 
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Solving the transfer equation with the line source function of Equation 100 is quite straight- 
forward, provided that one includes the contribution of the continuum source function to obtain a 
proper boundary condition deep in the atmosphere. A detailed treatment appears in Appendix B. 
We start with the general source function of Equation 8: 


1 + r 

V 


( 8 ) 


where = dr^/dr^ as in Equation 7, and is given by Equation 100. We assume that = b^(t^) 
and parameterize it by 


(l +A3 , (101) 

corresponding to the outwardly increasing chromospheric temperature, where 

A. and c. are constants for the j"’ line; and is the optical thickness of the region of integra- 
tion. In this study we can assume that the atmosphere is semi -infinite, i.e. , = od. Equation 101 

exhibits temperature-dependence through the exponential term 

S, , A, exp(-c, . 

1 ) ] L j 0 j J 

We shall show that this term can be ignored in the final transfer equation. Physically, this is 
because mechanical dissipation, which enters the equation through the 

C exp{-/3. 

term, predominates over the influence of the continuum for strong, collision-dominated chrom- 
ospheric lines except when r^. is very large. 

In Appendix B, we show how to transform transfer Equations into the second-order differential 
equation 


dr ^ 
0 


(Jv -sj 


( 102 ) 


where = 3 (dT^^dT ^)2 ^ by using the Eddington approximation. We can solve Equation 102 by the 
standard method of Gaussian quadrature; then the equation becomes 


dr^2 


(Ji -Si) 


(103) 
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The frequency division points "i" depend on the line profile and the order of the quadrature. For 
the Doppler -broadened line cores to which the non-LTE theory of Chapter n applies, the line pro- 
file is Gaussian and both the division points and the weighting functions x. are based on the zeros 
of Hermite polynomials. These division points and weights can be obtained from Trumpler and 
Weaver (1953). Table 7 gives x. for a three-point quadrature which we use for the calculations 
here. 

Defining such that 


s 


iz 




(104) 


we can write Equation 103 in the form 


fill 

dr^2 


J, -S. 2^ a. J.-8.Ce'^"o-S. r.Sj (l+Ae‘‘^"o) 


where we have used Equations 8, 100, and 


IT du 

J V 




(105) 


(106) 


and where the subscript j has been dropped, as we are here treating a single line. We obtain the 
a, , like the x. , once we know the line profile. They satisfy the normalization condition 



(107) 


and are tabulated for a Gaussian profile and three-point quadrature, along with the x. in Table 7. 
We have merely used the values given by Jefferies and Thomas (1960). 

We said earlier that 

» r. (Sj A) (108) 

would be satisfied for the strong lines of the report. If c ^ /3 and C » r. (s^ a) , condition 108 is 
fulfilled. That c - p follows from assuming that the region for large mechanical dissipation in a 
spectral line is roughly the same as the region where the temperature will rise owing to this dis- 
sipation, although the position of the rise might be at the bomidaries of this region, according to 
the temperature plateau arguments discussed in Chapter II, while the mechanical dissipation per 
\init mass could remain essentially constant throughout. If we examine the values of (c, sj for 
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the strong lines treated in Chapter IV, and choose the value of A to correspond to a maximum 
chromospheric of ^ 10® °K, we find that condition 108 is most nearly violated by Call. For 
Can, c = 4,2 X 10-1®, Sj = 2.5 x 10"®, and A ^ 10^. If we use = 10“®, which accords with 
Jefferies and Thomas (1960), we obtain Sj A ^ 2.5 x lO'i^. Thus, c » Sj A. However, for 
the condition becomes marginal and it is actually violated for r^. Since is most directly in- 
volved in the equation for jg (which, according to Equation 106 and the value of = 0.0051116 
plays the least important role in determining / dv) and since the problem does not even arise 
for the other lines, we accept approximation 108. Also, we may have overestimated the r. values 
for Call, Athay and Skumanich (1968) consider = 10 "i® a better choice. 

Equation 105 now becomes 


Ji 

dr^ 



a. J. - 8. Ce 

j j 1 




8 . 



(109) 


which is the form of the transfer equation we solve. This equation is solved in Appendix B, where 
the general solution to Equation 109 for a semi -infinite atmosphere is given by 


Ji 



( 110 ) 


where k. , Q. , and l . are given by 




( 111 ) 


Qi = 



( 112 ) 



0 , 


(113) 


respectively. The boundary conditions used to derive the above are j; - , as - <» and 

J. = X;-» (dj./dTp) for = 0, with the latter obtained from Krook (1955). 

Since 8; ~ i for all "i", we can make the approximation 5. = i in Equations 110, 112, and 113 
without significantly affecting the results. However, we cannot do this in Equation 111. If s. = i 
is introduced there. Equation 111 cannot be solved for n real, non-zero values for k. without 
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violating the condition 



1 . 


Consider Equation 111 further. For the three-point quadrature we adopt here. Equation 111 
becomes a cubic in making it sixth order in with k^ = ±(k.2)i/2 . We can solve this equa- 
tion for the three positive k. values appropriate to our values of (aj.x, ) already discussed, sup- 
plemented by calculated values of s. = + Alternatively, we can approximate 5. in a 

suitable manner and use values of k. already available. As this approximation illustrates the 
very small role of the continuum in the treatment of the strong spectral lines of this report, we 
discuss it now. 


The dependence of k. on r. , through S. = (i + r , can be determined from Thomas' (1960) 
work giving the dependence of the k. on parameter e, which plays the same formal role in the 
Thomas study as r. does here. For the three-point quadrature, Thomas' work shows us that k^ 
and k^ are practically independent of e and that kj ~ 0. l Thus, kj and k^ are practically 

independent of the r. values here. If we make the approximation r. - , Table 8 shows that k^ 


Table 7 

Values of x^, a. , and k. Used in Solving 
Equation of Transfer. 



Subscripts i , j for Different 
Quadrature Points. 

1 

2 

3 

X. 

1 

1.43210 

0.29078 

0.0069014 

a . 
1 

0.817657 

0.177232 

0.0051116 

k.(7. =10-®) 

0.664478 

0.045328 

0.95418 X 10-'' 


Table 8 

Values for k. Adopted from 
Jefferies and Thomas (1960). 




r . 
1 



lO"* 

1Q-® 

10"® 

>^1 

0.664575 

0.664487 

0.664478 


0.045736 

0.045365 

0.045328 

*"3 

0.94547 X 10 -3 

0.30148 X 10-3 

0.95418 X 10-" 


and k^ are indeed almost independent of 7. , and that kj o. 1 . Since r. < lO'^ for our 

lines where i = 1,2, and is probably somewhat larger for 1=3, we adopt the approximation 
Tj - = 10’® for our final calculations, and use the k, directly from Table 8 for 7^ = 10"®. 

This completes our outline of how we solve the transfer equation in the semi -infinite atmos- 
phere. In Chapter II, we stated that a rough criterion for LTE-R, where (Te), is ETo » 1. 

This is also our criterion for in applying the semi-infinite-atmosphere treatment where 

^ ('i'e (■'"o ■*“)) ^ niust be satisfied. We justify this by Equation 110 where, for k_^ /3 and 

n the order of quadrature, Ji when k_^ » i. In NETPRF, Chapter V, Thomas shows that k_^ 

scales somewhere in the range of £ through If we say k^ ~ e, then a rough criterion for 

j. - Sj becomes » i. 
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Now it is necessary to determine , the fraction of the mechanical energy of Equation 96 
that goes into the line transition. We have solved transfer Equation 102 on the assumption 
that the transfer equations for the individual spectral lines can be solved independently of one 
another. This is equivalent to saying that D. can be determined for each line. Noting that 
D. = A. defines the fractional part A. of that goes into the line, we proceed to develop 
a way to estimate A. . 

Between two levels in an atom designated L and u, the net number of upward transitions 
due to collisions per unit volume and time is 

Nu = nLCLuCNCB)^^ , (114) 


where the net collisional bracket (NCB)lu is given by Equation 41, and c^, is the collision rate, 
for which we shall use the general formula in NETPRF: 


1. 25 X 10" 






-(hO,,,/KT^ 


(115) 


Equation 115 was suggested by Burgess, Seaton, and van Regemorter, and is based largely on work 
by van Regemorter (1962). The g term is a weakly varying dimensionless parameter, whose value 
for chromospheric conditions is 2 ; 0*2 for singly ionized metal lines. 

We wish to use Equations 114 and 115 to derive an ejcpression for A. . Noting that D. = A. , 
we can use Equation 114 to write 


A. 




"Li Cluj (NCB)njj (h^-)Luj 


j N 


(116) 


2^ N„i (h^)LUi (h^)LUi 


where we have noted that 


D. = Njj. 


and 




z: 

i 


D; 
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for N spectral lines. Equation 116 is perfectly general, so far. Now, substituting the van Rege- 
morter ejq)ression of Equation 115 for in Equation 114, we derive 


A. 


'^Lj ( S ^Lu) j ^ 




IN 

^Lu)i ^ 


(h-)LUiATe 


(117) 


which is the ejqjression we use for A. and from which we obtain Dj . We shall see in Chapter IV that 
the -dependence of A. is not strong in those chromospheric regions where A. is not negligible. 

The final subject we consider here is the method for calculating the chromospheric tempera- 
ture distribution. We restrict ourselves to the problem without the continuum. In Chapter IV, we 
shall see that the H~ continuum dominates the energy balance up to some height h 2 ; 500 km, and 
we shall exhibit there a simple method for calculating the temperature distribution in these regions, 
considering the continuum alone. At present we consider Equations 96 and 97 for the lines only. 

For the one line of Equation 97, whose source function is given by Equation 98, we can easily cal- 
culate the temperature when the line source function can be written in the form 


= 


j + e. (Tj 


1 + s. 


(118) 


where s. is a general sink term equal to e + 77 , or whichever is predominant. Equation 118 is 
general enough to include Equations 31 and 35 of Chapter II, and we shall see in Chapter IV that 
Equation 118 is probably adequate for the lines we shall study in detail. 

If we equate the source fxmctions of Equations 98 and 116 and solve for (T^), we obtain 


B . (t ) 

V J \ el 



(119) 


From our solution of the transfer equation, we know 


1 , 


T 0 dv 


Thus, we can calculate the (t^^.) 's, if we know ch) , which is the same as knowing (h) . 
Therefore, to -apply the method developed here, we can assume the mass density and chemical 
composition, and either calculate or assume the ionization, in order to derive (h) . This is 
done in Chapter IV. 


Once we have (h) for a region dominated by one line, the calculation is complete when we 
have a convergent solution; but not if more than one line is important. Then, the problem of con- 
sistency between the different (h))'s and Equation 96 arises. This problem arises in the 

region 500 < h < 1250 km, in Chapter IV. 
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Chapter IV 

APPLICATION OF THE METHOD TO THE SOLAR CHROMOSPHERE 


A. Introduction 

Preliminary to our calculations, we shall review the current status of radiative -energy loss 
estimates for the chromosphere, to establish that these estimated losses do not exceed the values 
for mechanical flux given in Chapter n. Section F. We find that energy -loss estimates to date 
are unreliable by as much as an order of magnitude, but nevertheless the computed values of in 
Chapter n are adequate. 

Having established this, we shall consider the important question: wMch collisional processes 
remove mechanical energy from the system ? This can be answered by our deriving net radiative 
losses for atoms that can be approximated by only two discrete levels. For these cases, the local 
radiative loss must just balance the local gain in mechanical energy which has no other outlet. 

The question becomes more complicated if we must consider additional levels or the continuum. 
This latter case will arise with hydrogen. The result of these energy-loss estimates will be that 
the H ion continuum, Ha, the resonance doublets of CaH and MgH, Lya, and, for sufficiently high 
temperatures, Lya of Hell, should be considered as important energy sinks for the mechanical 
energy dissipated in the chromosphere. 

Given the above energy sinks, we shall proceed to the temperature calculation. We find that 
our calculation covers three chromospheric regions: one extending from h = 0 to h = 500 km, a 
second covering 500 < h < 1250 km, and a third covering the Lya region. These three regions 
corre^ond to the domination of the energy balance by the H" ion, the ionized metals and H a, and 
Lya, respectively. We shall calculate in each region, using the method of Chapter HI for the 
outer two regions, and developing a different method for the H' -dominated region. 

With regard to chemical composition and mass-density distribution we use values comparable to 
those for the PSC-model low chromosphere. PSC uses the chemical composition given in the first 
edition of Allen (1955). We shall assume the composition given in the second edition (Allen, 1962), 
which differs negligibly from values given in the first edition. The 1962 values are reproduced in 
Table 9. The mass density we assume for the region 0 5 h < 500 km is given by 

p(h) = p(0) , Hj = 130 km , (120) 

where we obtain p(0) by assmning a mean molecular weight of unity and all hydrogen is assumed 
to be neutral at h =0, so that p(0) = (0)/N^ = 6.28 x 10"® gm cm"^, where rij (O) is from 

Table 1 and is Avogadro's number. For the region 500 < h < 1250 km, we use 

p(h) = p(500) ^ = 145 km , (121) 
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Table 9 


Table 10 


Relative Abundance by Number of Chemical 
Elements from H through Fe, Based 
on Allen (1962). 


Element 

Atomic 

Number 

log 10 
Relative 
Abundance 

H 

1 

12.00 

He 

2 

11.16 

Li 

3 

3.0 

Be 

4 

2.4 

B 

5 

2.8 

C 

6 

8.48 

N 

7 

7.96 

0 

8 

8.83 

F 

9 

5.4 

Ne 

10 

8.44 

Na 

11 

6.22 

Mg 

12 

7.46 

A1 

13 

6.28 

Si 

14 

7.47 

P 

15 

5.53 

S 

16 

7.22 

Cl 

17 

5.4 

Ar 

18 

6.62 

K 

19 

4.88 

Ca 

20 

6.22 

Sc 

21 

2.91 

Ti 

22 

4.82 

V 

23 

3.78 

Cr 

24 

5.38 

Mn 

25 

5.10 

Fe 

26 

6.90 


p(h) from Equations 120 and 121 Compared 
with p(h) Obtained from PSC, 


h 

p(h) from 
120 and 121 
(gms cm~^) 

p{h) from 
PSC and 122 
(gms cm"3) 

0 

6.28 X 10"9 

6.28 X 10"S 

100 

2.92 (-9) 

2.66 (-9) 

200 

1.35 (-9) 

1.16 (-9) 

300 

6.27 (-10) 

5.23 (-10) 

400 

2.90 (-10) 

2.42 (-10) 

500 

1.34 (-10) 

1.16 (-10) 

750 

2.26 (-11) 

2.08 (-11) 

850 

1.20 (-11) 

1.11 (-11) 

950 

6.00 (-12) 

5.86 (-12) 

1000 

4.25 (-12) 

4.30 (-12) 

1100 

2.12 (-12) 

2.32 (-12) 

1200 

1.09 (-12) 

1.30 (-12) 

1250 

7.73 (-13) 

9.70 (-13) 


where p(h) is obtained from Equation 120. 
Equations 120 and 121 give values for p(h) for 
the low chromosphere that compare satisfac- 
torily (see Table 10) with those obtained from 
the PSC-model chromosphere of Table 1, 
assuming that 


p(h) - 


(h) + (h) 


( 122 ) 


For h > 1250 km, these scale heights are too low, and another representation for p(h) is neces- 
sary. We return to this when we calcxilate (h) for h > 1250 km. 


B. Dominant Energy Loss Mechanisms 

1. Estimates from Outside Sources 

In this section we first review some energy-loss estimates in the existing literature. Partly 
on the basis of this review, we can then limit the number of loss mechanisms to be studied further. 

Osterbrock (1961) provides a rough estimate of the radiative energy loss from the outer solar 
atmosphere. He assumes that begins to rise at = 0.003 and estimates that H "losses are 

completely dominant in the low chromosphere. His final value for H" loss varies linearly with the 
assumed where the temperature rise begins. Curtis (1965) places this rise at ^5000 = 0.0001; 
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this makes Osterbrock's figure appear a gross over-estimate. A large factor of uncertainty is 
characteristic of his estimates , which we shall now discuss at greater length. 

Although Osterbrock refers to his section on energy loss by the title "Observation Data," he 
makes what amounts to a theoretical estimate of the loss, which proves to be the largest. He 
derives a formula for the net H ' losses for a gray atmosphere whose temperature distribution is 
given by Equation 10, where, for the zero-th order approximation, q(r) = 2/3, By merely ejq)res- 
sing the excess of emission over absorption, with the emission produced in LTE, he derives 
Equation 123, below. His formula for the outward flux due to H" loss in the low chromosphere is 

(H‘) = 4<t ( ST) T , (123) 

where t^ is the minimum temperature, 5 t is a mean value for the temperature difference between 
the minimum temperature and the maximum temperature at the top of the region where the H ' loss 
dominates chromospheric losses, and t is the optical thickness of this region. He chooses 
Tp = 4600°K, ST = 200°K, and t = 0.003, and obtains (H") = 1.5 x 10’ ergs cm"’ sec"^. Multi- 
plying by 2, to account for energy radiated back down into the photosphere, he estimates the total 
H ■ loss from the chromosphere to be 3.0 x 10 ’ ergs cm sec " 

There are several uncertainties in t 7F^ (H"). The choice of ST and t could significantly affect 
the estimate. Also, Osterbrock did not allow for the possibility that the source function for this 
emission might depart from the Planck function owing to a departure of the b _ from unity. The 
uncertainty in r could decrease 7 tf+(H") by more than an order of magnitude. This quantity 
and ST too, both depend critically on the temperature structure of the low chromosphere and the 
extent of the low chromosphere over which H" dominates the net radiative loss. Osterbrock, in 
choosing the van de Hulst-model chromosphere, bases his estimate on a currently unacceptable 
model , as we pointed out in Section D of Chapter II. 

By choosing a radiation temperature of 6900°K, based on NRL rocket observations reported by 
Morton and Widing (1961), Osterbrock estimates the Lyman continuum loss at 3 x 10’ ergs cm"’ 
sec"', obtained by simply integrating the Planck function at T^ = 6900°K over the frequency range 
beyond the Lyman limit. In the higher-lying hydrogen continua, in which the chromosphere is 
optically thin, he computes the total radiation loss using formulas derived by Parker (1953), based 
on the work of Menzel (1937), for the emission per unit volume due to the free bound continua of 
hydrogen; thus. 


II. Ug 

i(Ba, Pa, •••,€) = 8. 94 x 10" , (124) 

where n. is the proton number density. Again using the unacceptable van de Hulst model chromo- 
sphere, and also noting that 

7TF^(Ba, Pa, • • ■ , C) = ( j dh , (125) 
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Osterbrock computes the total loss in these hydrogen continua as 5 x 10^ ergs cm ” 2 sec”i. The 
use of the PSC model chromosphere of Chapter n would increase the value of t 7F^ (Ba, • • • ) by at 
least an order of magnitude, because the PSC values of (n. ,n^) , given in Table 1 (with ~ 
assumed for h > 500 km), lie so much higher than the van de Hulst values. This conclusion follows 
even though (h) lies lower in the van de Hulst model for a given h and the latter model extends 
over a larger height range. However, there is an even more fundamental objection to using Oster- 
brock's tty ^ (Ba, • ■ • ). His value is the total loss in these continua. But most of this total loss is 
probably scattered photospheric light or radiation produced by what, ultimately, has its source in 
photospheric light that becomes redistributed in frequency. Thus, the net radiative loss to be 
balanced by mechanical dissipation is much less than if derived from Equation 125 alone. 

In treating losses in the chromospheric lines of the hydrogen Balmer and Lyman series. Call 
H and K, MgH \2795, \2803, and Sill M808, M816, Osterbrock does rely on observational data. 
These empirical estimates rely on actual intensity measurements taken of Lya emission, and 
emission in the central emission cores of the above mentioned ionized metal lines, as well as 
emission of Balmer-series lines observed in eclipse. On the basis of the above data, Osterbrock 
makes the following net loss estimates: 8 x 10® ergs cm”^ sec”^ from all Balmer lines, 8 x 10® 
ergs cm ”2 sec " ' from all metal lines, and 2x10® ergs cm sec from Lya. He examines data 
available in 1961 on helium emission and the euv spectrum in general, and concludes that helium 
and other losses are comparatively negligible. 

In assuming that the above losses are net radiative losses, Osterbrock implies that the source 
term in the source functions of these lines is the collision term «B^ . This is true for Lya 

and the ionized metal lines, but certainly not for the Balmer lines, nor for the neutral metal lines 
(which seem to have been included in his total for all metal lines, although his paper is unclear on 
this last point). Here again is significant source of error in his loss estimates. 

In all fairness to Osterbrock, he required only rough numerical estimates for his work, though 
his estimated maximum error of a factor of 3 or so in his final net chromospheric loss of 4 x 10" 
ergs cm”2 sec”^ is certainly optimistic. Also, his estimates are made for the entire chromo- 
sphere, and do not specify where the losses occur locally, which reduces their value. We have 
discussed Osterbrock’ s estimates at some length because the literature contains many references 
to them, and they are extremely uncertain. For example, see the discussion following the Light- 
hill (1967) article in the Fifth Symposium on Cosmical Gas Dynamics. 

Athay (1966) has provided the most recent and comprehensive treatment of net radiative losses 
from the chromosphere, and we shall now consider his work in some detail. He estimates the net 
H" loss from the entire chromosphere in exactly the same way as Osterbrock. Using observational 
data in choosing reasonable values for St and r in Equation 123, and performing the analysis for 
successive regions of Ah = 100 km, Athay concludes that the net H” loss lies between an order of 
magnitude below the Osterbrock estimate of 3 x 10^ ergs cm ” ^ sec ” ' and a value slightly above 
this number. Our criticism of Osterbrock' s net H” loss estimate applies equally to Athay' s 
estimates, except that Athay obtains his values from data upon which the PSC-model chromosphere 
is based. 
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Athay develops the following method for treating the losses in spectral lines. He combines 
the condition of statistical equilibrium for level L with an energy-conservation equation to obtain 
an expression for the net radiative loss in a series having lowest -level L. This is done by starting 
with transfer Equation 6 and deriving an equivalent form of oiur energy -conservation Equation 90 
for the case of a spectral line series with no continuum. His conservation equation for the radia- 
tive flux gradient, or local net radiative loss, in the series is 


dF 

dh 


k>L 


"k AuL(NRB)kL - 


(126) 


where is a suitably weighted mean energy for all the levels involved and we follow our 

regular notation. In particular, (NRB)^^ the net radiative bracket of Equation 39 and the 

Einstein coefficient for spontaneous emission. From the condition of statistical equilibrium for 
level L, the definition of the departure coefficients and the relation among the collision rates 
and given in Chapter II, Section C, for a Maxwellian electron- velocity distribution, Athay 
derives the following relation: 

^ Hk A,^(NRB),^ = . (127) 

k>L k^L i <L 

Here (NCB)j^,^ is the net collosional bracket given by Equation 44a. Now suppose that we are deal- 
ing with either a resonance series, for which A^^ = 0, or a subordinate series which is in radiative 
balance with all lower levels, making (NRB)j^; = 0. The second term on the right-hand side of 
Equation 127 is then zero. If we assume, with Athay, that these conditions apply, then substituting 
Equation 127 in Equation 126 yields 

= <Kl) 

ki‘L 


This is the equation Athay uses. 

He notes that the (NCB)j^^'s are generally much larger than the corresponding (NRB)^l'®- 
From a practical point of view, this means that the ratio \/\ need not be known as precisely 
as the term dv/s^j^j that appears in (NRB),^l> according to Equation 39, where 

~ B ) . Thus, although the b^'s are not, in general, independent of the radiation field, an 

estimate of the (NCB)^^ based on rough solutions of the transfer equation or even cruder proced- 
ures will not be as imreliable as an estimate of the corresponding (NRB)^^^ if the latter is obtained 
in an equally rough fashion. He makes use of this property of the (NCB)j^^ to estimate losses in 
various spectral series, using Equation 128 and rough values for the ratios obtained from 

various sources. 
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Table 11 


Collisional Excitation and Ionization Rates for Hydrogen with 
= 10^0 Computed by Henze and Reported by Athay (1965b) 


Transition 


T 



1 X 10 * 

2.5 X lO"* 

4 X 10“ 

8 X 10“ 

1-2 

0.0050 

5.2 

30.0 

1.3 X 10 2 

1-3 

7.0 X 10"5 

0.27 

2.2 

13.0 

1-4 

1.0 X 10"5 

0.063 

0.57 

3.7 

1-5 

3.3 X 10"® 

0.025 

0.24 

1.6 

1 - k 

6.1 X 10"® 

0.18 

2.6 

28.0 

2-3 

5.4 X 10 2 

2.5 X 1Q3 

3.8 X 103 

5.3 X 103 

2-4 

18.0 

1.5 X 10^ 

2.7 X 102 

4.5 X 10 2 

2-5 

3.3 

35.0 

67.0 

120.0 

2 - k 

21.0 

3.8 X 102 

8.1 X 102 

1.6 X 103 



For the Lyman and Balmer series of atomic hydrogen, Athay (1965b) uses the results of some 
collisional excitation-rate calculations by Henze, who in turn used Seaton's (1962) weak-coupling 
method to evaluate collision rates C^u for several transitions arising on the ground or the first 
few excited levels of hydrogen. Some of these results are given in Table 11. Judged by the table, 
the following estimates by Athay seem not xmreasonable: 


c. , » ; c 


z: 


12 -Ik 

k>2 


(129) 


and: 


^23 » ^2k 

k>3 


(130) 


where inequality 130 is restricted to ^ 25000°K. Using inequalities 129 and 130, Equation 128 
becomes, for the Lyman series. 


- n, C,,(NCB)j2 ’ (131) 

and for the Balmer series, 

"2C„(NCB)23 (hv)^3 , (132) 

subject to the condition that the (NCB)j^y's of Equations 131 and 132 are sufficiently predominant. 
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Since (NCB)^^^ = (l , Equation 131 is valid provided that b^ » b^, which condition cer- 

tainly prevails for the upper chromosphere, where Athay assumes the bulk of the net chromospheric 
loss occurs in the Lyman series. Though he does not justify Equation 131 on theoretical grounds by 
estimating terms like (NCB)j 3 (hi^j^gat lower heights in the chromosphere where bj ~ b^ > bj 

prevails and where is considerably larger than in the upper chromosphere, nevertheless the 
observed predominance of Lya over the other members of the Lyman series gives some support 
to his estimate ejqiressed by inequality 129. It is evident from the relevant b^ 's that (NCB) ^3 
= (1 -b 3 /b 2 ) is the dominant (NCBj^u for the Balmer series. We shall discuss the variation of 
these b^'s in the chromosphere presently. Equations 131 and 132 state that the net radiative 
losses in the Lyman and Balmer series of hydrogen are dominated by the corresponding collisional 
terms for Lya and Ha, respectively. Whether the net radiative loss actually occurs in that spec- 
tral line cannot be determined until a transfer problem is solved. 

Athay's conclusion that, for h > 1000 km, the condition bj » b^ exists, is based on his analysis 
(Athay, 1965a) of strong resonance lines in which the effect of continuum processes on the net flux 
in the line is ignored. He shows that for those lines where the further condition e » -q B*/B^ (t^) 
is met, one obtains 



(NRB),j • 


(133) 


Using a result due to Avrett and Hummer (1964), which is a refinement of Thomas' (1960) earlier 
rough estimates of the variation of (NRB) ^ with t^, and which says that, for Lya, 

(NRB)^. ~ V (134) 


where 


Tg ~ 0. 1 T? 1 


Athay concludes that, at heights above -0 - 3x 103 in Lya, bj » bj is valid. This follows from 
Equations 133 and 134 where one takes v = 3 x 10 "3, although the "»" condition becomes marginal 
below Tg 2 ^ 10 2 . Also, since Athay discusses Lya losses beginning as low as h = 1000 km, using 
the result bj » b^, which gives (NCB)^^ 1 in Equation 131, he is introducing a contradiction 

with the PSC model chromosphere of Table 1 which shows > 7 x lO^ in Lya for all h < 1300 km. 
Thus, it is debatable whether his treatment of the Lyman series losses is valid for h < 1300 km. 
Essentially the same conclusion can be drawn from a study of b^, • • • , b 3 by Pottasch and Thomas 
(1959), These authors evaluate the b^'s under the assumption of detailed balance in the Lyman 
lines, but variable opacity in the Lyman continuum. Using = 10“* °K, = 10**, and = 6000°K 

for the Balmer continuum, they obtain a significant increase in b^ over b^ in the range 10 > 

> 1. If we examine (h) in the PSC model chromosphere, we find that 10 > > 1 in the 

region around h = 1300 km. 
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Thus, Athay's application of Equation 131 for Lya with (NRB)^^ ~ 1 is somewhat unsatis- 
factory. However, he does not rely on Equation 131 alone for net loss estimates in Lya, but 
observes that Hinteregger's (1964) measured flux in Lya of 2.7 x iQs ergs cm'^ sec certainly 
makes Lya an important source of net energy loss from the chromosphere. Unfortunately, these 
measurements do not tell us where in the chromosphere the losses occur. 

To treat the Balmer series we need and b^. Here, Athay states that (NCB) 23 - 1 for the 
region of interest (500 < h < 1250 km), which implies that b^ » b^ in this region. This is in 
complete accord with the Pottasch and Thomas study already mentioned. Athay then estimates 
Balmer series net losses with 


-357- = ^^2 ^23 

which is essentially Equation 132 with (NCB )23 = 1. Using the afore-mentioned due to Henze 
and taking values for from the PSC-model chromosphere, Athay derives the values for 

the net Balmer series loss— see Table 12. He estimates = 4 x lO''^ ergs in Equation 135, 

intermediate between Ha and the Balmer continuum. 

Table 12 includes Athay's report of the 
total radiative output in the Balmer continuum 
observed during the 1952 Khartoum eclipse 
expedition which provided the data for the 
PSC-model chromosphere. A comparison of 
these values with the calculated values of 
(dFg^/dh) Ah reveals that the calculated net loss 
in the Balmer lines series is roughly half the 
total radiation observed in the Balmer con- 
tinuum. Athay also observes that the total 
radiation observed in Ha at Khartoum is 
4x10® ergs cm’^ sec which is more than 
than twice the calculated net loss for the 
Balmer series and approximately equals the 
total radiation observed in the Balmer con- 
tinuum. Examination of these statements reveals that much of the observed radiation in Ha and 
the Balmer continuum must be scattered photospheric light, which is what we would expect for 
radiative transitions with the source function controlled by photoionization and not by collisions. 

Athay ignores helium losses, and in particular Lya of HeH at ^304, claiming that the observed 
energy radiated in these lines is only about 5 percent of the observed energy radiated in hydrogen 
Lya. Since the hydrogen Lya value represents, effectively, a net loss, and since the various 
helium lines’ observed loss must be equal to or greater than the net loss from helium, he con- 
cludes, reasonably enough, that helium losses are negligible. 


Table 12 

Computed (dFgjj/dh)Ah and Observed Total 
(Not Net) Losses in Balmer Continuum, 
Both from Athay (1966). 


Ah 

1 dh / 

(ergs cm ^ sec ') 

Total Observed 
Balmer Continuum 
Losses 

500- 700 

6 X 10® 

1.3 X 10® 

700- 900 

4 X 10® 

8.3 X 10® 

900-1100 

4 X 10® 

5.7 X 10® 

1100-1300 

3 X 10® 

3.6 X 10® 

1300-1500 

1 X 10® 

2.1 X 10® 
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Finally, Athay ignores the neutral and ionized metal lines, claiming that the former are com- 
paratively weak while the latter have (NCB) « 1.0 for a lesser optical depth than for the 
hydrogen lines. His conjecture for the (NCB)j^jj for the ionized metal lines is rather speculative 
and, for this and other reasons we shall discuss presently, we do not accept his decision to ignore 
the ionized metal lines. 

Zirker (1960) also derived an energy-conservation equation from the equation of radiative 
transfer, for the purpose of estimating net radiative losses from the solar chromosphere. He 
obtains the radiation -field terms for spectral lines and continua at a given frequency, and j^, 
respectively, by solving two coupled transfer equations, using the Eddington approximation. He 
then uses these values of and in the energy-conservation equation to obtain the dF/dh for 
those spectral regimes treated. His results include a quantitative treatment of helium, so we 
include them here. 

He estimates net radiative energy losses from material at the center of a uniform hydrogen- 
helium atmosphere satisfying the following conditions: thickness = 10^ km, total particle density 
= 10 '1 cm "3, hydrogen -helium abundance ratio = 10, and Balmer continuum intensity fixed by a 
background "photosphere" which produces black-body continuum radiation with a radiation tem- 
perature of 6000°K. These conditions roughly approximate upper chromospheric conditions of the 
PSC-model chromosphere, as given in Chapter II, Section D. 

Zirker's calculations give the net radiative loss per gram at the center of the above atmos- 
phere as a function of electron temperature for the following radiators: the Lyman and Balmer 
continua and Lya for hydrogen, the resonance continua of Hel and Hell, the Balmer continuum of 
Hell, and Lya of HeH. His results show the Lyman continuum of hydrogen predominant for 7500 
< < 25000°K, Lya predominant for 25000 < < 50000°K with the Lyman continuum still non- 

negligible and Lya of Hell predominant for > 50000°K. The predominance of the Lyman con- 
tinuum of hydrogen for < 25000°K is contrary to the observationally based conclusion of 
Osterbrock that the Lyman continuum is negligible. Since Zirker's study cannot completely 
duplicate the actual chromosphere, where net losses in the Lyman continuum and Lya depend 
strongly on the actual chromo^heric opacity variation for these transitions, we do not believe 
that the suggested predominance of the Lyman continuum for < 25000°K applies to the 
chromosphere. 

A brief digression from chromospheric loss estimates is now in order, to consider the net 
radiative losses from the corona. We must include these losses in our discussion, as sufficient 
mechanical energy must be available to balance them. In particular, if this mechanical energy 
has its source in the turbulent motions near and above the top of the hydrogen convection zone, 
as we are assuming for the mechanical energy that heats the chromosphere, then enough energy 
must be left over after the shock waves have passed through the chromosphere to balance the 
coronal losses. We have already mentioned two difficulties which the accretion hypothesis (heat- 
ing from an external, mechanical source) encounters when applied to the corona in Chapter H, 
Section F. It is very likely, therefore, that the coronal net radiative losses must be added to 
those of the chromosphere for an estimate of the total outer atmospheric net radiative losses to 
be balanced by the mechanical source in the lower photosphere. 
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Billings (1966) estimates net coronal energy losses from radiative, conductive, and convection 
mechanisms. The latter must be added, when they are non-negligible, if they represent net losses 
from the outer atmosphere. Billings concludes that the net radiative losses from the corona 
exceed net losses from convection (solar wind) and outward conduction by about one order of 
magnitude. He bases his net radiative -loss estimates on a review of several observational studies 
of the solar ultraviolet radiation. These studies treat emission lines whose intensities lie so far 
above the background continuum that they effectively yield net radiation losses. Lya losses are 
not included in his estimate, as Billings accepts the hypothesis we have already reviewed that Lya 
losses occur in the upper chromosphere. He concludes that the net radiative losses from the 
corona lie in the range 10® - 10 6 ergs cm'^ sec'L 

There is one other important source of energy loss from the corona which, according to 
Billings' estimates, is also of the order of 10® ergs cm sec "L This is the electron thermal 
conduction down to the upper chromosphere which we discussed briefly in Chapter H. As this does 
not constitute a net radiative loss from the corona, it is omitted from estimates of the same. In 
fact, one current conjecture (de Jager and Kuperus, 1961) is that the upper chromosphere is heated 
largely by conduction down from the corona, producing the intense Lya emission observed. We 
do not include any heating of the Lya region due to conduction, but it could be very important, 
Kuperus and Athay (1967) argue that thermal conduction down from the corona might even be the 
source of energy for the spicules. 

We devote no further attention to the corona here, as we shall not carry our temperature 
calculations out beyond the upper chromosphere. Theoretical estimates have been made of 
coronal radiative energy losses, assuming an optically thin corona. Pottasch (1963, 1964) has 
estimated losses due to oxygen, carbon, silicon and magnesium, concluding that for > 25000°K 
these elements produce radiative losses exceeding those of hydrogen and helium. By neglecting 
dielectronic recombination (which causes a higher fraction of the net radiative loss to come out 
of the corona in the form of emission lines rather than free-bound continua) it is possible that 
Pottasch may have erred in estimating at what frequencies the corona will lose energy. For our 
purposes, we shall not rely on theoretical estimates of coronal radiative losses, but accept as an 
upper limit Billings' higher value of 10® ergs cm“^ sec'L 

It is now possible for us to offer a rough estimate of the net radiative energy loss from the 
entire outer solar atmosphere. We do this to see if the estimates of the mechanical flux given 
in Chapter H, Section F, provide enough energy to balance these net radiative losses. We see from 
an examination of the estimates of Osterbrock and Athay that there is much uncertainty regarding 
the lower limit of the net radiative loss from the chromosphere, because of the scattered-light 
problem, but that we can accept their higher estimates to provide a number that should serve well 
as an upper bound. Thus, the chromospheric net radiative losses should not exceed 4x10’ ergs 
cm sec'i. Billings' estimate for the upper limit to the coronal net radiative losses, 10® ergs 
cm" 2 sec"i , does not alter the 4 x 10’ ergs cm"^ sec"^ value when applied to the entire outer 
atmosphere. Thus, our criterion for accepting the proposed turbulent generation of acoustic waves 
is that the resulting mechanical flux satisfy F^ ^ 4 x 10’ ergs cm"’ sec"’. In Chapter H, 
Section F, we showed that F^ is probably in the range 2.9 x 10® ^ F^ 2.7 x 10’ ergs cm"’ sec"’. 
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The problem we consider for, the^remainder of this section is the choice of spectral lines and 
continua to include in our final temperature calculations. We shall base this choice partly on the 
work of Osterbrock, Athay, and Zirker already discussed, although this work is not always ade- 
quate for our needs. 


2. First-Hand Estimates 

The most straightforward procedure for estimating net radiative losses in a collision- 
dominated line is to solve transfer Equation 6 for the case where the line source fimction is 
given by the general form 


, (136) 

where the sink term is s = e + t). This is the same as our general source function of Equation 25 
except for the neglect of source term tjb*. Temperature dependence enters primarily through the 
Planck function (T^) and, to a lesser degree, through the sink term s. NETPRF and the source 
function series of papers by Thomas and Jefferies give numerous examples of the solution of the 
transfer equation with a source fimction of the form of Equation 136, subject to the approximations 
of a constant sink term s and a single exponential representation for the Planck function. Given 
the solution for / <p^ dv. Equation 136 is then used to get s^. Finally, as in Thomas (1960), the 

(NRB)^l of Equation 39 is evaluated and the net radiative loss is obtained from 

dF 

dh = "u Aul (NRB)ul , (137) 

provided that is known. Having for the line and the population for the lower level, we 
can obtain n^, from the relation 


= B^(n„/nJ . 


So one procedure for determining the predominant energy-loss mechanism is to evaluate Equation 
137 for all suspected lines and retain those lines for which dr/dh predominates in some region of 
the atmosphere. 

The procedure we use for treating the lines and most continua here is simpler than solving 
the transfer equation. We rely on the A. of Equation 117, derived in Chapter HI; A. is the frac- 
tional amoimt of energy that goes into the transition L - U by inelastic collisions and does not 
necessarily give the energy lost from the atmosphere in a corresponding U - L radiative transi- 
tion. This occurs only in the restricted case of a two-level atom. Otherwise, the net radiative 
loss must be obtained by solving a transfer problem. For example, when we conclude that Ha is 
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a transition for which A. should be evaluated for the low chromosphere, we must remember that 
the mechanical energy going into the 2-3 transition in hydrogen will not, in general, equal the 
net energy loss in Ha. In relying on A. , we replace the problem of estimating net radiative loss 
mechanisms by the problem of estimating net collisional gain mechanisms for levels (L, U). This 
is the spirit of the Athay (1966) study, reviewed in the last subsection. Our only departure from 
this procedure is in the treatment of H", where the small optical thickness of the chromosphere 
in the continuum permits a simple estimate of the net radiative losses. 

We begin with the lowest part of the chromosphere, where H" is important. We proceed along 
two lines. We estimate the excess of emission over absorption in the lowest region 0 < h < 500 
km, including the non-LTE effect due to the departure of the continuum source function from the 
Planck fimction. We also estimate the mechanical energy dissipated in this region, according to 
our assumptions about shock dissipation of acoustic waves, to see if the mechanical energy dissi- 
pated roughly equals the energy loss due to H" and any other significant sources of radiative loss. 

To estimate the excess of emission over absorption in H", we first calculate the b _ whose 

H 

value determines the departure of the source function from the Planck function. If we assume 
that H“ dominates the continuum opacity and ignore the free-free absorption, it is easy to show, 
as in PSC, that the continuum source fimction is given by 


Sc = —bJtJ . 

H 


(138) 


The dominance of the continuum opacity by H" for 0 < h < 500 km is justified presently. That 
we can neglect the free-free absorption, compared to the bound-free absorption, is apparent from 
the estimates of a.t k - 4700 A: 


T, (”K) 

5040 

6300 

^ V f b 
J vf f 

24.8 

10.1 


We obtain these estimates from the following equations, which PSC derives using the Planck 
function and the Saha equation: 


4. 14 X lO'i® 


2nhv~ 


"e "hT, 


X” - 

-3/2 ^ H 


V b f 


(139) 


= 4, 14 X 10' 


S'nViV 
3c ^ 


n n„T e -'■"a „ 

e H e IV f f 


(140) 


The units for are cm^; for they are cm'* dyne Using the from Chandrasekhar 

(1958) and the from Chandrasekhar and Breen (1946) gives the we exhibit. Newer 
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calculations of the make our conclusions stronger. 4700 A is a representative line-free 
point of the solar spectrum near the maximum observed intensity. For z 16000 A, where 
does become large, 



5800°) dv 


is very small compared to 



so our neglect of the free-free absorption is justified. 

Induced emission, also, can be ignored in our further work with H“, owing to the fact that 
(hv) > 0.75 ev, the ionization potential for H‘, and hence the factor exp{-hv/kT^} is negligible in 
the chromosphere. 

With these preliminaries, we can go on to calculate b _. From the condition of a statistically 

H 

steady state, we obtain 


n. D 


_ + n D _ b + 

e , H " H , H "2 


I, 


" b,(tJ 


hv vbf 


a... . dv 


n. D 


_ + n„ D 
e , H " H , H 


+ 2n 


c ^ 


" J.(T.ff) 


hi^ 


“vbf 


(141) 


where we use the notation of PSC. The units for the collision rates per particle are cm^ sec~*. 
Thomas has shown that the term for electron collisional excitation n d _ is orders of magnitude 

° e , H 

less than the term for the associative detachment reaction n„ D _ under photoimheric conditions. 

H , H 

This conclusion remains valid in the low chromosphere, when we use more recent values for the 
electron-negative hydrogen ion cross section (lower values) reported in Tisone (1966), or even 
the older, higher cross section used by Thomas in making his estimates of n d _ for the 

e . H 

photosphere. Frisch (1966) argues in her thesis that the concentration of is maintained in LTE, 
giving b =1. Thus, Equation 141 becomes 


n„ D + 4 tt 

" H, H 




hv <^vbf 


n„ D _ + 2 t7 
” H, H 


•'v , 


J.(Teff) 


hv v\> f 


a . , dv 


(142) 
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Table 13 


Values for Departure Coefficients b^_ in Lowest Region of the Chromosphere. 


-r 


•'-0 

- 


0. 95 X 10® sec ^ . 


h (km) 

rijj (cm" 3) 

“ H. H 

(sec“^) 

T, (°K) 

r b„(tj 
•'*'0 

(sec" 1) 

b 

H 

0 

3.78 (101 S) 

5.70 (106) 

5050 

1.34 (10®) 

1.06 

100 

1.60 (15) 

2.08 (6) 

5250 

1.50 (6) 

1.18 

200 

6.99 (14) 

0.95 (6) 

5460 

1.66 (6) 

1.37 

300 

3.15 (14) 

0.41 (6) 

5680 

1.84 (6) 

1.65 

400 

1.46 (14) 

0.19 (6) 

5910 

2.02 (6) 

1.94 

500 

6.95 (13) 

0.09 (6) 

6150 

2.24 (6) 

2.24 


This equation gives us b We exhibit the results of our calculation in Table 13 above. To derive 

H 

b , we have used the PSC model for the low chromosphere for , the already quoted value 

of = 1.3 X 10 “9 cm3 sec"i (Dalgarno, 1967; Schmeltekopf, et al. , 1967), and the value of the 

integral 


{‘■el -X - , , 

- 1.5x lO^T^e H , (143) 

given in PSC, which uses the Chandrasekhar (1958) values for For ^ , we choose 

= 5800°. 

Table 13 shows that at h = 200 km the associative detachment reaction rate n„ d has 

H H H 

decreased to become equal to the photoionization rate 



2rr 



a 


vh i 




Beyond 200 km, the photoionization rate dominates. These conclusions are similar to Schmelte- 
kopf's conclusions, which we gave in Chapter II, Section F. 

In the low chromosphere, we would ejqiect a combination of the photoionization and our as- 
sumed local value of mechanical energy input to fix the local electron temperature until a 
height is reached where becomes a minor source of continuum opacity. We use this to calculate 
the temperature in the lowest region 0 < h < 500 km in the next section. 
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With these b we can now write the excess, E, of emission over absorption in H : 

H 


E 



S^K^di^ - 477 



J,. K di7 . 

V V 


( 144 ) 


If we use Equation 138 for and ignore the frequency-dependence of the nearly grey H bound- 
free absorption coefficient, Equation 144 becomes 

E 2aK^’_ [2T,‘* - (b^.)T/„] , (145) 

where cr is again the Stefan-Boltzmann constant and K *_ = n *_ a _ for the LTE population of n 

_ H H H H 

We obtain K *_ (h) by noting that, at h = 0, H completely dominates the opacity and b _ ^1, so 

H H 

n„ n * . From the Bilderberg Continuum Atmosphere or the opacity tables of Cox and Stewart 

H 

(1964a, b) for physical conditions at h = 0 in the Sun, we find K •_ ~ 10"^ P ~ 6.31 x 10"^^ cm“^. 
We use the h = 0 value of p of the PSC model to obtain the latter value. In addition, the PSC model 
gives us values of n (h). Thus we obtain 


K *. (h) 

H 


K *. (O') 

H 


n *_ (h-) 

H 

n *- (O') 

H 


(146) 


The results of the calculation of excess emission E appear in Table 14. 

It is interesting to compare the results of this calculation with calculated values of the shock 
dissipation D(h) in the same chromospheric region. Table 15 gives D(h) obtained by using (h) 
from the PSC model, p(h) from Equation 120, from Equation 84, and Equation 81 for D(h) itself. 
We see that D(h) ^ E(h) at h = 400 - 500 km and D(h) ~ E(h) ath = 300 km, all of which is con- 
sistent with the idea of H" dominating the energy loss out to h ~ 500 km— though we must point out 
that this is no proof, considering that every estimate upon which this coincidence is based is quite 
rough, particularly the estimate of D(h) . Furthermore, the huge excess of E(h) over D(h) as we 
proceed in toward the limb is not unreasonable when we consider that the (h') values from the 
PSC model chromosphere are almost certainly too high there. We conclude that the mechanical 
energy input and the excess emission in H ' are roughly equal in the region 0 £ h < 500 km. The 
total mechanical dissipation / D(h) dh in this range is approximately 10® ergs cm“^ sec“^. If we 
ignore our very large estimates of E(h) below 200 km, the J E(h') dh over this region yields a like 
value. This is roughly equal to Athay's lower estimate and is a good order of magnitude below 
Osterbrock's estimate. It suggests that the net radiative loss from the entire outer atmosphere 
may lie between 10® and 10’ ergs cm”^ sec“* , well below the oft-quoted Osterbrock value. 

We must consider the other possible energy_ sinks for the mechanical energy in the region 
0 < h <■ 500 km, to check our conjecture that H' does dominate the net radiative losses there. 

We must consider the hydrogen continua and the spectral lines from hydrogen and the neutral 
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Table 14 


Values of Excess Emission E(h) in H" from Low Chromospheric Region 

0 < h < 500 Ian. 


h (km) 

K * (cm"i) 

H 

T^ (°K) 

b 

H 

E(ergs cm "3 sec-i) 

0 

6.31 X 10" 

5050 

1.06 

0.714 

100 

1.87 (-11) 

5250 

1.18 

0.401 

200 

6.66 (-12) 

5460 

1.37 

0.173 

300 

2.64 (-12) 

5680 

1.65 

0.063 

400 

1.09 (-12) 

5910 

1.94 

0,031 

500 

4.71 (-13) 

6150 

2.24 

0.018 


Table 15 

Calculated Shock Dissipation D(h) for Height Range 
0 5 h < 500 km Using D(h) = 1.025 x 10® . 


h (km) 

P (gm cm " ®) 

T, CK) 

V, 

D(ergs cm ® sec *) 

0 

6.3 X 10"’ 

5050 

0.0 

0.0 

100 

2.9 (-9) 

5250 

0.22 

0.0158 

200 

1.4 (-9) 

5460 

0.39 

0.0452 

300 

6.3 (-10) 

5680 

0.53 

0.0548 

400 

2.9 (-10) 

5910 

0.63 

0.0440 

500 

1.3 (-10) 

6150 

0.71 

0,0303 


metals. The temperatures below 500 km are too low for any helium transitions to be important 
and we shall show, in our treatment of the ionized metals to follow presently, that the optical depth 
in the lines of the more abundant ionized metals becomes so large below 500 km that detailed 
balance can be ejqiected for them there. 

First, consider the hydrogen continua. It will be convenient to consider these continua for the 
entire low chromosphere, rather than restrict ourselves to the region 0 < h < 500 km. 

In the lower part of the low chromosphere, where o f h < 500 km, the departure coefficients 
will all lie near unity, with the extreme variation from unity given by bj = 3.4, according to 
the PSC model of Table 1. Actually, may be slightly greater than bj here, but the difference 
is not important for our purposes. Since the b,^ for k > 2 all lie below this value, and since all 
the b|^ decrease rapidly as height h decreases, we expect all the (NCB)^,^ connecting the boimd 
levels of hydrogen with the continuum to rapidly approach zero with decreasing height in this 
region, at least xmtil such height is reached where bj b^ decreases below unity. 

Frisch (1966) has calculated bj for this low region, assuming bj = b^ and including the con- 
tribution of H" and the Paschen continuum to the opacity and an Osterbrock (1961) type of 
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mechanical dissipation term. She gives values for ranging from 0.78 to 1.00 at h = 0, depending 
on initial conditions for the model. These are minimum values for in the region 0 < h < 500 
km; the thickness of the sub-region for which bj < 1 never exceeds about 100 km in her models, 
except in one extreme case where it is about 200 km. These b^ < 1 will make the (NCB)^^ negative 
for b,^ = 1. However, these negative values seem to occur over a thin sub-region only, the very 
sub-region near our h = 0 point where we know that all our estimates— D(h), E(h), and (h) —are 
rather unreliable. Thus, we ignore the departures of (NCB)^,^ from zero in this sub-region, and 
exclude the hydrogen continua from our energy balance in the lower region, 0 < h < 500 km. 


In the next region, 500 < h < 1250 km, we are faced with a different situation. Here the 
(NCB)^^’s for the lower levels of hydrogen are better approximated by unity than by zero. It be- 
comes necessary to calculate the quantities k = i, 2, • • • up to some 

value of k where this quantity becomes negligibly small. We do this for the first three hydrogen 
continua, using (NCB)^^ = 1, the from the previously referenced Henze work for = lO^, 

= lO'*, Table 1 for values of and , the Boltzman equation and b^ 0. 2 b, for Hj, and the 
previously referenced Thomas and Pottasch study and Tables 1 and 16 for the behavior of b^, b^ , 
and bj as functions of and of height h. Table 16 gives the results of this calculation. 


Table 16 


Calculated Values of (NCB)^,^ (hv)^,^(ergs cm‘3 sec'i) 

for the First Three Hydrogen Continua for Selected Heights 
in the Range 500 < h < 1250 km. 


h (1cm) 

(ergs cm sec'i) 

C2^(NCB)2^ 

(ergs cm'^ sec“ i) 

"3 ^3,<(NCB)3,,(hv)3^ 
(ergs cm“3 sec“^) 

750 

1.6 X 10"2 

1.1 X 10"3 

2.4 X 10"3 

1000 

3.0 X 10-3 

7.3 X 10-‘> 

1.6 X 10'^ 

1250 



4.2 X lO”'* 

4.9 X 10"* 

1.0 X 10"3 


Table 16 shows that the amount of mechanical energy delivered to the chromosphere in the 
first three hydrogen continua is not completely negligible compared to the mechanical energy input 
given by J D(h) dh over the region 500 < h < 1250 km. We have already noted that this integral 
has a value of '^lO® ergs cm'^ sec“‘. We readily see that the mechanical energy required by each 
of the first three hydrogen continua for 500 < h < 1250 km is of the order of 10 s ergs cm’^ sec“i. 
We have somewhat overestimated the requirements of these three continua by using = 10‘*, 
making c^,^ too large. However, this could be balanced by the additional requirements of the higher 
continua, which must be included tmtil (NCB)^,^ becomes a small fraction of unity. The hydrogen 
continua could easily take up half of the mechanical energy jD(h) dh which we estimate to be 
available. 


The Balmer continuum deserves further treatment. The observed radiation in the Balmer 
continuum is far from negligible, as examination of Table 12 reveals. Since Table 12 is confined 
to the region 500 f h < 1500 km, we should ask if the Balmer continuum losses below 500 km are 
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of comparable magnitude. Athay (1966) remarks that the observed radiation in the Balmer con- 
tinuum actually reaches a maximum at 400 km (below which good data were not available from the 
1952 eclipse), so this observed radiation is of the order of the mechanical energy dissipated in the 
lower region 0 < h < 500 km, according to Table 15 for D(h). The comparatively low values of the 
mechanical energy estimated in Table 16 to ionize hydrogen from the second bound level suggest 
that much of this observed Balmer continuum radiation is scattered photospheric light. However, 
this still leaves the possibility of more absorption of Balmer continuum radiation than is emitted 
in the chromosphere, so we must assure ourselves that there is no appreciable net transfer of this 
energy into the thermal modes. 

Thomas in PSC evaluated the net radiative (or rate) bracket for the Balmer continuum (NRBC) ^ 
for the two cases of detailed balance in the Lyman continuum and negligible Lyman continuous 
opacity, where the atmosphere is assumed transparent in the Balmer continuum but is locally 
opaque in the Lyman lines. In both cases, he finds (NRBC) ^ equals a negative number. Physically, 
this is due to the predominance of photoionization from level 2 over collisional excitation between 
levels 2 and 3, an effect we mentioned when discussing the source function for photoionization - 
dominated lines. The negative value of (NRBC )2 shows that there is no net radiative loss from the 
atmosphere in the Balmer continuum, but rather that the atmosphere actually does absorb more 
radiant energy in the Balmer continuum than it emits. 

We now ask if this could be an important source of thermal energy. We proceed by determining 
whether the net flux of photospheric radiant energy into the Balmer continuum is small compared 
to the local mechanical dissipation, and thus ignorable. We start by exhibiting several useful rela- 
tions derived by Thomas in PSC. For the Balmer continuum, it is easy to show from the equations 
in PSC, Chapter IV, that the flux gradient under conditions of detailed balance in the Lyman lines 
can be approximated by 


- 8.2 X 10'^ (NRBC) 2 , (147) 

where a three-level atom plus continuum approximation has been used and no transfer problem has 
been solved. Under the same conditions, using values of departure coefficients bj and for the 
low chromosphere taken from Pottasch and Thomas (1959), and also assuming transparency in the 
Balmer continuum, Thomas obtains the two expressions for (NRBC ) 2 mentioned before. These are 

(NRBC)^ = -A . (148) 

for detailed balance in the Lyman continuum, and 

(NRBC>2 = -(2 +A) , (149) 

for transparency in the Lyman continuum. The quantity A is given by 

A ~ 0.4 63 , (150) 

with all parameters defined in Chapter II, Section C. 
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If we take t 

e 

150, we obtain 


= 6000°K and further reduce Equation 


49.0 b (NEB) 


6000 °K 


(151) 


Since the temperature dependence of A is not too strong for 
5000 < T^- < 8000°K, we need only determine bg and(NRB )3 2 
for Ha, and we can evaluate the flux gradient in the Balmer 
continuum in the low chromosphere with Equations 147-149, 
using Equation 151 for A, 

Consider , first the lower region 0 < h < 500 km. Here, 
^ 100, as shown in Table 17, which also shows that this 
lower region satisfies detailed balance in the Lyman con- 


tinuum. 


Using the result b^ < b. 


to obtain b, ~ 1 for 


Table 17 


Height Dependence of 
and from PSC Model 
Chromosphere . 


100, and using Thomas' (1960) values for (NRB)^ 


OHa 

as a fimction of 


h (km) 

'^UyC 

^OHa 

750 

1 X 103 

34.5 

850 


28.5 

950 


23,4 

1000 

2 X 102 

21.7 

1100 

96.0 

18.7 

1200 

35.0 

17,0 

1250 

19.8 


1300 

6.6 





100, (NEB), 


= 0 . 0012 ), we obtain the value A ~0.06. Using this value and combining Equations 148 and 147, we 
obtain 


~ - 3.0 X 10"“' ergs cm"^ sec”^ , (152) 

0 £ h A 500 km . 

Comparing this with the values for E(h) for H" and D(h) for the local mechanical energy input 
(Tables 14 and 15, respectively), we see that dFg^,j,/dh is more than two orders of magnitude below 
either and can surely be neglected in the region 0 < h < 500 km. 

Now consider the higher region 500 £ h < 1250 km. Repeating the type of calculation done for 
the lower region, we find that 


dF 


Ba 1 C 


dh 


10 ^ ergs cm ^ sec 


(153) 


500 £ h £ 1250 km , 

Vi^here we use intermediate values of (NRB )32 as functions of from Thomas (1960), Equation 
151 with 2: 6000°K for A, bj 1 from Pottasch and Thomas (1959), and either Equation 148 or 
149 for (NRBC)^. It does not matter much which equation, 148 or 149, is used as Equation 151 
gives A ~ -1, so (NRBC)^ is of order unity in either case. Equation 153 is uncertain to more than 
a factor of two anyway, as we are in the region of intermediate Lyman continuous opacity where 
all quantities are estimated at values between those for optically thick and thin atmospheres. 
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We are now faced with two problems; 1) the roughness of the above analysis, difficult to im- 
prove without solving a transfer problem, and 2) the magnitude of dFg^,^/dh, which says that the 
net flux of photospheric radiant energy into the Balmer continuum is of the same order as the 
mechanical energy dissipated in this region. We see the latter by noting that /D(h) dh for 500 < h 
< 1250 yields ^lO® ergs cm'^ sec“i, while f (dFg^j^/dh)dh also gives ~106 ergs cm'^ sec'i in 
this region. 

Without going into detail on where this dFg^jj,/dh energy flux goes, we can make one simple 
observation concerning energy balance. The region currently under discussion has ^ Terr 
according to the PSC model of Table 1, and the photospheric radiation field absorbed by the chro- 
mosphere is well described by a Planck function at T 2 ; in the continuum. Thus, the energy 
of the Balmer continuum photons absorbed in this j region cannot heat the gas beyond 

T^f f without violating the second law of thermodynamics. As far as competing with D(h) as a 
source of thermal energy, the Balmer continuum can be ignored in this upper region, provided 
that the condition > T^f j is actually satisfied there, as the PSC model suggests. 

This completes our special study of the Balmer continuum. Let us consider the lower region 
0 < h £ 500 km again, regarding the mechanical-energy input into bound-bound transitions of 
hydrogen and the neutral metals. 

First consider the hydrogen transitions. The mechanical energy input into Lya is negligible 
because of the great optical thickness of the lower region in the Lyman continuum, so that , 

giving (NCB)j 2 ~ 0- However, the Ly/3 transition may take up much mechanical energy. If we use 
the Table 11 value for at = 10 and = 10 and the (h = 200 km) from Table 1, and 
allow that b^ at h = 200 km may be as large as 2.0 with bj £ 1, we see that rij (NCB)j 3 
is more than an order of magnitude larger than the total D(h) at 200 km. On the other hand, 
should really be recalculated in this region, since the actual temperature is more like 5000° than 
10,000°. Also, Frisch’s (1966) models show that bj (h £ 200 km) may be as low as unity, depending 
on the amount of mechanical energy dissipated in the low region which is poorly known. Thus, it 
is quite possible that our above estimate at h = 200 km is several orders of magnitude too large. 
Below and above h 2 ; 200 km, c^g (NCB)jg (hi^)jg should drop off from the value around h = 200 
km because b^ < 1 below and (h) decreases above. Finally, C^^is almost an order of magni- 
tude less than Cjg so the members of the Lyman series above Ly^ are comparatively unimportant 
mechanical-energy sinks. Unfortunately, we cannot be sure from these estimates that Ly/3 is com- 
pletely negligible in the region 0 < h < 500 km. 

The local net gain of mechanical energy by Hais C^g (NCB)gg (hi^) ^g. Even for (NCB)^ = 1, 
which is too large for the low region, we see that c_, (NCB)„, (hv),, << lO"^ ergs cm“^ sec'^ , 
considering our use of Cgg at = 10''. Thus, the Ha collisional excitations are small enough 
compared to D(h) to be ignored. The same argument shows all of the Balmer and higher hydrogen- 
series transitions to be negligible compared to D(h). 

To estimate the mechanical-energy input into transitions of the neutral metals, we cannot 
assume LTE to be valid in the low chromosphere, as Zirker (1956, 1958) has shown in his analysis 
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of several lines of neutral metals (Na, Mg, Al, Ca, and Fe) for which he provides estimates of the 
height h at which (h) equals unity for each line. The necessity for some non-LTE effect arises 

when we compare 1) the inference from the Saha equation that the ionization of the metals should 
leave such a small LTE ratio of neutrals to ions (lO”^ - 10“^ for 0 < h < 500 km) that the chro- 
mosphere should be optically thin in the metal lines, with 2) Zirker's analysis which shows definite 
self -absorption effects for these lines. Zirker uses data from the same 1952 Khartoum eclipse 
that the PSC-model chromosphere is based on. Unfortunately, the data were of poor quality for 
h < 500 km, and many of his results do not directly affect the region 0 < h < 500 km, where the 
neutral metals could play an important role. 

The rough estimates of radiative losses by Osterbrock and Athay, giving the total observed 
radiative flux in the lines of all neutral and singly ionized metals as a few times 10® ergs cm” 2 
sec”* , are not too helpful either, particularly when one considers that each is based on different 
assumptions, different data, and even different spectral regimes. Furthermore, the observed 
fluxes from the neutral metal lines are certainly not net radiative losses from the chromosphere, 
but arise primarily from scattered photospheric light, if we judge by their similarity to the Balmer 
lines, i.e., most have photoionization-dominated source terms. On the other hand, to examine the 
100 strongest, or 1100 strongest such lines described by Athay, even by several "prototypes" of 
similar excitation potential, arising from atoms of similar fractional abundance, etc., would be 
extremely time-consuming. And it should be pointed out that the total observed loss from the 1100 
mentioned was observed to be twice that of the strongest 100, so we cannot restrict ourselves to 
the few strongest when discussing these lines. 

We proceed, therefore, by first noting that the total radiative flux from the chromosphere in 
the neutral lines is not likely to exceed 10® ergs cm"^ sec”^, if we judge by the observations of 
Osterbrock and Athay. Then we recall that lines with photoionization-dominated source functions 
scatter more energy than is received from collisional excitation, as mentioned above. How much 
more, we could estimate by examining the predominance of the photoionization source term -qB* 
over the collision source term in the line (Te ) • However, since the mechanical dissipation 
D(h) of Table 15 jdelds a value for /D(h) dh over 0 £ h < 500 km of order 10® ergs cm' 2 sec"* , 
it follows that at least the bulk of our mechanical energy must be available for collisional excita- 
tions in other energy sinks than all the neutral metals taken together. 

This concludes our consideration of the other possible mechanical-energy sinks in the region 
0 £ h < 500 km. We cannot discount the possibility that the cumulative contribution of the hydro- 
gen continua, the hydrogen series, and the neutral metals could equal /D(h) dh over this region. 

Our earlier conjecture that only H” is an important mechanical energy sink in the lower region is 
weakened by the above considerations. 

Nevertheless, we shall ignore the contribution of all mechanical energy sinks other than H” in 
this lower region, justifying this procedure on the generally reasonable correspondence between 
E(h) and D(h) in Tables 14 and 15, at least above 200 km. However, we must keep in mind the 
reservations mentioned here. 
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This concludes our consideration of the region 0 < h < 500 km. What about the chromosphere 
for h > 500 km? Provided that does not exceed 50000°K, we can ignore helium, according to 
the Zirker (1960) study, so we consider the lines of neutral hydrogen and the abundant singly 
ionized metals. 

First, consider the hydrogen lines. We have already noted that, for ^ 25000°K and chro- 
mospheric conditions at h > 500 km, the condition of Equation 130 

^23 ^ ' ^2k ’ 

k>3 


holds, so in the upper region of the low chromosphere the 2-3 transition is the only Balmer 
collisional transition that could be a significant energy sink. The Paschen and higher hydrogen 
series can be ruled out because of the \'s. For k > 3, if ^ , (NCB)^ 2: ^.nd, if b^ 

» + j (j = 1. ••■) , then the number density . In either case, (NCB);^ 

(hv)i^ k + j negligible when compared to this quantity for the Balmer case with (k, j ) = (2, 1). 
To see this, we use again the table of Henze's provided by Athay (1965b), and the non-LTE 

Boltzmann relation (with the b^ ’s), and assume that the behavior of b^ (^Lyc) all k > 3 mimics 
that of b 3 from Pottascfi and Thomas (1959), subject to the condition b^ i b^,.j . This covers the 
case of hydrogen except for the Lyman series, which we now examine. 

Consider once more inequality 


c 


1 2 



k >2 


which says that the 1-2 collisional transition is the dominant energy sink in the Lyman series as 
long as (NCB) ^ Where in the chromosphere does this mechanical energy input occur? To 
answer this question, we compare the quantity (NCBj^^ (ht') for this transition with the 

analogous ejqjression for the singly ionized metals, which are treated next. When this comparison 
is made for the region of h < 1250 km, we find the condition (NCB)jj 2 ; 0 to hold, by examining the 
PSC-model low chromosphere and the Pottasch and Thomas (1959) study of departure coefficients 
bj and bj . The PSC model gives ^ 8000°K and r(LyC) ^ 10 for h < 1250 km. The Pottasch 
and Thomas study gives bj < b^ for r(LyC) > 10, 6000°K, = lO^, = 10“*, and detailed 

balance in the Lyman lines. Thus, (NCB)^^ 2: ^ should be a good assumption below 1250 km and 
above the H~ region. From this it follows that ® below 1250 km, so we 

ignore this transition there. 

Consider the singly ionized metals. Because of their low ionization potentials, we would 
ejq)ect and shall assume the more abimdant metals to be singly ionized in the region 500 < h 
< 1250 km. Zirker (1956) has studied this problem. From his determination of optical depths in 
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the resonance lines from the singly ionized metals BaH, Sell, Srll, Till, observed in the 1952 
eclipse, he calculates the neutral hydrogen density at h = 500 km, assuming that all these metals 
are in the singly ionized state above this height (ionization potentials: Ba, 5,21 ev; Sc, 6.54 ev; Sr, 
5.69 ev; Ti, 6.82 ev) and that the Allen (1955) relative abundances apply everywhere. He obtains 
excellent agreement with the value for the neutral hydrogen density reported in Table 1 for the 
PSC model, which is based on continuum data. This supports the assumption that above 500 km, 
all metals with I.P. 8.0 (Ca, Mg, Si) will be singly ionized in the region 500 < h < 1250 km. Ad- 
ditional support for this assumption is provided in a study by Athay and Zirker (1962), which shows 
that for = lO^i, = 10“*, most of the calcium will be singly ionized, but that, by the time 
= 15000°K is reached, most of the calcium will be doubly ionized. They obtain these results by 
approximating the Call ion by three botmd levels (4S, 3D, 4P) and the continuum. Their results 
are independent of optical depth in the center of the resonance doublet. In view of the above result, 
when we examine the temperature structure of the PSC-model chromosphere of Chapter n we can 
reasonably assume that all calcium is singly ionized in the region 500 < h < 1250 km and, for 
h > 1500 km, calcium is at least doubly ionized. We can make no statement for 1250 < h < 1500 
km, except that the transition from predominantly singly to doubly ionized calcium occurs there. 

What transitions in the singly ionized metals should we e3q>ect to be strongest in the solar 
spectrum? From considerations of chemical abundance and observations of strong chromospheric 
lines, we conclude that the resonance transition in Mgll is by far the strongest, and that the 
resonance transitions in Call and Sill might also be important. This conclusion is supported by 
the spectra of Purcell, et al. (1963) for Mgll, Dunkelman and Scolnik (1959) for Call, and Detwiler, 
et al. (1961) for SHI, These three studies span the spectral range 4000A - 500A and show that the 
resonance lines of the above three ions and Lyaof hydrogen dominate the line spectrum. Also, in 
the extreme ultraviolet spectra down to 250A reported by Zirin, et al. (1963) and Austin, et al. 
(1966), there are emission lines of highly ionized elements characteristic of the corona, but no 
line comparable in intensity to those just mentioned. 

The disk spectra of Mgll, Call, and Sill exhibit the emission cores characteristic of the singly 
ionized metals and illustrated, with a central reversal, in Chapter H, Section C, The central 
reversal is also present. Mgll has a particularly strong emission core. 

In addition to assuming that all magnesium, calcium, and silicon are singly ionized in the 
region 500 < h < 1250 km, we can show that it is reasonable to assume that the ground state 
population can be approximated by the total population of the ion. Since < 10,000° in this region, 
according to Table 1, we consider the most unfavorable case (t^ = 10,000°). Using partition func- 
tions from Allen (1962) and estimating the fraction of the total population lying in the ground state 
for each ion, we find that this fraction is greater than one-half for all three ions. Table 18 shows 
this. Non-LTE effects favoring the gro\md will enhance this result. Table 18 is derived from 
the relation 


- h. 

n “ U ’ 


(154) 
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Table 18 


Ratio n^/n of Ground State to Total Ion 
Populations Using Partition Functions 
at T = 10000° from Allen (1962). 


where and u are the statistical weight of 
the ground state and the partition function of 
the ion, respectively. Of course, = (2S^ + i) 
^ (2Lp + l) = 2(Xp + l) for the one-electron 
states. 

Now let us consider which of the impor- 
tant mechanical energy sinks discussed for 
500 < h < 1250 km should be included in the 
final temperature calculation for this region. 
We shall not treat the hydrogen continua or 
the collisional transition for Ly^S further, even 
though they may absorb half the mechanical 

energy. Thus, an uncertainty of a factor of 2 is introduced into the amount of mechanical energy 
D(h) available for the other four. This uncertainty does not affect the final results to a first ap- 
proximation. Table 19 contains the necessary atomic parameters for the ions we consider. The 
first three columns need no ejqjlanation, except that refers to the principle line of the multiplet. 
in column 4 is the Einstein coefficient for spontaneous emission, calculated from 


Ion 

Ground 

State 

Partition 

Function 

U 

no/n 

Mgll 

3S 

2.04 

0.980 

Can 

4S 

3.48 

0.575 

Sill 

3P 

6.0 

1.0 


(\il) t 


0.667 


10 ® 


(gf)T 


(155) 


Table 19 


Atomic Parameters for Selected Ions. See Text for Description. 

(Lya) = ^ ® • 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 







Relative 




Ion 

Transition 

Terms 

(A) 

(\l) T 

(sec“^) 

A = i;: A,, 

(fLu)T 

Element 

Abundance 

(log 10 

“o 

(cm 2) 

To0^> 

e 







scale) 




Hydrogen 

2g . 2p0 









(for Ha) 

2p0 . 2g 

6562 

0.441 X 10® 

0.428 

0.6408 

12.00 




2p0 _ 2p 









Call 

2g . 2p0 

3933 

1.51 X 10® 

0.700 

1.05 

6.22 

3. 1 X lO'io 

Tgj X 10"“- 

2. 5x lO"'* 

T 1 / 2 
A e 

MgH 

2g . 2p0 

2795 

2.56 xlO® 

0.600 

0.90 

7.46 

1.45x 10" 

Tjij X 10"®- 

9. 2x 10"® 

J 1 /2 

Sin 

2p0 _ 2p 

1816 

8.10 xlO® 

0.445 

0.67 

7.47 

7. 5 X 10"ii 

X 10"®- 

2.5x 10"® 

X 1/2 
^ e 
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where (gf is the product for the multiplet, g^ is the statistical weight of the 

upper level, and is the wavelength in microns (Allen, 1962). Subscript T refers to the entire 
multiplet. The statistical weight g^ for the k-th level is given by g^ = ( 2 S^ + i) (x,, + l) for the 
k-th level in Mgll, Call, and Sill and by g^ = 2 k ^ for the k-th level in hydrogen. Column 5 with 
A = (Su/Sl) Aul^o^ self-ejqjlanatory. In column 6 f^u is obtained from Allen using f^u = (gf)^/gL . 
Column 7 follows from Table 9 of relative abundances of the elements. 

Column 8 is the theoretical absorption coefficient at line center for zero radiative damping 
given by 


4 X 2\ ^ 

^gj_ ^JL -^0 y 877^/2 


(156) 


where m is the mass of the absorbing ion. Equation 156 is derived in Allen (1963). It reduces 
further to 


a 


0 


= 4.06 X 10'® 


A 



(157) 


by our definition of A = (su/El) (a„l • We use to estimate the optical thickness of the chro- 
mosphere in the ionized metal lines. By assuming that the relative chemical abundance is constant 
with height, and recalling our assumption that all magnesium, calcium, and silicon are in the 
ground state of the singly ionized ion for 500 < h < 1250 km, we can use 


0 j 
"Oi 


(158) 


to estimate the ratio of optical thicknesses in different lines. Equation 158 follows immediately 
from the variation of with . Subscripts "i" and "j" refer to different lines. If 

we know eh) for one line, we can find it for the others. 

If we examine Table 1 to estimate (h) for Lya, we find that a reasonable fit to the few data 
points is obtained by the equation 


r^j(Lyct) = 10^3 X 10'® (3 X 10® - h(cn))] (I59) 

This is illustrated in Table 20. The fit is adequate for our purposes in the region 500 < h < 1250 
km, though we are not justified in using Equation 159 for h > 1250 km, where we shall actually 
treat Lya itself. The scale height corresponding to Equation 159 is 145 km, the value we chose 
for ^(h) in the region 500 < h < 1250 km in Section A of this chapter. 
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Table 20 


Comparison of '^'^Lya ^ from PSC Model 

Chromosphere with Our 

T (h) = lO^^XlO”® (3XlO®-h)] 

from Equation 159. 


h (km) 

PSC T (Lya) 

Equation 159 

t-q 

750 

1 X 10 7 

5.6 X 10® 

1000 

2 X 10® 

1 X 10® 

1200 

3 X 10 = 

2.5 X 10® 

1250 

2 X 10® 

1.8 X 10® 

1300 

7 X lO'' 

1.3 X 10® 


For the number density (h) , which 
appears in Equation 158 for , we use 


-(h-500)/H. 

"Lj ~ "i (500) ~ ^ 


= 145 km , (160) 


where (500) = 8.07 x 10^3 corresponding to 
Equation 120 for no hydrogen ionization and 
mean molecular weight unity in the region 
0 f h < 500 km, and (n^j/oj) is obtained 
from the relative abundances of Table 10. 


We should mention that there is some inconsistency between our adopted scale height of 145 
km for the region 500 < h < 1250 km and Zirker's (1956) value of 200 km at a height of 500 km, 
slowly increasing to a value just under 300 km at a height of 1500 km. (This latter value is highly 
uncertain due to scatter in the data.) Zirker's results are based on his study of lines from singly 
ionized metals. Even if we assume that his results give Hj = 200 km throughout the region 500 
< h < 1250 km, we note that this gives densities at 1250 km which are 6 times as great as those 
obtained with = 145 km in Equation 160. We use the smaller scale height for greater consis- 
tency with the PSC model of the low chromosphere, noting that our previous conclusions drawn 
from Zirker's work concerning single-stage ionization of the metals do not strongly depend on 
scale height. 


In Table 19, the last column gives values of e computed from 


£ 


21 


V’ 


g 1 


- h v / k 




(161) 


which is given in NETPRF and is derived by combining the definition of e in Equation 27, 
from Equation 115, from Equation 155, and Equation 43 connecting and for a Max- 
wellian electron-velocity distribution. We use = IQi', = 10'', and g = 0.2 to derive e for 
singly ionized metals. 

To estimate the fractional part of the mechanical energy going into the j-th energy sink (which 
is equivalent to estimating the quantity A. (h) for the four lines we consider here), substitute 
Equation 160 for rij^.(h)into Equation 117 for A. (h) , use g = 0.2, and approximate the net- 
collisional-bracket by (NCB)j^yj = 1 to derive the equation to be solved for A. (h) , which is 


(h) 


LUj 


-(hv) 

e 


LU 


/kT^ 



-(hi/) 

e 


, ■ /kT„ 


(162) 
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The g = 0,2 value is very good for the three singly ionized metals and is somewhat too large 
for Ha, introducing a temperature -dependent error of about a factor of 2 which will give over- 
estimation of A. (h) for Ha (van Regemorter, 1962). Since the values of A. (h) for Ha are less 
than or equal to the values for the other principle mechanical-energy sinks, this error is not 
serious. 


Consider the assumption (NCB)j^y = 1, For Ha, this demands that bj « b^. Although bj 
< b^ for 500 < h < 1250, we offer no proof that bj « b^ for this region. The Pottasch and Thomas 
(1959) study does not apply in this region as it was done for = lO'*. Athay (1966) reports that 
(NCB)j 2 -= 1 is justified here, and offers some estimates of bj made by Henze which support this 
view. We assume that (NCB )32 = 1 in our calculation of A. (h)'s. 


Regarding the singly ionized metals, Athay (1966) notes that the product C^u Mgll 

resonance doublet in particular may exceed the same product for Ha by as much as two orders of 
magnitude at h = 500 km. He also notes that higher opacities in singly ionized metals at a given 
height h should make the b^^ 's of the upper and lower state approach each other higher in the at- 
mosphere, which is reasonable. Thus, allowing that (NCB)^^^ = 1 where is sufficiently small for 
all these lines, we would indeed expect (NCB)^^^ to approach zero (or depart from unity) higher in 
the atmosphere for these ionized metal lines than for Ha. How then do we justify assuming (NCB)j^jj 
= 1 for ionized metals over 500 < h < 1250 km without a detailed separate study of these lines, 
when we have every reason to believe that this may be completely wrong for the lower part of this 
region? We do so because our results obtained in the next section with this assumption tell us 
that, for any reasonable value of the mechanical-energy input obtained from the weak-shock theory, 
(h) in this region lies well below the PSC values, and that letting (NCB)^^ approach zero mid- 
way through the region will not raise the temperature. Thus, as we shall argue presently, either 
more mechanical energy is being dissipated in this region than the weak shock theory gives, or the 
metals are not mostly ionized in this region, or severe departures from spherical symmetry are 
encountered in this region. Some combination 

of these effects might also be encountered. The Table 21 

point is that this conclusion is independent of . . „ 

Values of Fractional Apportionment Coefficient A. 

how the (NCB) varies in the region 500 < h fQj. Predominant Lines in the Region 

< 1250 km. 500 < h < 1250 km. (h) is from Table 1, 


So we perform the calculation of the A. ' s 
using Equation 162 and the temperature distri- 
bution of the PSC model of Table 1, and exhibit 
the results in Table 21. These are interesting. 
The temperature dependence of the A. 's for 
Call and Mgll for the temperature range con- 
sidered is rather weak, with variation of less 
than a factor of 2 for Call and effectively no 
variation for Mgll. Thus the entry of into 
Equation 162 does not affect the grosser results 
of Table 21, except for the case of Ha. The two 


h 

T 

e 

^Ha 

^Call 

^Me 1 1 

i I I 

500 

6150 

0.0154 

0.4389 

0.5417 

0.0040 

750 

6650 

0.0483 

0.3851 

0.5602 

0.0064 

850 

6840 

0.0705 

0.3573 

0.5652 

0.0070 

950 

7040 

0.1017 

0.3311 

0.5590 

0.0082 

1000 

7150 

0.1200 

0.3195 

0.5522 

0.0083 

1100 

7400 

0.1660 

0.2861 

0.5385 

0.0094 

1200 

7700 

0.2218 

0.2502 

0.5178 

0.0102 

1250 

7901 

0.2547 

0.2420 

0.4922 

0.0111 
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dominant metal lines also greatly dominate Sill, which is not smprising at these temperatures. 
Finally, the lower region, where the assumption (NCB)j^^ = 1 is most in error for the metal lines, 
is also the region where is smallest, reducing the possibility that the energy delivered to the 
two dominant metal lines is a small fraction here. 

We conclude that Mgll, Call, and Ha are important energy sinks for the region 500 < h < 1250 
km. Sin may not be negligible at some height above 1250 km but it certainly is below 1250 km. 

The final region where we calculate is h > 1250 km. Here, Lya of hydrogen is important. 
We showed earlier that we could exclude Lya from the region of the singly ionized metals where 
the PSC model says ^ 8000° and T(LyC) ^ 10 for h < 1250 km. At some higher h, where bj > 
so that 0 < (NCB) < 1, the product (NCB)^^ (hz^)i 2 becomes predominant. If we use 

n, = IQi", n = lO'i T = 10'», and c. = 5 X 10“ 2 from the Henze table for these conditions, we 

l''e*^e'12 ^ 

obtain (NCB)^^ ergs cm'^ sec'S which is even larger than D(h) in this upper 

region, according to our weak-shock theory. However, uncertainties in our , n^, and 
(NCB)^j in the upper region could account for this. We merely note the support given our as- 
sumption that Lya of hydrogen should dominate the energy balance in the region immediately 
above, where the singly ionized metals are important. 

Our further assumption that we may ignore Lya of Hell, \304, is supported by the calculations 
of Zirker (1960) already mentioned, using the admittedly questionable PSC-model upper chromo- 
sphere of Table 4 for which (h) < 50000°. Further support is provided by the disk euv observa- 
tions reported in Athay (1966) showing Lya of Hell to be negligible. On these grounds, we include 
only the Lya line of hydrogen in our temperature calculations for the upper region h > 1250 km. 
Our results support the assumption that Lya is formed in a region for which < 50,000°. 

C. The Chromospheric Temperature Calculation 

Our calculation of the temperature distribution is divided into three parts, corresponding to 
three different height regimes in the chromosphere. In what we call the "low chromosphere," we 
make separate calculations for two subregions: 0 < h < 500 km, which we call the H“ region, and 
500 < h < 1250, where we calculate the temperature distribution from the Mgll and Call resonance 
doublets. The third region for which we calculate the temperature distribution lies at heights above 
1250 km. We extend this last calculation up to small t^'s in Lya. 

1. Temperature Distribution from in 0 < h < 500 km 

This calculation employs the values of b _ obtained in the previous section (Table 13) and is 

H 

therefore subject to the assumptions upon which these b _ values are based. These assumptions 
are 1) that we can ignore the free-free H“ opacity and induced emission; 2) b^ 1; 3) the PSC 
model chromosphere of Table 1 gives the best available representation of (n^,T^) for this region; 

3) D _ = 1.3 X 10"9 cm^ sec"' is the best available value for the associative detachment 
reaction-rate constant; and 4) Equation 138, which states that = ^b (T^.), is a good repre- 

sentation of the continuum source function. We have already discussed assumptions 1,2, and 3. 
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This leaves assumption 4, which is rigorously true only if H” is the dominant source of 
opacity. We have already noted that this may not be true for h -?,400 km, because of the Paschen 
continuous opacity. In particular, Frisch (1967) claims that above 400 km the Paschen continuous 
opacity is predominant; she offers this as the reason why the Cayrel mechanism is not effective 
in heating the low region above 400 km, since for bg = 1 the of Equation 163 (below) approaches 
the Planck fimction (t^) as the E 3 term becomes predominant. 

To determine the effect of the Paschen continuum on our continuum source fimction s , we 
note that, with both H" free-bound and Paschen continuous opacity considered, can be written 


s 


B.(Te) 


b 

H 


1 + Ro 


‘^3^3. fb 


(163) 


where R, ,, = (n * a, /n * a \ . If we assume that b = 1 for 0 < h < 500 km, and evaluate 

3 , fb V 3 3 , fb' jj- 3 - - » 

R 3 from PSC, which gives 


4 X 10*5 T3/2 


' 3 , fb - 


(164) 


we can evaluate the quantity (l +R 3 fb)y/j^b^_ +R 3 = S^/B^ (T^) to determine the effect of the 

Paschen continuous opacity on the source function in this region. Using the PSC model of Table 1 
for once again, and b _ from Table 13, we compute R 3 jj, and (l +R 3 +^ 3 , fb) > 

giving the results in Table 22. For comparison we include ^b -)”*> which equals S^/B^ (T^) when 
R 3 equals zero. 

We see immediately that, for the PSC model, we can conclude that H" continues to dominate 
the opacity out to 500 km, and that the continuum source function is well represented by 
Sc = (b^-)-‘B^(Tj • 

Table 22 


Values of R 3 and Related Quantities to Illustrate 
Dependence on Paschen Continuous Opacity. 


h (km) 

10 * 3 n 

e 

(cm" 3) 

T, (“K) 


b 

H 

1+R3.fb 

b + R,hf 

H" 


0 

6.36 

5050 

3.32 X lO""* 

1.06 

0.943 

0.943 

100 

5.04 

5250 

1.37 (-3) 

1.18 

0.848 

0.847 

200 

4.65 

5460 

4.72 (-3) 

1.37 

0.731 

0.730 

300 

4.58 

5680 

1.47 (-2) 

1.65 

0.610 

0.606 

400 

4.63 

5910 

4.29 (-2) 

1.94 

0.526 

0.515 

500 

4.72 

6150 

0.119 

2.24 

0.474 

0.446 
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We now use the b _ (h) of Table 13 along with the approximation = (b 
culate (h) in the region 0 < h < 500 km. Our conservation equation for this region becomes 

d(tJ = 2o-k;_ ; (165) 

this is exactly like Equation 145, which we used 
to estimate the excess of emission over ab- 
sorption in H", except that here we identify the 
left-hand side with the local mechanical-energy 
input, d(t^) is always given by Equation 84. 
Using the same K *_ and b _ as in this last 

H H 

section, and T^ff = 5800 again, we calculate 
the (h) given in Table 23, which includes 
(h) from the Table 1 PSC model for com- 
parison. Commensurate with the imcertainties 
in d(t^) , the agreement is quite good. 

This temperature calculation offers only a 
consistency check on the temperature distribu- 
tion of the PSC model and does not prove this the best model available for the low chromosphere in 
the light of current observations. Nevertheless, in spite of our disbelief in a boundary temperature 
as high as 5000°, this is as reasonable a distribution as we know to date for 0 < h < 500 km, and it 
leaves unchanged the PSC conclusion that there is an appreciable rise of the order of 1000° in the 
first 500 km. 

2. Temperature Distribution from Call and Mgll in 500 < h < 1250 km 

Here, we apply the method of Chapter III to calculate the temperature distribution in this upper 
portion of the lower chromosphere. We have already shown that Mgll and Call are the predominant 
mechanical-energy sinks among the singly ionized metals studied in Section B of this chapter. Sill 
was shown to absorb <1.0 percent of the mechanical energy as long as < 8000°. Ha was shown 
to be an important mechanical-energy sink, but because the source function for Ha is of the photo- 
ionization type discussed in Chapter II, we shall restrict the temperature calculation to the reso- 
nance doublets of Mgll and Call, for which we assume a two-level atom source function of the form 
of Equation 30, 

= TT-l . (30) 

for which the solution to the transfer equation is particularly simple. 


Table 23 

Values of (h) for 
0 < h < 500 km. 


h (km) 

K (cm- 1 ) 

H 

b 

H 

T^ (“K) 

(PSC)(“K) 

0 

6.31 X 10“" 

1.06 

4980 

5050 

100 

1.87 (-11) 

1.18 

5090 

5250 

200 

6.66 (-12) 

1.37 

5320 

5460 

300 

2.64 (-12) 

1.65 

5655 

5680 

400 

1.09 (-12) 

1.94 

5985 

5910 

500 

4.71 (-13) 

2.24 

6285 

6150 
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Table 24 


Calculated Values of Mechanical Dissipation Term (2/Aj )(Dj/rij^p for Call 
and Mgn with Local Values of Selected Representations C. e ^ ° ^ 
from Equations 166 and 167. 


h (km) 


A. n, . 

- 1 i"!. 

Call 

500 

0.000 


750 

1.91 X 10-1“ 

850 

2.55 X 10-1“ 

950 

3.17 X 10-1“ 

1000 

3.40 X 10-1“ 

1100 

3.71 X 10-1“ 

1200 

3.83 X 10-1“ 

1250 

4.19 X 10-1“ 





1.71 X 10-*2 
1.57 X 10->o 
2.56 X 10'‘® 
3.28 X 10-10 
3.53 X 10-10 
3.85 X 10-10 
4.01 X 10-10 
4.07 X 10-10 


0.000 

1.86 X 10-11 

2.70 X 10-11 
3.59 X 10-11 
3.94 X 10-11 
4.68 X 10-11 
5.30 X 10-11 

5.70 X 10"ii 


7.50 X 10-1 = 

1.16 X 10“' 1 

2.57 X 10-11 
3.83 X 10-11 

4.29 X 10" 11 
4.94 X 10-11 

5.30 X 10- 11 
5.42 X 10- 11 


As a preliminary to the calculations, we substantiate the claim, made in Chapter HI, that we 
can represent the line source fimction ‘ Sf in the form 




1 


J <* di/ + Ce 

j 1/ 




( 100 ) 


for both the Mgll and Call resonance doublets. In terms of the parameters of Section B of this 
chapter, this means that Equation 99, 


C. 

J 



(99) 


must be fulfilled, where the "j" refers to a particular line. We can easily show that this represen- 
tation is possible. Table 24 compares (2/A.)(D./n, .) with c exp (-/?. for Mgll and Call. The 

A. are obtained from Table 19. The D. are obtained from Equation 81 and the A. of Table 21. Fi- 
nally, the n, . are obtained from the assumptions of constant chemical composition with height, us- 
ing the relative chemical abundances by number in Table 9, and a density scale height of 145 km, 
where we assume all hydrogen to be neutral at h = 500 km and all magnesium and calcium to be in 
the ground state of the first ionization stage throughout 500 < h < 1250 km, to obtain absolute values 
of n^^.. The values of and /3. used for the Call and Mgll resonance doublets, respectively, are 

For Call For Mgll 
C 4.2x10-1“ 5.7x10-11 

/3. 10-2-1 lO'S-'i 
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Substituted in Ce ° , these give 



Call 


4. 2 X 


10- 10 


^-(10 2- 7)^jj(CaII) 


(166) 


ajid 



Mgl I 


5. 7 X 


10" 11 


-(10 ‘‘)to(MeII ) 
e 


(167) 


Table 24 amply justifies our choice of (c. ,>SJ for each line. Equations 166 and 167 are now used 
to calculate the temperature distribution. 

First we solve the transfer equation for each line in the manner outlined in Chapter HI, to 
derive 




a. j. 
1 1 


for each line. We use Equations 166 and 167 in the source function, which is of the form 100 in the 
transfer equation. We shall also solve the transfer equation for a second set of c. e for each 

line, for reasons that appear soon. This second set is 



"] 


Ca I I 


4. 2 X 


10"5 


^-(10 ^l)Tg(Call) 


(168) 


and 



Mgl I 


5.7 X 


10"! 1 


^-(10 “')Tp(Mgll) 


(169) 


The actual solution of the transfer equation follows the method outlined in Chapter III and 
developed in detail in Appendix C. The results are 



, -8- . 286t^ 

2. 497 x 10"® - 1.822 X 10 o + ... 


2.830 X 10 


-7-.019ST- „ ^„-6- . 000868T. 

® - 0.8772 X 10 0 


- 1.2929 X 10 for Call, corresponding to Equation 166 ; (170) 
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2.513 X 10'^ - 7.895 X 10 + 


- 1.187X lo"’ _ 3,3529 X lo‘® 


+ 3.2396 X 10 ^ . 00004 iTg Mgll, corresponding to Equation 167 


(171) 




2.497 X 10"® - 1.660 x 10 ’’ •^®®'^o 


^ ,«"6'.0195r« ^ ,^-6- . 000868TA 

- 2. 586 X 10 ° - 8.772 x 10 ° 


+ 9.263 X 10 * .000041 g Call, corresponding to Equation 168 


(172) 


and 




2.513 X 10"^ - 8.098 x + 


- 1.164X 10’®'-°'®^^° - 3,3530 X 


^ 3.4629 X 10 ® oooo41tq Mgll, corresponding to Equation 169 


(173) 


To obtain a temperature distribution from these solutions, rewrite Equation 119 


b.(tJ 


1 + S -/3 

Ce 


+f 


djy 


where the sink term s = e + 77 becomes s 2 ; £ for the Call and Mgll resonance doublets. To obtain 
Tg (h), we use the (h) relations for these doublets derived in Section B and given in Table 19. 
The results of the temperature calculation are given in Table 25. 

To interpret Table 25, note that T„ (l) and T. (l) both correspond to the set of 

Call Hgll 

parameters represented by Equations 166 and 167, i.e. , the initial set obtained from all the as- 
sumptions and estimates made so far. On the other hand, the T„ (2) and (2) corre- 

'Call ®MgII 

spend to the set of parameters represented by Equations 168 and 169, which differ from the first 
set in that the coefficient c has been increased by one order of magnitude. 
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Table 25 


Results of (h) Calculation for 500" < h < 1250 Icm Based on Call and Mgll 
Resonance Doublets. T (1) and T (1) Correspond to Equations 

^Call 

170 and 171. T (2) and T (2) Correspond to 

Call ®MgII 

Equations 172 and 173. 


h (km) 

^Ocall 

T, (1) 

® C a I I 

T, (2) 

Call 


Te (1) 

® Mg I I 

Te (2) 

® M g I I 

600 

1.38 X 103 

4360 

5783 

1.12 X W* 

5412 

7029 

750 

4.91 X 10^ 

4413 

6024 

3.99 X 1Q3 

5624 

7476 

1000 

87.3 

4512 

6279 

7.10x 102 

5364 

7038 

1250 

15.5 

4501 

6271 

1.78 X 102 

5187 

6745 


We first note that the temperatures obtained from the first set of parameters are too low, 
particularly for Call. Even if our value of 6285° at 500 km, obtained from the H ' calculations, is 
somewhat high, it is not reasonable, in the light of everything we have done so far, for at that 
height to be much below 6000°K. Retaining our spherically symmetric model for the moment, we 
see three possibilities; 1) The parameter e may have been overestimated. Equation 119 with 
s = e « 1, shows that, if e were smaller, the term [(i + e)/e] • c exp(-,3Tg} would be increased 
without change of parameter C. 2) We may have underestimated the local mechanical dissipation, 
owing perhaps to a breakdown of the weak-shock theory. 3) Our assumption that all magnesium 
and calcium atoms are singly ionized in this region may be wrong. Both alternatives 2 and 3 
would require that the parameter C be increased. To verify this, see Equation 99. We consider 
these possibilities in hirn. 

Since the values of (n^,x^) = (l0“, 10'') , which we used to estimate the values of e for these 
lines, give slightly smaller values than would be obtained from any combination of (n^,x^) taken 
from the PSC model of Table 1, this would at least partially compensate for any errors in Equa- 
tion 161, used to obtain e. Also, if we look at Equation 170, which is our 



a. T. 
1 J a 


for Call for the first set of parameters (C,/3), we note that there is a monotonic outward decrease 
of the radiation field, which does not appear in the other solutions. Equations 171-173, until some 
distance out in the atmosphere has been reached. This lack of an increase in 



a. T. 

I J J 
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is not compatible with the observed emission 
core in the line profile of the H and K lines 
seen on disk spectra. Increasing e will not 
change this, but a sufficient increase in the ^ 
parameter C will, as illustrated by the solution. 
Equation 172, for these lines, where only C 
has changed among all the parameters. 

On these grounds, we reject the possibility 
that an error in e is responsible for our low 
Te (h) values obtained with the first set of 
parameters. 


Table 26 

Values of Local Mechanical Dissipation D(h) 
for the Eegion 500 < h < 1250 km. 


h (km) 

D(ergs cm ^ sec 

h 

D 

500 

0.030 

1000 

0.0024 

750 

0.010 

1100 

0.0013 

850 

0.0058 

1200 

0.00074 

950 

0.0033 

1250 

0.00055 


We must now consider which of our alternatives, 2) or 3), offers the more likely explanation for 
the increase in parameter C needed to raise the temperature. To determine this, it is useful to 
obtain the mechanical energy input D(h) and jD(h)dh for this region. D(h) from Equation 81 is 
given in Table 26 for the PSC model chromosphere in the region 500 < h < 1250 km. Table 26 
shows that jD(h)dh for the entire region yields approximately 5 x lO® ergs cm'^ sec“i. From this 
estimate and the contribution from the lower H“ region, we find that the total mechanical energy 
dissipated below 1250 km is only about 1.5 x 10 ^ ergs cm ~2 sec"^ The fact that our estimates of 
the total efflux of mechanical energy from the photosphere give more than an order of magnitude 
more mechanical energy than this value shows that there might be enough mechanical energy 
available to substantially increase the amount of chromospheric heating by shock dissipation. Also, 
the fact that the shock strength parameter 77 ^ becomes a significant fraction of imity in the region 
500 < h < 1250 km (which makes the weak-shock theory of doubtful validity) indicates that D(h) 
might be substantially larger than the initial values we have computed. Finally, we recall the 
arguments of Zirker in favor of single ionization for these metals in this region. From these 
considerations, we conclude that present evidence favors an increase in local mechanical-energy 
dissipation over a decrease in the ionization of the metals in this region. In any event, either or 
both of these effects Avill increase the parameter C, and this in turn will increase (h), as we 
readily see from t (2) and t ( 2 ) in Table 25. 

^Call ®MgII 


The manner in which (h) varies with h presents a further problem. For Call, (h) in- 
creases with h at least out to 1000 km, though the increase is slow. But for Mgll, (h) begins to 
drop somewhere between 750 and 1000 km. This is not a realistic temperatme distribution, at 
least not in the outer part of the region. Possibly our single e3q>onential representation of the 
mechanical dissipation in Equation 100 is inadequate, and an additional exponential term or terms 
may be required to account for an increase in mechanical dissipation in the smaller region. 

An increase in /3 would place the mechanical dissipation relatively farther out and change the shape 
of the temperature distribution. Such an increase in J3 for the Mgll doublet is reasonable, since it 
is doubtful that the assumption (NCB) = 1 is valid for this doublet in the lower part of the region 
in question owing to the large optical thickness, and beginning to decrease (NCB)nj at a higher 
point in the atmosphere has the effect of increasing /3 . 
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Also, we do not obtain the same distribution (h) from Mgll and Call with either set of 
parameters (C,/3). This means that our original conservation Equation 94 is not strictly satisfied. 
We could develop an iterative program to obtain a convergent temperature distribution in this two- 
line region. This should be done when we have more confidence in the mechanical-heating term D(h) 
and good estimates for (NCB)^^^ as a function of for the lines. Under present circumstances, 
this would be an imrewarding numerical exercise. 

What can we conclude about (h) in this region from these calculations? Obviously, we 
must set bounds to the temperature in the region 500 < h < 1250 km. We cannot find the tempera- 
ture gradient without a better knowledge of the (NCB)n,'s. As a lower limit, we would choose 
6000°, in the light of our H' study. As an upper limit, we would select 7500°, on the grounds that 
this value, calculated for the Mgll doublet, is the peak temperature reached anywhere in the 
region. On the other hand, the fact that our analysis does not adequately cover the outer, low 
region does not eliminate the possibility that may reach higher values. These limiting values, 
along with all our calculated values, are plotted in Figure 4 at the end of this section. 

A "plateau effect" appears in the distribution for Call, but since this depends ultimately on 
our choice of D(h), which is poorly known, this might be accidental. Such an effect is barely dis- 
cernable in the Mgll distribution of X^ (h). 


3. Temperature Distribution from Hydrogen Lya for h > 1250 km 

We do not extend the calculations for the region of the singly ionized metals beyond h = 1250 
km, because we have no reason to believe that the homogeneous low chromosphere extends beyond 
this height. Indeed, we have already noted the evidence of House (1961) which suggests that 
spicules may originate as low as 700 km, so some of our previous calculations for an assumed 
homogeneous region between 700 and 1250 km may be subject to modification on this account. We 
showed at the end of Section B of this chapter why we should expect Lya to become the dominant 
energy sink above the region of the singly ionized metals. We now proceed to calculate x^ (r^) 
from Lya. 

The results we now give for x^ from Lya pertain to heights above 1250 km. We do not 
attempt to relate to h because our knowledge of this upper chromospheric region is too poor. 
We shall refer briefly to the PSC model of the upper chromosphere for estimating parameters for 
the Lya temperature calculation. 

We start the Lya calculation with a study of the Lya source function. Using Equation 161 for 
e, with the mean Gaunt factor g given by 0.034 for x^ = lO"* and 0.124 for x^ = 10® from the van 
Regemorter (1962) study, we obtain for e; 


n 


e 


10 '° 


10 




1. 28 X 10"L 1. 28 X 10'® 


= lO" 


(174) 
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e 


T 


e 



1.48 X 10"^ 1.48 X 10"® . 


(1V5) 


We make the temperature calculation for both e ~ 10'® and 10"^. The value 10"® corresponds 
more closely to the upper chromospheric model of PSC. 

Evaluating parameter -q from Equation 33 gives tj = 4.3 x 10"®, where 6900°K, 

= 6000°K, and the dilution factor w^j equals 1/2. Thus, we conclude that, even for = IQn and 
e = 10 "®, 17 » e is a good approximation. 

So far, we have reproduced the results of Morton and Widing (1961). The approximation 

» 77 B* is more marginal for ovir purposes than in their case, because they overlooked the 
possibility that may be as low as 10^ at the base of the Lyo. region, and considered only 
> 2 X for which eB^ (t^) >> 778 * is liberally fulfilled. Using t_.j = 6900° again, we obtain 77 B* 
as 0.1574 • exp(-Yj)y/vj using Equations 33 and 34 and evaluating the atomic constants. This yields 
77 B* = 8.06 X 10"i®. The worst possible case for (t^) corresponds to = (l0‘^ 10“'), 

since lO^i , lO^ at the base of our Lya region. For these values, = 2.04 x 10"^^ 

Thus, at worst, eB„ (T^) 2 . 477 B*. Indeed, for = 1.2 x 10“*, we already have 14 . 0178 *. 

Therefore, we adopt the approximation eB^ (x^) » 178 * , and our results show that this is reasona- 
ble. The source function for Lya is indeed of the form given by Morton and Widing and discussed 
in Chapter II, namely. 


^ ( 35 ) 

When the expression for of Equation 35 and the form in Equation 100 are equated, we again ob- 
tain the equation for in the form of Equation 119, only now the approximation s = 77 is used for 
the sink term. 

Proceeding exactly as with Mgll and Can, we first obtain (C,/S) from Cj exp {-/3. x^} 

- > where "j” now refers to Lya. For Lya, = 0.059. We obtain 8 ^ from Equation 

81, where we assume that Dj = d in this upper region. To calculate d, we need the density scale 
height for h > 1250 km. We approximate this by the emission scale height for Paa of HeH at A4686, 
reported in PSC as 2.22 x 10® cm. We also need p(h = 1250 km), so that we can use 

-<h-12S0)/H, 

/ 0 (h) = p(1250) e ' ® , 83 = 2220 km , (176) 


to obtain the absolute value of p(h) at a given height. For p(1250), we use Equation 131 which 
estimates p(h) for the low chromosphere above 500 km. We obtain p(1250) = 8.0 x 10" gm cm"®. 
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Table 27 


Estimates of Neutral Hydrogen Densities for the Upper Chromosphere. 


h (km) 

Te 

(Table 4) 



nj^j (estimate) 

o 

11 

4) 

n =1010 

e 

1500 

23440 

9.5 X 10^ 

9.5 

10®-10'1 

10®-10’ 

2000 

36550 

43 

0.43 

lO'i-lO® 

10®-10i 

3000 

42130 

20 

0.20 

10® 

10 ^- 10 ’' 


Corresponding roughly to the two values of the shock strength = 0.5, 1, there are two different 
values for the mechanical energy input D, one an order of magnitude higher than the other. 

To derive n^, we solve the Saha equation for rij* ~ [n (neutral)]* . Then we estimate values 
for the departure coefficient bj in the upper chromosphere (using the Table 4 model of the upper 
chromosphere) and, from these, estimate The Saha equation for hydrogen can be written 


n 


1 


4. 11 X 10“‘® 


X 3/2 


157, OOO/T^ 
e ® 


(177) 


If we solve for n^*, we obtain the results tabulated in Table 27, assuming = n^ for the calcula- 
tion. There we also give our estimates of bj obtained from the early Thomas (1948b) and the 
previously mentioned Thomas and Pottasch (1959) studies of bj. From Table 27, we conclude that 
lies in the range 10® - 10^ in the upper chromosphere. 

There is considerable uncertainty in the values of the parameters (C,/3) to be used to represent 
the term ( 2/A.) D,/n, . for Lya. These uncertainties arise from the high approximate values we 
must use for the mass density p(h), the neutral hydrogen particle density and the shock 
strength parameter One way to assure a reasonable estimate is to demand that the mechanical 

Table 28 


Temperature Distribution as a Function of (Lycc). 


Tq (Ly-oc) 

n^ = 10 1 1 (cm 0) 

= 10 1° (cm ®) 

Cj = lO'ii 

11 

O 

o 

L 

Cj = 10"ii 

Cj = lO'io 

104.5 

9,865 

12,205 

12,205 

16,000 

10 3 

12,225 

16,035 

16,035 

23,368 

10^ 

11,995 

15,640 

15,640 

22,480 

10 

11,928 

15,527 

15,527 

22,151 

0 

11,851 

15,397 

15,397 

21,970 
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energy input corresponding to the chosen values of (C,/3) must equal the net radiative loss inferred 
from the observed Lya line profile. The net flux per c.p.s. and steradian is given by Jcexp {-^0 
between limits 0 andco. Multiplying by 4 v|Ai/| = 477 c|AAi /\2 ^ where izV\.| is the line width at half 
maximum, gives an estimate of the net flux in ergs cm” ^ sec”' . 

If we represent the expression c exp {-/Jt-q} for Lya by 


L'^ jLya( 1 ) 


10-11 


(178) 


L jLya(2) 


= 10”'o 


(179) 


we can estimate the net flux for each set (C,p) in the manner just prescribed. For (C,/S) = (10”", 
10”^") of Equation 178, we obtain 1.28 x 10 ® ergs cm”^ sec For (C,/3) = (10”'° , 10 ”^") of 
Equation 179, we obtain 1.28 X 10® ergs cm”^ sec”'. In both cases, we have used |AA.| = 10”® cm 
= lA. The Hinteregger (1964) observed flux in Lya of 3 X 10 ® ergs cm”^ sec”' is spanned by the 
above values. Since these values also agree in order of magnitude with estimates of ( 2 /A.)Dyn^. , 
based on the values of p(h), 0 ^, and 77 ^ already considered, we accept Equations 178 and 179 for 
the temperature calculation. 



Figure 4 — Graph of (h) obtained from H,CaII, (1,2), Mgll (1,2), 

and T,,(Tq) from Lya. 

The temperature calculation proceeds exactly as before for Mgll and Call. The results are 
given in Table 28 and exhibited graphically in Figure 4. A discussion of all resiilts, including 
these for Lya , is now in order. 




Chapter V 
CONCLUSIONS 


In the lowest region 0 £ h < 500 km, we have calculated the temperature distribution taking 
into account both the shock dissipation of acoustic waves and the departure of the continuum source 
function from the Planck function. The contribution of this latter, non-LTE effect is determined 
by the departure from imity of the coefficient b _ (h) , whose values we have calculated for the PSC- 
model chromosphere. 

Can the departure of b^_ from unity cause the temperature rise in the low chromosphere with- 
out any mechanical dissipation, i.e., in radiative equilibrium? This is the question raised by Cay- 
rel (1963), who concluded that it is possible in the low chromosphere. In the strictest sense this 
situation does not seem to be met, because of the incompatibility of radiative equilibrium and depar- 
tures from LTE in the H” ionization; this is demonstrated below. 

The statistical equilibrium condition for the bound level of H', expressed by Equation 141, can 
also be written in the form 


4t7 



By (Te) 


dv 



"h- °e,H- 






(180) 


The energy-balance equation for H~ ionization can be written 



where the bars over the brackets mean the enclosed quantities are suitably averaged over energies 
from (hvg) to Cayrel, in effect, equated the left-hand sides of these two equations to zero, which 



corresponds to radiative equilibrium, as Equation 181 shows. He then solved the resultii^ equa- 
tions for b^- and the temperature in the chromosphere corresponding to a suitable dilution factor 
for the photospheric background radiation field. The b^- enters the equations through b^-= ny/njj*-, 
after which rijj- cancels out of both of the equations Cayrel used. 

However, Frisch (1966) concluded that b„ =1 for the low chromospheric region under discus- 

2 

Sion. When b^ = 1, the only time Cayrel's equations are valid is for bj^_ = 1, as indicated by either 

or both of Equations 180 and 181. Thus, Cayrel's procedure, which yields b^,- > 1 ("^l.S), is invalid 

in the low chromosphere if b = 1. It does not seem possible {is not possible if b„ = l)for depar- 

2 2 
tures from LTE in the H“ ionization to exist in the low chromosphere in radiative equilibrium. 

Therefore, the temperature rise cannot be produced solely by non-LTE ionization of H" in radia- 
tive equilibrium. 

It is interesting to consider the relative contribution of the radiation field and mechanical- 
energy input in determining the temperature distribution in the light of these conclusions. Recalling 
that we have used the symbol D, without subscripts, to represent the mechanical energy input in ergs 
cm"^ sec"*, we can write 


D = 



"h- (' " \) 






(182) 


for the mechanical input, where b„^ = 1 has been assumed. It follows directly from Equation 182 
that, as b„- increases ( bjj- > 1), D increases, and vice versa. So in that part of the low chromo- 
sphere where H" dominates the continuum opacity, the non-LTE effect corresponding to a depar- 
ture of b„- from unity and the mechanical heating of the gas go hand in hand. This is exactly what 
we would expect on physical grounds. We can envision a cycle of collisional and radiative ioniza- 
tion (with the associative detachment reaction H + H" - Hj + e dominating the collisional ioniza- 
tion) followed by collisional and radiative recombination, where the collisional ionizations minus 
the collisional recombinations, representing the increase in mechanical energy, are just balanced 
by the excess radiative emission, obtained from the radiative recombinations minus the radiative 
ionizations. 


A comment on the Frisch (1967) model of the low chromosphere is now in order. We have al- 
ready noted her claim that the Paschen continuum opacity predominates over the H" bound-free 
opacity for h ~ 400 km. As she assumes bj = 1 for the Paschen continuum, the departure of the 
continuum source function from the Planck function is less under these conditions than in our 
calculations. 


We have estimated the ratio Paschen continuum opacity to the H" bound-free opac- 

ity in Table 22, for the PSC model chromosphere, and the results show that this ratio reaches a 
maximum of 0.119 at h = 500 km. Why is there such a discrepancy between Frisch's and our results? 

It is not easy to answer this question, but we can make the following observations. Frisch states 
that, above 400 km, the Paschen continuum opacity predominates, but she does not say how she 
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arrives at this. She arrives at her final models by solving a system of rate-process and energy- 
balance equations plus an equation of state. If we use the (P^jT^) values from her final models, 
compute n^(p^,T^) and take Equation 164 to compute R3 we derive for her models I, II, and HI, 
respectively: (n^ = 3.98 x 10^°, Rg = 0.0045), (n^ = 2.67 x 10“ , R3, fb = 0.142), and (n^ = 9.6 
X 10^° , R3 fb = 0.269) at those heights nearest to 400 km among her tabulated values. Ignoring the 
first set (which is for her zero mechanical dissipation case), we conclude that her statement on the 
Paschen continuum opacity still does not seem consistent with her results, which predict conditions 
similar to those of the PSC model. This conclusion is true even though her (h) lie below the (h) 
of PSC in Table 1, and we see from Equation 164 that Rg fb varies inversely with n^. Her results 
and Equation 164 support the conclusion that the H" bound-free opacity dominates the lower region 
0 £ h £ 500 km. Thus, we retain the conclusion that the H" bound-free opacity is dominant. 

The reasonable conformity between our calculated temperature distribution and that of the PSC 
model of Table 1 tends to confirm the conclusion that rises at least 1000° from h = 0 to h = 500 
km. This suggests a lower limit of - 6000°K at h = 500 km, a value lying somewhat below our cal- 
culated value and that of Table 1 (6150° and 6285°, respectively). We should note, however, that we 
have established only that the PSC model offers a reasonable description of the H" region in view of 
our results. We have not, in any sense, proved that the PSC model is correct. Nor have we consid- 
ered the possibility that the type of consistency test we applied to the H~ region of the PSC model 
might not apply equally well to another model. So the possibility that the temperature at h = 500 km 
may lie below 6000°K does exist, though this higher value seems more probable. 

The results of the temperature calculation in the next higher region above 500 km present sev- 
eral problems. The numerical problem of convergence of the two temperature profiles obtained 
from MgH and CaH is easy to solve, with convergence assured in the two-line case, but the solution 
would offer no useful information until a more detailed study of the (NCB)^^^ 's has been completed. 

So we shall consider this problem no further here. This leaves the following problems: 

I. In that portion of the chromospheric gas found in the region 500 < h £ 1250 km (and as long 
as we retain the assumption of spherical symmetry there, this means all the gas in the region), we 
concluded that either 1) an order of magitude more energy (-^5 x 10® ergs cm"^ sec'M is needed 
to satisfy (min) ~6000°K at h = 500 km for the two lines treated than is provided by rough esti- 
mates of the available mechanical energy given by the weak- shock theory, or 2) the magnesium and 
calcium are predominately neutral in this region, in conflict with Zirker's conclusions reported in 
Section B of this chapter. 

II. The net radiative loss from the MgH doublet corresponding to the first set of the heating 
parameters (C,/S) agrees to within a factor of 2 with Osterbrock's (1961) estimate of 2.5 x 10® ergs 
cm"^ sec"^, thus placing it about an order of magnitude above the observationally based Osterbrock 
value for the second set of (C,/3), which represent more heating. 

HI. The net radiative loss from the CaH doublet is significantly larger in calculated than in 
observed values, even for the lower values of the heating parameters (C,/3). The flux in H and K is 
estimated in the Osterbrock study as ^10'' ergs cm"^ sec"^, while the lower values for (c,/3) cor- 
respond to > 10® ergs cm" 2 sec. 
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IV. The temperature gradients obtained from the line calculations are wrong— at least in the 
outer part of the region. 

V. Possible departures from spherical symmetry have so far been ignored. Fredga's (1966) 
spectroheliograms of the Sun in Mgll \2802.7 confirm the idea that most of this emission comes 
from essentially the same regions as the CaKj emission. Thus, severe departures from spherical 
symmetry, apparently associated with the magnetic-field structure, are virtually assured. 

Let us consider these problems in turn. We have already mentioned the possibility that the 
weak-shock theory, valid only for values of the shock- strength parameter 77^ < 1 (how much less 
is difficult to establish), might break down for values of approaching unity in this region. H 
this happens, we have underestimated the mechanical heating, and the correct value would be in the 
right direction to produce the higher temperatures. 

Zirker's (1956) conclusion that the metals should be singly ionized in this region seems well 
founded, the chief source of error in his study being the data from the 1952 eclipse. Of the two 
alternatives, 1) mechanical dissipation in the regions producing emission from the singly ionized 
metals significantly exceeds 10® ergs cm"^ sec"' or, 2) those metals are predominately neutral 
there, the former seems more probable. 

Consider the problems posed by the apparent discrepancy between the calculated and observed 
net radiative losses from Mgll and Call. Let us examine the procedure for estimating net radiative 
losses from the observed emission cores, as outlined by Clearman (1953). It is this work to which 
Osterbrock and Athay ultimately refer in discussing Mgn. The procedure is to estimate an equiva- 
lent width for the emission core as high as the adjacent continuum and multiply this by the black- 
body emission at a temperature characteristics of this continuum. Clearman used =5000°. If he 
had used =5800°, this would have increased the predicted flux by more than 4 times in the Mgll 
resonance doublet. So in the Mgll case, at least, the flux estimate is sufficiently uncertain to raise 
doubts about the reliability of the quoted numbers, though it is unlikely that the observations could 
give a flux value of greater than 10® ergs cm"^ sec"' for the average value over the solar surface. 

In the Call case, no errors of the type just described could possibly explain the vast discrepancy 
between the calculated and observed fluxes. Here we must examine the possibility that the region 
producing the H and K emission may be optically thin compared with our estimates (Zirker, 1968), 
or that our two-level atom model has failed. 

We could at least partially remove these discrepancies by forcing (NCB)^^^^ to approach zero 
higher in the atmosphere, corresponding to a more realistic treatment of the relevant departure 
coefficients . Such an approach, which would force the mechanical heating farther out for the 
line in question (requiring a change in the mechanical-energy input distribution or a different dis- 
tribution of the other mechanical-energy sinks), would also raise the temperature in the outer layers. 
Thus, we might also be able to reduce or change the unrealistic, negative temperature gradient in 
those layers. 

Finally, the probable departures from spherical symmetry indicate a higher local value for 
the mechanical heating, where the Mgll and Call emission is strong. If 10 percent of the gas in 
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this region emits most of the observed flux, the higher temperatures obtained from the second set 
of heating parameters (c,/S) are obtained there, and the total mechanical heating will not signifi- 
cantly exceed the quoted fluxes from Mgll, though the Call emission will still be too large to agree 
with observations. Unfortunately, the behavior of the other 90 percent of the gas is still unknown. 

To conclude the discussion of the Mgll and Can doublets and the temperature distribution cal- 
culated from them, let us ask: what model for the region 500 < h < 1250 km is roughly consistent 
with the results of our calculations? The preceding remarks demand a two- component model for 
this region, with one component (producing most of the emission) covering about 10 percent of the 
surface. In this component, shock dissipation, greatly influenced by the stronger magnetic fields 
in the spicules (which we identify with this component), could produce the necessary mechanical 
energy to balance the observed emission. The temperature in this component should be 7000° 

± 1000°K. Nothing can be said about conditions in the other component— a serious problem. 

Several changes occur in the chromosphere where the transition to the Lya emitting region 
takes place; 1) The temperature will begin to rise rapidly above the 6000° - 8000°K values charac- 
teristic of the singly ionized metals region. Such a rapid rise begins in the PSC-model low chro- 
mosphere around 1250 km. 2) Associated with this rapid rise, (Lya) will decrease rapidly owing 
to hydrogen ionization. Thus, the condition of detailed balance in Lya, valid in the singly ionized 
metals region, becomes invalid, and net radiative losses in Lya become large. 3) The dominance of 
Lya is furthered by second-stage ionization of the metals Mgll and Call when exceeds 10,000°. 

4) Helium emission should be observed above 1250 km. Significant helium emission in Hel and 
\4686 of Hen (whose resonance line \304 cannot be observed below the ozone layer of our atmosphere) 
are observed at and below 1250 km, as shown by the Khartoum data reported in PSC. 5) There is 
direct evidence, already discussed, for departures from spherical symmetry above 1500 km, in addi- 
tion to the strong possibility of a departure below this height. 

The results of the Lya calculation support this picture, even to the extent of fully justifying the 
neglect of the contribution of helium emission to the energy balance in the region producing strong 
Lya emission. We recall that a temperature of 50,000°K was required (Zirker, 1960) for helium 
emission to predominate. Even the effect of strong departures from spherical symmetry in the Lya 
emission (which could concentrate the mechanical heating of our spherically symmetric study by 
more than an order of magnitude) could not raise the temperature to this value, in view of the re- 
sult showing an increase of ~6000° corresponding to an order-of -magnitude increase in mechanical 
heating. The greatest uncertainty in the Lya calculation, other than that caused by ignoring depar- 
tures from spherical symmetry, is the estimate of the mechanical heating. The final temperature 
values should be considered as "ballpark" figures only, offered to illustrate the effects of varying 
the mechanical heating and the electron densities. 

What is the stability of the temperature gradients obtained from these calculations? We have 
already noted that the negative dT^ /dh found at small for all the calculations above h = 500 km 
(quite pronounced for Mgll) conflicts with the expectation that, when one radiator becomes an in- 
effective loss mechanism, the temperature should rise to a higher value at which a new radiator 
can become effective, provided that we maintain a mechanical input per unit mass of material that 
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is approximately constant. The quantity D(h)/p (h) actually increases slightly in the ionized metal's 
region, so we certainly do not believe in the negative dT^/dh. 

Ignoring this aspect of the results and considering Mgll and Call to represent a single tempera- 
ture plateau (or perhaps two temperature plateaus), we ask if the results for these ionized metals 
and the other radiators are reasonable as regards general stability. The general -stability criterion, 
given in Chapter n, Section D, is 


dE„ , , dE . 

out in 

> 


dT. 


dT. 


(45) 


where it is convenient to express E in units of ergs gm i see”* . 

Starting the stability study with the H" region, 0 < h < 500 km, we find that E.^ =D(h)/p(h) in- 
creases rapidly with height h and hence with increasing , as the shock strength parameter in- 
creases rapidly in this region and D(h)a [77^ (h)]^. What can we say about E^^^ ? For the free-bound 
continuum transitions of the H" ion. 


n^Oi A,l (NRB),l (hv),Ldk/p(h) , (183) 

in the notation of Chapter n. Section C. We could numerically solve equation 183 for e^^^ directly, 
or we could use the estimates of the excess of emission over absorption in H" from Table 14. The 
latter procedure is more direct. Noting that the results of Table 14 are expressed in units of ergs 
cm"^ sec" ‘ , requiring that we divide them by p (h) to get the net radiative loss per unit mass, and 
recalling that, for 0 < h < 200 km, the values of Table 14 are greatly overestimated, we see that 
dE^^j/dT^ > 0 is certainly satisfied for the assumed density-scale height of 130 km. Since the big 
jump in 77 ^ (hence D(h)/p(h) = E.^) occurs for the range in h where there are many uncertainties in 
all the estimates, the possibility that condition 45 will be violated there is not too disturbing. Be- 
yond h 200 km, dEj^/dT^ does not increase rapidly and with dE^^^ /dT^ > 0, condition 45 can prob- 
ably be fulfilled. Since the estimate of the variation of with h is rough (as is the estimate of 
excess emission) it would hardly be worth a more quantitative study at this point. 

Then there is the question of the Mgll and Call region. Again, we shall not attempt to calculate 
exactly where we expect this region to begin. We simply note that beyond h = 500 km, the Paschen- 
continuum opacity becomes increasingly important. This means that at some height above 500 km, 
according to our estimates, H” will cease to dominate the continuum source function and the tem- 
perature determined thereby. Also, at some such height, the (NCB)^^^^,^ for Call and Mgll will be- 
gin to depart from unity. 

If the region 500 < h < 1250 km is a region of singly ionized magnesium and calcium, then the 
form of E^^j we must use is 
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(184) 


Since and (hi^)uL constants for each line, the variation of over 500 < h < 1250 km is 
determined by the variation of , p , and (NRB) Qualitatively, we may state that (NRB)jjj^ in- 
creases from zero to unity through the region, and that decreases faster with increasing height 
h than as (NCB)^ (1 - b^/bj^) increases. Thus, since andp(h) have the same scale height if 
the metals are singly ionized, the ratio riy/p(h) decreases with height, somewhat offsetting the ef- 
fect of the increasing (NRB)^^. However, does not increase rapidly through this region. From 
the variation in over the region, E.^ = D(h)/p(h) only increases by about a factor of 2. Thus, only 
a small increase inE^^^ is necessary to satisfy condition 45. Again, without making a quantitative 
study, we note that condition 45 is probably satisfied by MgH and Call. 

Could Lya. possibly play a significant role at some height below 1250 km, where we expect it 
to predominate? We argued that Lya should be negligible below 1250 km on the grounds that the 
optical depth Tg(Lya) is greater than 10= below this height, according to the PSC-model chromo- 
sphere of Table 1. For such optical depths, we noted that detailed balance is expected for Lya. Un- 
less this conclusion is in error, or unless the PSC model for (Lya) is in error so that the upper 
part of the region below 1250 km loses Lya radiation, we can stand by the conclusion that Lya losses 
and the associated higher temperatures are confined to the region above 1250 km. 

Finally: is the calculated Lya region thermally stable? According to Table 27, where (h) is 
given for three heights (in a highly questionable upper-chromosphere model), n^(h) appears to be 
almost constant with height h after the bulk of the hydrogen has been ionized. With p(h) almost con- 
stant in the Lya region (scale height chosen = 2200 km), and (NCB)j^„ = 1 (a good approximation for 
the outer part where hydrogen is largely ionized), the variation in E^^^j from Equation 184 is deter- 
mined largely by the variation in n^. To obtain this, knowing ~ , would require a more de- 

tailed knowledge of the variation of b, and b^ through the Lya region than we have. We merely note 
that, since here E^^^ is roughly constant with height (y^ is assumed constant.), it may be possible 
for dE^^ydX^ > 0 to be fulfilled, thus satisfying condition 45. 

It is too early to appraise the method offered here for calculating the chromospheric tem- 
perature distribution, as compared with the standard approach, namely: assvuning a distribution 
of electron densities and temperatures for various chromospheric models and computing line pro- 
files from the resulting source functions as outlined in NETPRF, after which a comparison of the 
calculated with observed profiles is made to select the best model. Our approach has been to 
make a reasonable guess of the mechanical energy input (which still depends on temperature), and 
construct a source function from this to solve the transfer equation. From the solution a tempera- 
ture distribution can be calculated. 

The former, standard approach has recently been applied to the Mgll and Call doublets by 
Athay and Skumanich (1968). Since their results are also based on a two-level atom, but were ob- 
tained after a much more thorough series of numerical calculations, it will be interesting to com- 
pare their results with those we hope to obtain from a more detailed calculation. In addition, the 
results of the new calculation should yield profiles good enough to compare with observed profiles. 
Only then will it be possible to offer a reasonable estimate of the utility of calculating the tempera- 
ture as described in this report. 
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Appendix A 

The General Line Source Function 



In this appendix we derive, in detail, the line source function of Equation 25, Chapter II, Sec- 
tion C. We begin by applying the very general condition of a statistically steady state expressed 
by Equation 16 


dn^ 

^ = 0, (16), (A- 1) 


where, of course, is the number density of the k-th energy level. (We shall use a double refer- 
ence number for equations that appear in the body of the report.) From Chapter II, Section C, we 
also have the following completely general expression for the line source function S^, which is 


S. = ^ • B, (T„) 


(18), (A-2) 


where and are the profile coefficients for spontaneous emission and absorption, respectively, 
and 


B. (T, J 


2hu3 \-i 

c2 7 


(19), (A-3) 


with the same notation used throughout the report. Since, for reasons given in Chapter II, Section 
C, the profile coefficients j „ and vary in a similar way across the Doppler core of resonance 
lines, we restrict our treatment here to determining the detailed form for (Te*)- 1''°^ ^ detailed 
treatment of the profile coefficients, see Chapter III, Part 1 of NETPRF or the first paper in the 
Source Function Series (Thomas, 1957). These references also outline the derivation of a general 
expression for (T^^), so this appendix is intended to fill in some details for readers unversed in 
non-LTE spectral line theory. 


The object of this derivation is to find a general expression for the ratio in Equation 

A-3. This expression must include the general collisional and radiation field dependence of n^/n^. 
We obtain such a ratio by applying the SSS condition A-1 to both levels (u, L) of the spectral line. 
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This gives 
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for the upper level u, and 
22 + % V(NRB)ul + ••• 

^>L 


+ 22 ^ C.^(NCB);, + 

j <L 


I A^L^i "e (NRB)^j^d/< 
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22 (NRB)Li " 22 (NCB),^, 

1>L 

_ '^i'U _ 


(A- 5) 


for the lower level L. All quantities appearing in Equations A-4 and A-5 are defined in Chapter n, 
Section C. In particular, the net radiative brackets (NRB) are given by Equation 39, the net colli- 
sional brackets (NCB) by Equation 44, and (any L and U where L < u ) by Equation 23. Note that 
Equations A-4 and A-5 are equivalent to Equation 24 where k becomes equal to u and L, respectively 
We retain the convention of including collisional processes that have an upper level in the continuum 
in the middle term on the right-hand sides of Equations A-4 and A-5. 


It is useful to introduce several new quantities at this point. Though this appears to compli- 
cate matters at first, these quantities have physical significance, and their introduction greatly 
facilitates our derivation and the subsequent physical interpretation. First, we replace the 
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radiation terms involving the continuum by 
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A^„n, (NRB)^jjdK - 
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(A- 6 ) 


(A-7) 


noting that is simply used to represent p(rad; C, k) - R(rad; k, c) in the notation of 

Equation 21. Next, we introduce four quantities which couple U and L to levels other than L and u, 
respectively. These are 


^ + Jlc ’ (A- 8 ) 

•t>L 

■l!‘U 

1>L 


S 


3 


22 Auj (NRB)^jj 

j <u 
j/L 


+ 22 (NCB)ij.e + Juc 

-£>u 


S 


4 


22"iCju(NCB).„ 

i<u 

j?^L 


+ 22 , 

-t>u 


(A- 10) 


(A- 11) 


where and S 3 represent terms that depopulate levels L and u, respectively, and S^ and repre- 
sent terms that populate L and U, respectively. These "indirect" processes (not coupling U and L 
directly) represent a complete set of all indirect processes except for Rj,„ and which we in- 
troduced in equations A -6 and A-7. To include these also, we define 


= + S 4 


(A- 12) 


^CL ^ 


(A- 13) 
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- <!p- 


folloAving the notation of NETPRF. Now we have included all indirect processes in and 

. We introduce no more new quantities. 

If we combine Equations A-4, A-6, A-10, A-11, and A-12, we can write 


[^3 ^ ^UL ^LU )lU ' 


(A- 14) 


In a like manner, combining Equations A-5, A-7, A-8, A-9, and A- 13 gives 


Fl = "l [Si + C^u(NCB)lu] - "uA„l(NRB)ul • 


(A- 15) 


We express A- 14 and A- 15 in a slightly different form by noting that 


Aul 


"uAul "u®UL 


h 


1 


J 0 dv 

J V 


(A- 16) 


for '/'y = , which we assume throughout, and 

"lClu(NCB)lu = "l^lu “ "u^UL • (A-17) 

Equation A-16 follows from Equations 21 and 36. Equation A-17 follows from Equations 22 and 40. 
Substituting Equations A-16 and A-17 in Equations A- 14 and A-15 gives the equations to solve for 
and : 



” 

- 


“ 


"u 

1 


IT 4> dv + C„, 

•J V UL 

" "l 


T di^ + C 

^LU 


(A- 18) 



AuL + Fyy 

1“ 

It 0 dv + c„, 

•J V UL 

+ "l 


J. dU + Cyy 


_ 

- 


L 

- 


(A- 19) 


These two equations express only direct processes between levels u and L except for the terms S^, 
s^, Fy, and Fy . In addition, the net brackets for the direct processes no longer appear. We solve 
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for and by determinants. This gives, for the ratio nj^/n^ ; 


, + + A„l(Fj. + F„) + C„^(F,+F„) 

Blu j h dz. (F„ + Fj + (Fu + F J + S ^ F„ 


+ S3F, 


(A-20) 


We obtain a long expression for by substituting Equation A-20 in Equation A-3, Now recall 

the relations among the three Einstein coefficients , and expressed by Equations 38 and 

reproduced here; 


'uL 2hv^ 


SlBlu = guB„ 


(38), (A-21) 


Also recall that, for a Maxwellian electron-velocity distribution, we have shown that 


C = C — 

^UU '-UL 


Xlu 


hi^ 

kT. 


(43), (A-22) 


Using Equations A-21 and A-22 in the expression obtained for B^ (T„) by substituting Equation A-3, 
we arrive at 


(1 B.(tJ 


B (T 

V \ e 


2h!>3 gj_ 5j Fj. 


2 s A F "f" K 

&u "UL 


1 + 


Hi (i-e-^Lu) 

/S3 F^ g^ S F„ 
\Aul Fl + F„ - g„ F^+F„ 




We need only consider the third term in the denominator. 
First, note that 


^3 F^ 
Aul Fi.+F„ 


Fl gu ^ /gL ^ F„ \ 
Fu Sl \®u A(jl FL'^'F^y 


(A-23) 


(A-24) 
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where the term in parentheses is identical, except for a sign change, to the second part of the 
third term in the denominator. Thus, 


/ S 3 ^ gi, ^ \ ^ 

V^UL ^UL 

Sl Si F„ ( f ^ g„ 83 \ 

Su ^UL ^ V U ^ 1 / 


(A-25) 


Now consider the defined quantities we introduced in Equations 27 - 29 in Chapter n. Section 
C. We reproduce them here; 




(27), (A-26) 


gL Si /F, g, 83 

^ gu ■^UL V^U ®L S] ) ’ 


(28), (A-27) 


2hv3 Sl Sj F„ 

VB - 2eA F+F 

c-' &U UL u 


(29), (A-28) 


Using these definitions and Equation A-25, we can write Equation A-23 for (T^,j) in the desired 
form: 


jj. + eB^ (tJ + 7,B* 
{”^ex) ” 1 + 6+77 


(25), (A-29) 


The derivation is completed. 

A brief discussion of the physical significance of the terms appearing in Equation A-29 has 
been given in Chapter n. Section C. For a more complete discussion of the physical interpretation 
of Equation 29, see Chapter in of NETPRF. 
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Appendix B 

Special Solutions to the Transfer Equation 


In Chapter HI, where we derived the method for making the calculations, we deferred to this 
appendix the details of deriving and solving the basic differential equation of the problem, Equa- 
tion 105. We do this now, starting with the transfer equation in the form 


dr 

V 


Operating on both sides of Equation B-1 with / • • • dfi, we obtain 

Att Jv ’ 


(6), (B-1) 


(B-2) 


provided that is not a function of m = cos0 , Now we introduce a new auxiliary quantity K„ , de- 
fined by 


K 


V 



I,. (m)M^ d/i . 


(B-3) 


With this definition. Equation B-3, we can operate on Equation B-1 with J/j. ■ ■ ■ 6/j.to obtain 


^ _ J_ 

47T 


(B-4) 


where, again, we assume that is not a function of m. 

Now we introduce the Eddington approximation we use here, by replacing the in the integral 
of Equation B-3 by its mean value, namely 1/3. This immediately gives 

- 3K^ . (B-5) 

This approximation ignores the anisotropy of the radiation field near the boundary where - 0. 
Since we have applied the subsequent equation to strong spectral lines whose behavior a.t Z 1 is 
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our primary concern, this approximation is adequate for the calculations we have performed. 
Combining Equations B-2, B-4, and B-5, we obtain 


1 

3 


= J. - s. . 


withdr^ = 3'^^^ xJdTg, Equation B-6 becomes 


^^0 


(Jv-S.) . 


(B-6) 


(102), (B-7) 


which is the equation that we solve. The form of the source function is given in Chapter III by 
combining the general source function for the line and continumn of Equation 8 with detailed ex- 
pressions for the individual line and continuum source functions of Equations 100 and 101, respec- 
tively. For reference, we write 


Sc 

1 + 


(8), (B-8) 




1 


J 0 dv + Ce 

J If ^ 




(100), (B-9) 


S. = Sj(l+Ae"""°) 


(101), (B-10) 


where = dr^/dr^ and we assiune that only a single line is important at frequency If we define 
Sj, = (1 + r^)"‘, as in Equation 104, and replace all frequency points by quadrature points "i", 
we obtain 




J.-S. a. J. - S. - S. r, S, (l + Ae 

1/ 1 lll\ / 


4_. 

) 


(105), (B-11) 


for an n -point quadrature. In Chapter III, we gave physical arguments for dropping the last term 
S . r . Sj Ae This gives the equation we solve: 


dM, 


dr 


2 

0 


J.-S. y~' a. J . - 5 . Ce - S . r . S, 

J 1 1 / J J 1 111 


(109), (B-12) 
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To solve Equation B-12, a second-order linear differential equation, we obtain the general 
solution to the homogeneous equation and add a particular solution to the complete inhomogeneous 
equation, thus giving the general solution to the inhomogeneous equation. We then apply the bound- 
ary conditions of our problem to obtain the solution used for our calculations. The homogeneous 
equation corresponding to Equation B-12 is 




dr 


2 

0 



(B-13) 


If we assume a solution of the form 


where 


and 


= s. 


Bii 




8 . 






(B-14) 


(B-15) 


(B-16) 


and substitute Equation B-14 in Equation B-13, we obtain the following equation for the k.'s: 



(111), (B-17) 


Now we need a particular solution to the inhomogeneous equation. Assume for this that 


Ji, 



(B-18) 


Substituting Equation B-18 into Equation B-12 yields the following two relations; 


n 


P - S, 


= 0 


(B-19) 


’A 
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(B-20) 


n 

S^2^a.Q.+S^C 



The solution to C- 19 is P; =1; this is evident if we recall that 5 j = (1 + r . )"i and 


n 



We can obtain an explicit expression for Qj from Equation B-20. Summing a. q. over "i", and noting 
that the summation index is arbitrary, we can write 



(B-21) 


Finally, substituting Equation B-21 in Equation B-20, we obtain 


Qi = 



(112), (B-22) 


Thus, the general solution to the inhomogeneous Equation C-12 is the sum of and J^p , and is 
given by 


Ji 





(B-23) 


The +L. and _L. are determined by the boundary conditions. 

So far, we have not introduced any restriction on the extent of the atmosphere. We do this now 
by recalling that, for all the spectral lines for which we actually solve Equation B-12, we showed 
that the chromosphere could be treated as a "semi- infinite" atmosphere. The boundary conditions 
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for such an atmosphere can be stated: Jj ^ Sj as tq"”- Thus _Lj = 0 in Equation B-23. The 
are then obtained from the boundary condition at t,, = 0. Krook (1955) has provided the outer bound- 
ary condition we use. This is 


Ji 


d 


dJi 


(B-24) 


Applying condition Equation B-24, to the solution, Equation B-23, with _Lj = 0 gives for the +L. : 




+ S.-i 


Qi 




= 0 


(113), (B-25) 


The solution of Equation B-12 becomes 


Ji 


1+8 






;) 


(109), (B-26) 


for a semi-infinite atmosphere, where everything has now been determined. 
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